|
|
A047966
|
|
a(n) = Sum_{ d divides n } q(d), where q(d) = A000009 = number of partitions of d into distinct parts.
|
|
126
|
|
|
1, 2, 3, 4, 4, 8, 6, 10, 11, 15, 13, 25, 19, 29, 33, 42, 39, 62, 55, 81, 84, 103, 105, 153, 146, 185, 203, 253, 257, 344, 341, 432, 463, 552, 594, 747, 761, 920, 1003, 1200, 1261, 1537, 1611, 1921, 2089, 2410, 2591, 3095, 3270, 3815, 4138, 4769, 5121, 5972, 6394, 7367, 7974, 9066, 9793, 11305, 12077, 13736, 14940
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Number of partitions of n such that every part occurs with the same multiplicity. - Vladeta Jovovic, Oct 22 2004
Christopher and Christober call such partitions uniform. - Gus Wiseman, Apr 16 2018
Equals inverse Mobius transform (A051731) * A000009, where the latter begins (1, 1, 2, 2, 3, 4, 5, 6, 8, ...). - Gary W. Adamson, Jun 08 2009
|
|
LINKS
|
|
|
FORMULA
|
G.f.: Sum_{k>0} (-1+Product_{i>0} (1+z^(k*i))). - Vladeta Jovovic, Jun 22 2003
|
|
EXAMPLE
|
The a(6) = 8 uniform partitions are (6), (51), (42), (33), (321), (222), (2211), (111111). - Gus Wiseman, Apr 16 2018
|
|
MAPLE
|
with(numtheory):
b:= proc(n) option remember; `if`(n=0, 1, add(add(
`if`(d::odd, d, 0), d=divisors(j))*b(n-j), j=1..n)/n)
end:
a:= n-> add(b(d), d=divisors(n)):
|
|
MATHEMATICA
|
b[n_] := b[n] = If[n==0, 1, Sum[DivisorSum[j, If[OddQ[#], #, 0]&]*b[n-j], {j, 1, n}]/n]; a[n_] := DivisorSum[n, b]; Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Dec 06 2016 after Alois P. Heinz *)
Table[DivisorSum[n, PartitionsQ], {n, 20}] (* Gus Wiseman, Apr 16 2018 *)
|
|
PROG
|
(PARI)
N = 66; q='q+O('q^N);
Vec( sum(e=1, N, D(q^e)-1) ) \\ Joerg Arndt, Mar 27 2014
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|