OFFSET
1,3
COMMENTS
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
The asymptotic density of the occurrences of 0's in this sequence is 1 - Product_{k>=1} (1 - 1/prime(k)^(A000005(k)+1)) = 0.28458561874647424574... . - Amiram Eldar, Nov 11 2025
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..10000
FORMULA
Multiplicative with a(prime(k)^e) = e! * binomial(A000005(k), e). - Andrew Howroyd, Nov 10 2025
EXAMPLE
The prime indices of 90 are {1,2,2,3}, with choices:
((1),(2),(1,1),(3))
((1),(1,1),(2),(3))
((1),(2),(1,1),(1,1,1))
((1),(1,1),(2),(1,1,1))
so a(90) = 4.
MATHEMATICA
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Table[Length[Select[Tuples[Select[IntegerPartitions[#], SameQ@@#&]&/@prix[n]], UnsameQ@@#&]], {n, 100}]
PROG
(PARI) a(n) = { my(f=factor(n)); prod(i=1, #f~, my([p, e]=f[i, ]); e!*binomial(numdiv(primepi(p)), e)) } \\ Andrew Howroyd, Nov 10 2025
CROSSREFS
For divisors instead of constant partitions we have A355739 (also the disjoint case).
For prime factors instead of constant partitions we have A387136.
For initial intervals instead of partitions we have A387111.
For strict instead of constant partitions we have A387115.
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
