%I #8 Jan 03 2018 17:55:04
%S 1,-9,54,-45,-1404,6264,1890,-76050,187272,-8865,-1540944,6200280,
%T -1621620,-51195330,109055700,42125400,-868685040,2946297888,74093670,
%U -21584605122,44912353824,-17376284250,-302040439920,1069478852112,249392931480,-7095191496489
%N a(n) = Pell(n)*A109041(n) for n>=1, with a(0)=1, where A109041 lists the coefficients in eta(q)^9/eta(q^3)^3.
%C Compare the g.f. to the Lambert series of A109041:
%C 1 - 9*Sum_{n>=1} Kronecker(n,3)*n^2*x^n/(1-x^n).
%H G. C. Greubel, <a href="/A209453/b209453.txt">Table of n, a(n) for n = 0..1000</a>
%F G.f.: 1 - 9*Sum_{n>=1} Pell(n)*Kronecker(n,3)*n^2*x^n/(1 - A002203(n)*x^n + (-1)^n*x^(2*n)), where A002203(n) = Pell(n-1) + Pell(n+1).
%e G.f.: A(x) = 1 - 9*x + 54*x^2 - 45*x^3 - 1404*x^4 + 6264*x^5 + 1890*x^6 +...
%e where A(x) = 1 - 1*9*x + 2*27*x^2 - 5*9*x^3 - 12*117*x^4 + 29*216*x^5 + 70*27*x^6 - 169*450*x^7 + 408*459*x^8 +...+ Pell(n)*A109041(n)*^n +...
%e The g.f. is also given by the identity:
%e A(x) = 1 - 9*( 1*1*x/(1-2*x-x^2) - 2*4*x^2/(1-6*x^2+x^4) + 12*16*x^4/(1-34*x^4+x^8) - 29*25*x^5/(1-82*x^5-x^10) + 169*49*x^7/(1-478*x^7-x^14) - 408*64*x^8/(1-1154*x^8+x^16) +...).
%e The values of the symbol Kronecker(n,3) repeat [1,-1,0, ...].
%t A109041[n_]:= If[n < 1, Boole[n == 0], -9 DivisorSum[n, #^2 KroneckerSymbol[-3, #] &]]; Join[{1}, Table[Fibonacci[n, 2]*A109041[n], {n, 1, 50}]] (* _G. C. Greubel_, Jan 02 2018 *)
%o (PARI) {Pell(n)=polcoeff(x/(1-2*x-x^2+x*O(x^n)),n)}
%o {A002203(n)=Pell(n-1)+Pell(n+1)}
%o {a(n)=polcoeff(1 - 9*sum(m=1,n,Pell(m)*kronecker(m,3)*m^2*x^m/(1-A002203(m)*x^m+(-1)^m*x^(2*m) +x*O(x^n))),n)}
%o for(n=0,40,print1(a(n),", "))
%Y Cf. A109041, A205973, A209452, A209454, A204270, A000129 (Pell), A002203.
%K sign
%O 0,2
%A _Paul D. Hanna_, Mar 10 2012
|