|
|
A092269
|
|
Spt function: total number of smallest parts (counted with multiplicity) in all partitions of n.
|
|
65
|
|
|
1, 3, 5, 10, 14, 26, 35, 57, 80, 119, 161, 238, 315, 440, 589, 801, 1048, 1407, 1820, 2399, 3087, 3998, 5092, 6545, 8263, 10486, 13165, 16562, 20630, 25773, 31897, 39546, 48692, 59960, 73423, 89937, 109553, 133439, 161840, 196168, 236843, 285816, 343667, 412950, 494702, 592063, 706671
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
G.f.: Sum_{n>=1} x^n/(1-x^n) * Product_{k>=n} 1/(1-x^k).
Asymptotics (Bringmann-Mahlburg, 2009): a(n) ~ exp(Pi*sqrt(2*n/3)) / (Pi*sqrt(8*n)) ~ sqrt(6*n)*A000041(n)/Pi. - Vaclav Kotesovec, Jul 30 2017
|
|
EXAMPLE
|
Partitions of 4 are [1,1,1,1], [1,1,2], [2,2], [1,3], [4]. 1 appears 4 times in the first, 1 twice in the second, 2 twice in the third, etc.; thus a(4)=4+2+2+1+1=10.
|
|
MAPLE
|
b:= proc(n, i) option remember; `if`(n=0 or i=1, n,
`if`(irem(n, i, 'r')=0, r, 0)+add(b(n-i*j, i-1), j=0..n/i))
end:
a:= n-> b(n, n):
|
|
MATHEMATICA
|
terms = 47; gf = Sum[x^n/(1 - x^n)*Product[1/(1 - x^k), {k, n, terms}], {n, 1, terms}]; CoefficientList[ Series[gf, {x, 0, terms}], x] // Rest (* Jean-François Alcover, Jan 17 2013 *)
b[n_, i_] := b[n, i] = If[n==0 || i==1, n, {q, r} = QuotientRemainder[n, i]; If[r==0, q, 0] + Sum[b[n-i*j, i-1], {j, 0, n/i}]]; a[n_] := b[n, n]; Table[a[n], {n, 1, 60}] (* Jean-François Alcover, Nov 23 2015, after Alois P. Heinz *)
|
|
PROG
|
(PARI)
N = 66; x = 'x + O('x^N);
gf = sum(n=1, N, x^n/(1-x^n) * prod(k=n, N, 1/(1-x^k) ) );
v = Vec(gf)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
More terms from Pab Ter (pabrlos(AT)yahoo.com), May 25 2004
|
|
STATUS
|
approved
|
|
|
|