%I #13 Oct 10 2024 05:35:17
%S 1,2,2,0,4,4,4,4,0,8,8,8,16,8,8,8,16,16,16,32,32,32,32,32,16,48,48,32,
%T 64,64,96,96,96,96,96,96,128,128,96,160,160,192,256,256,256,320,320,
%U 320,320,384,384,384,448,384,512,576,704,704,768,896,896,1024,1024,1024
%N G.f.: Sum_{k>=0} 2^k * x^(k^2) * Product_{j=1..k} (1 + x^j).
%H Vaclav Kotesovec, <a href="/A376945/b376945.txt">Table of n, a(n) for n = 0..10000</a>
%F a(n) ~ (1+r) * exp(sqrt((2*log(2)^2 + 8*log(2)*log(r) + 12*log(r)^2 + 4*polylog(2, 1/(1+r)) - Pi^2/3)*n)) / (2*sqrt((3*r + 2)*n)), where r = ((46 - 6*sqrt(57))^(1/3) + (46 + 6*sqrt(57))^(1/3) - 2)/6 is the real root of the equation 2*r^2*(1+r) = 1 (A273065).
%t nmax = 80; CoefficientList[Series[Sum[2^k * x^(k^2) * Product[1+x^j, {j, 1, k}], {k, 0, Sqrt[nmax]}], {x, 0, nmax}], x]
%t nmax = 80; p = 1; s = 1; Do[p = Normal[Series[2*p*(1 + x^k) * x^(2*k - 1), {x, 0, nmax}]]; s += p; , {k, 1, Sqrt[nmax]}]; Take[CoefficientList[s, x], nmax + 1]
%Y Cf. A273065, A306734, A376943, A376944, A376948.
%K nonn
%O 0,2
%A _Vaclav Kotesovec_, Oct 10 2024