OFFSET
1,2
COMMENTS
The asymptotic density of the occurrences of 0's in this sequence is 1 - Product_{p prime} (1 - 1/p^(A000041(p)+1)) = 0.13580468148150748566... . - Amiram Eldar, Nov 11 2025
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..10000
FORMULA
Multiplicative with a(p^e) = e! * binomial(A000041(p), e). - Andrew Howroyd, Nov 10 2025
EXAMPLE
The prime factors of 9 are (3,3), and the a(9) = 6 choices are:
((3),(2,1))
((3),(1,1,1))
((2,1),(3))
((2,1),(1,1,1))
((1,1,1),(3))
((1,1,1),(2,1))
MATHEMATICA
Table[Length[Select[Tuples[IntegerPartitions/@Flatten[ConstantArray@@@FactorInteger[n]]], UnsameQ@@#&]], {n, 30}]
PROG
(PARI) a(n) = { my(f=factor(n)); prod(i=1, #f~, my([p, e]=f[i, ]); e!*binomial(numbpart(p), e)) } \\ Andrew Howroyd, Nov 10 2025
CROSSREFS
Twice partitions of this type are counted by A296122.
For strict partitions and prime indices we have A387115.
For constant partitions and prime indices we have A387120.
A003963 multiplies together prime indices.
A120383 lists numbers divisible by all of their prime indices.
A289509 lists numbers with relatively prime prime indices.
KEYWORD
nonn,mult
AUTHOR
Gus Wiseman, Aug 26 2025
STATUS
approved
