OFFSET
1,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..10000
FORMULA
Totally multiplicative with a(p) = A000041(p) for prime p. - Andrew Howroyd, Nov 10 2025
Sum_{n>=1} 1/a(n) = 1/Product_{n>=1} (1-1/A058698(n)) = 3.88219666056645519287... . - Amiram Eldar, Nov 12 2025
EXAMPLE
The a(1) = 1 through a(7) = 15 ways:
(1) (2) (3) (2)(2) (5) (2)(3) (7)
(11) (21) (11)(2) (32) (11)(3) (43)
(111) (2)(11) (41) (2)(21) (52)
(11)(11) (221) (11)(21) (61)
(311) (2)(111) (322)
(2111) (11)(111) (331)
(11111) (421)
(511)
(2221)
(3211)
(4111)
(22111)
(31111)
(211111)
(1111111)
MAPLE
a:= n-> mul(combinat[numbpart](i[1])^i[2], i=ifactors(n)[2]):
seq(a(n), n=1..56); # Alois P. Heinz, Nov 10 2025
MATHEMATICA
Table[Length[Tuples[IntegerPartitions/@Flatten[ConstantArray@@@FactorInteger[n]]]], {n, 30}]
PROG
(PARI) a(n) = { my(f=factor(n)); prod(i=1, #f~, my([p, e]=f[i, ]); numbpart(p)^e) } \\ Andrew Howroyd, Nov 10 2025
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Gus Wiseman, Sep 05 2025
STATUS
approved
