|
%I #25 Aug 14 2018 00:25:33
%S 1,3,4,6,8,8,11,13,13,15,17,19,20,22,22,24,28,28,30,30,31,35,37,37,39,
%T 41,41,43,45,47,48,52,52,52,54,54,58,60,62,64,64,64,67,69,69,71,75,75,
%U 77,79,79,83,83,83,83,87,90,92,94,94,96,98,98,98,100,102,106,108,108,110,112,112,115,117,117,117,121,121,123
%N G.f.: 1/(1-x) * Product_{n>=1} (1 - x^(2*n))^2/(1 - x^(2*n-1))^2.
%C Equals the partial sums of A008441, where A008441(n) is the number of ways of writing n as the sum of 2 triangular numbers.
%H G. C. Greubel, <a href="/A199015/b199015.txt">Table of n, a(n) for n = 0..1000</a>
%F Lim_{n->infinity} a(n)/n = Pi/2.
%F a(n) = Sum_{k=0..n} Sum_{d|4*k+1} (-1)^floor(d/2). - _Michael Somos_ [see A008441]
%F G.f.: 1/(1-x) * Sum_{n>=0} x^n/(1 - x^(4*n + 1)). - _Michael Somos_ [see A008441]
%F G.f.: theta_2(sqrt(x))^2/(4*x^(1/4)*(1 - x)), where theta_2() is the Jacobi theta function. - _Ilya Gutkovskiy_, Apr 15 2018
%e G.f.: A(x) = 1 + 3*x + 4*x^2 + 6*x^3 + 8*x^4 + 8*x^5 + 11*x^6 + 13*x^7 + ...
%e where the g.f. equals the product:
%e A(x) = 1/(1-x) * (1-x^2)^2/(1-x)^2 * (1-x^4)^2/(1-x^3)^2 * (1-x^6)^2/(1-x^5)^2 * ...
%e Illustrate the limit a(n)/n = Pi/2:
%e a(10)/10 = 1.7, a(10^2)/10^2 = 1.58, a(10^3)/10^3 = 1.574, a(10^4)/10^4 = 1.5704, a(10^5)/10^5 = 1.57086, a(10^6)/10^6 = 1.570784, a(10^7)/10^7 = 1.5707972, ...
%t CoefficientList[Series[EllipticTheta[2, 0, Sqrt[x]]^2/(4*x^(1/4)*(1 - x)), {x, 0, 50}], x] (* _G. C. Greubel_, Aug 13 2018 *)
%o (PARI) {a(n)=sum(k=0,n,sumdiv(4*k+1, d, (-1)^(d\2)))}
%o (PARI) {a(n)=polcoeff(1/(1-x)*prod(m=1,n\2+1,(1-x^(2*m))/(1-x^(2*m-1)+x*O(x^n)))^2,n)}
%Y Cf. A008441.
%K nonn
%O 0,2
%A _Paul D. Hanna_, Nov 02 2011
|
