cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A205972 a(n) = Fibonacci(n)*A122859(n) for n>=1, with a(0)=1, where A122859 lists the coefficients in phi(-q)^3/phi(-q^3) and phi() is a Ramanujan theta function.

Original entry on oeis.org

1, -6, 12, -12, -18, 0, 96, -156, 252, -204, 0, 0, -864, -2796, 9048, 0, -5922, 0, 31008, -50172, 0, -131352, 0, 0, 556416, -450150, 2913432, -1178508, -3813732, 0, 0, -16155228, 26139708, 0, 0, 0, -89582112, -289893804, 938116056, -758951832, 0, 0, 6429943104
Offset: 0

Views

Author

Paul D. Hanna, Feb 04 2012

Keywords

Comments

Compare the g.f. to the Lambert series of A122859:
1 - 6*Sum_{n>=1} Kronecker(n,3)*x^n/(1+x^n).

Examples

			G.f.: A(x) = 1 - 6*x + 12*x^2 - 12*x^3 - 18*x^4 + 96*x^6 - 156*x^7 +...
where A(x) = 1 - 1*6*x + 1*12*x^2 - 2*6*x^3 - 3*6*x^4 + 8*12*x^6 - 13*12*x^7 + 21*12*x^8 - 34*6*x^9 +...+ Fibonacci(n)*A122859(n)*x^n +...
The g.f. is also given by the identity:
A(x) = 1 - 6*( 1*x/(1+x-x^2) - 1*x^2/(1+3*x^2+x^4) + 3*x^4/(1+7*x^4+x^8) - 5*x^5/(1+11*x^5-x^10) + 13*x^7/(1+29*x^7-x^14) - 21*x^8/(1+47*x^8+x^16) +...).
The values of the symbol Kronecker(n,3) repeat [1,-1,0, ...].
		

Crossrefs

Cf. A209452 (Pell variant).

Programs

  • Mathematica
    A122859:= CoefficientList[Series[EllipticTheta[3, 0, -q]^3/EllipticTheta[3, 0, -q^3], {q, 0, 60}], q]; Table[If[n == 1, 1, Fibonacci[n - 1]*A122859[[n]]], {n, 1, 50}] (* G. C. Greubel, Dec 03 2017 *)
  • PARI
    {Lucas(n)=fibonacci(n-1)+fibonacci(n+1)}
    {a(n)=polcoeff(1 - 6*sum(m=1,n,fibonacci(m)*kronecker(m,3)*x^m/(1+Lucas(m)*x^m+(-1)^m*x^(2*m) +x*O(x^n))),n)}
    for(n=0,40,print1(a(n),", "))

Formula

G.f.: 1 - 6*Sum_{n>=1} Fibonacci(n)*Kronecker(n,3)*x^n/(1 + Lucas(n)*x^n + (-1)^n*x^(2*n)).

A209451 a(n) = Pell(n)*A034896(n) for n >= 1, with a(0)=1, where A034896 lists the number of solutions to a^2 + b^2 + 3*c^2 + 3*d^2 = n.

Original entry on oeis.org

1, 4, 8, 20, 240, 696, 280, 5408, 21216, 3940, 57072, 275568, 277200, 1873816, 2585024, 4680600, 54616512, 81841608, 10976840, 530008720, 1919331360, 1235646880, 4474673184, 21605633376, 28253665440, 162655527004, 177341693872, 30581480180, 2953208968320
Offset: 0

Views

Author

Paul D. Hanna, Mar 10 2012

Keywords

Comments

Compare g.f. to the Lambert series of A034896:
1 + 4*Sum_{n>=1} Chi(n,3)*n*x^n/(1 - (-x)^n).
Here Chi(n,3) = principal Dirichlet character modulo 3.

Examples

			G.f.: A(x) = 1 + 4*x + 8*x^2 + 20*x^3 + 240*x^4 + 696*x^5 + 280*x^6 + ...
where A(x) = 1 + 1*4*x + 2*4*x^2 + 5*4*x^3 + 12*20*x^4 + 29*24*x^5 + 70*4*x^6 + ... + Pell(n)*A034896(n)*x^n + ...
The g.f. is also given by the identity:
A(x) = 1 + 4*( 1*1*x/(1+2*x-x^2) + 2*2*x^2/(1-6*x^2+x^4) + 12*4*x^4/(1-34*x^4+x^8) + 29*5*x^5/(1+82*x^5-x^10) + 169*7*x^7/(1+478*x^7-x^14) + 408*8*x^8/(1-1154*x^8+x^16) + ...).
The values of the Dirichlet character Chi(n,3) repeat [1,1,0,...].
		

Crossrefs

Programs

  • Mathematica
    A034896[n_]:= SeriesCoefficient[(EllipticTheta[3, 0, q]*EllipticTheta[3, 0, q^3])^2, {q, 0, n}]; Join[{1}, Table[Fibonacci[n, 2]*A034896[n], {n, 1, 50}]] (* G. C. Greubel, Dec 24 2017 *)
  • PARI
    {Pell(n)=polcoeff(x/(1-2*x-x^2+x*O(x^n)),n)}
    {A002203(n)=Pell(n-1)+Pell(n+1)}
    {a(n)=polcoeff(1 + 4*sum(m=1,n,Pell(m)*kronecker(m,3)^2*m*x^m/(1-A002203(m)*(-x)^m+(-1)^m*x^(2*m) +x*O(x^n))),n)}
    for(n=0,61,print1(a(n),", "))

Formula

G.f.: 1 + 4*Sum_{n>=1} Pell(n)*Chi(n,3)*n*x^n/(1 - A002203(n)*(-x)^n + (-1)^n*x^(2*n)), where A002203(n) = Pell(n-1) + Pell(n+1).

A209453 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.

Original entry on oeis.org

1, -9, 54, -45, -1404, 6264, 1890, -76050, 187272, -8865, -1540944, 6200280, -1621620, -51195330, 109055700, 42125400, -868685040, 2946297888, 74093670, -21584605122, 44912353824, -17376284250, -302040439920, 1069478852112, 249392931480, -7095191496489
Offset: 0

Views

Author

Paul D. Hanna, Mar 10 2012

Keywords

Comments

Compare the g.f. to the Lambert series of A109041:
1 - 9*Sum_{n>=1} Kronecker(n,3)*n^2*x^n/(1-x^n).

Examples

			G.f.: A(x) = 1 - 9*x + 54*x^2 - 45*x^3 - 1404*x^4 + 6264*x^5 + 1890*x^6 +...
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 +...
The g.f. is also given by the identity:
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) +...).
The values of the symbol Kronecker(n,3) repeat [1,-1,0, ...].
		

Crossrefs

Programs

  • Mathematica
    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 *)
  • PARI
    {Pell(n)=polcoeff(x/(1-2*x-x^2+x*O(x^n)),n)}
    {A002203(n)=Pell(n-1)+Pell(n+1)}
    {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)}
    for(n=0,40,print1(a(n),", "))

Formula

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).
Showing 1-3 of 3 results.