OFFSET
1,12
COMMENTS
This is an auxiliary function for computing A317145 with help of A064988. Note the similarity of the formula to that of A300385, with only difference being in the value of a(1) and that here we have multiplication (*) instead of addition (+) between primepi(p) and primepi(q).
From Gus Wiseman, Oct 09 2018: (Start)
Combinatorial interpretation is: In the poset of multiset partitions ordered by refinement, number of maximal chains from the n-th multiset multisystem (A302242) to the maximal multiset partition of the same multiset, assuming n is odd. For example, the a(315) = 10 maximal chains are
{{1},{1},{2},{1,1}} < {{1},{1},{1,1,2}} < {{1},{1,1,1,2}} < {{1,1,1,1,2}}
{{1},{1},{2},{1,1}} < {{1},{1},{1,1,2}} < {{1,1},{1,1,2}} < {{1,1,1,1,2}}
{{1},{1},{2},{1,1}} < {{1},{2},{1,1,1}} < {{1},{1,1,1,2}} < {{1,1,1,1,2}}
{{1},{1},{2},{1,1}} < {{1},{2},{1,1,1}} < {{2},{1,1,1,1}} < {{1,1,1,1,2}}
{{1},{1},{2},{1,1}} < {{1},{2},{1,1,1}} < {{1,2},{1,1,1}} < {{1,1,1,1,2}}
{{1},{1},{2},{1,1}} < {{1},{1,1},{1,2}} < {{1},{1,1,1,2}} < {{1,1,1,1,2}}
{{1},{1},{2},{1,1}} < {{1},{1,1},{1,2}} < {{1,1},{1,1,2}} < {{1,1,1,1,2}}
{{1},{1},{2},{1,1}} < {{1},{1,1},{1,2}} < {{1,2},{1,1,1}} < {{1,1,1,1,2}}
{{1},{1},{2},{1,1}} < {{2},{1,1},{1,1}} < {{2},{1,1,1,1}} < {{1,1,1,1,2}}
{{1},{1},{2},{1,1}} < {{2},{1,1},{1,1}} < {{1,1},{1,1,2}} < {{1,1,1,1,2}}.
(End)
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..12960
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000
FORMULA
If A001222(n) <= 2 [when n is one, a prime or semiprime], a(n) = 1, otherwise, a(n) = Sum_{p|n} Sum_{q|n, q>=(p+[p^2 does not divide n])} a(prime(primepi(p)*primepi(q)) * (n/(p*q))), where p ranges over all distinct primes dividing n, and q also ranges over primes dividing n, but with condition that q > p if p is a unitary prime factor of n, otherwise q >= p. Here primepi = A000720.
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 08 2018
STATUS
approved