|
%I #33 Feb 01 2024 08:51:55
%S 1,2,4,5,8,10,14,16,23,26,34,40,50,58,74,83,102,120,142,164,198,226,
%T 266,308,359,412,482,548,634,730,834,950,1094,1240,1416,1609,1826,
%U 2068,2350,2648,2994,3382,3806,4280,4826,5408,6070,6806,7619,8522,9534,10632
%N Sum of smallest parts of all partitions of n into distinct parts.
%H Alois P. Heinz, <a href="/A092265/b092265.txt">Table of n, a(n) for n = 1..10000</a>
%F G.f.: Sum_{n >= 1} (-1 + Product_k >= n} 1 + x^k).
%F G.f.: Sum_{n >= 1} n*x^n*Product_{k >= n+1} (1 + x^k). - _Joerg Arndt_, Jan 29 2011
%F G.f.: Sum_{k >= 1} x^(k*(k+1)/2)/(1 - x^k)/Product_{i = 1..k} (1 - x^i). - _Vladeta Jovovic_, Aug 10 2004
%F Conjecture: a(n) = A034296(n) + A237665(n+1). - _George Beck_, May 06 2017
%F a(n) ~ exp(Pi*sqrt(n/3)) / (2 * 3^(1/4) * n^(3/4)). - _Vaclav Kotesovec_, May 20 2018
%p b:= proc(n, i) option remember; `if`(n=0, 1,
%p `if`(i>n, 0, b(n,i+1)+b(n-i, i+1)))
%p end:
%p a:= n-> add(j*b(n-j, j+1), j=1..n):
%p seq(a(n), n=1..80); # _Alois P. Heinz_, Feb 03 2016
%t b[n_, i_] := b[n, i] = If[n == 0, 1, If[i > n, 0, b[n, i + 1] + b[n - i, i + 1]]]; a[n_] := Sum[j*b[n - j, j + 1], {j, 1, n}]; Table[a[n], {n, 1, 80}] (* _Jean-François Alcover_, Jan 21 2017, after _Alois P. Heinz_ *)
%Y Cf. A046746, A005895, A006128, A092269, A092319, A092316, A034296, A237665.
%Y Cf. A026832, A336902, A336903.
%K easy,nonn
%O 1,2
%A _Vladeta Jovovic_, Feb 14 2004
%E More terms from Pab Ter (pabrlos(AT)yahoo.com), May 25 2004
|
