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-10 of 14 results. Next

A122865 Expansion of chi(x) * phi(x^3) * psi(-x^3) in powers of x where chi(), phi(), psi() are Ramanujan theta functions.

Original entry on oeis.org

1, 1, 0, 2, 2, 1, 0, 0, 3, 0, 0, 2, 2, 2, 0, 0, 1, 2, 0, 2, 2, 1, 0, 0, 2, 0, 0, 2, 4, 0, 0, 0, 2, 3, 0, 2, 2, 0, 0, 0, 1, 0, 0, 4, 0, 2, 0, 0, 4, 2, 0, 0, 2, 2, 0, 0, 3, 0, 0, 2, 2, 0, 0, 0, 2, 1, 0, 2, 4, 2, 0, 0, 0, 0, 0, 2, 2, 2, 0, 0, 2, 2, 0, 4, 0, 1, 0
Offset: 0

Views

Author

Michael Somos, Sep 15 2006

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = 1 + x + 2*x^3 + 2*x^4 + x^5 + 3*x^8 + 2*x^11 + 2*x^12 + 2*x^13 + ...
G.f. = q + q^4 + 2*q^10 + 2*q^13 + q^16 + 3*q^25 + 2*q^34 + 2*q^37 + ...
		

Crossrefs

Programs

  • Mathematica
    phi[q_] := EllipticTheta[3, 0, q]; chi[q_] := ((1 - InverseEllipticNomeQ[q]) * InverseEllipticNomeQ[q]/(16*q))^(-1/24); psi[q_] := (1/2)*q^(-1/8)*EllipticTheta[ 2, 0, q^(1/2)]; s = Series[ chi[q]*phi[q^3]*psi[-q^3], {q, 0, 104}]; a[n_] := Coefficient[s, q, n];
    (* or *) a[n_] := If[n == 0, 1, Sum[Boole[Mod[d, 4] == 1] - Boole[Mod[d, 4] == 3], {d, Divisors[3n+1]}]]; Table[a[n], {n, 0, 104}] (* Jean-François Alcover, Feb 17 2015, after PARI code *)
    a[ n_] := If[ n < 0, 0, DivisorSum[ 3 n + 1, KroneckerSymbol[ -4, #] &]]; (* Michael Somos, Sep 02 2015 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x^2] QPochhammer[ -x^3] QPochhammer[ x^6], {x, 0, n}]; (* Michael Somos, Sep 02 2015 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^2 * eta(x^6 + A)^4 / (eta(x + A) * eta(x^3 + A) * eta(x^4 + A) * eta(x^12 + A)), n))};
    
  • PARI
    {a(n) = my(A, p, e); if(n <0, 0, n = 3*n+1; A = factor(n); prod(k=1, matsize(A)[1], [p, e] = A[k, ]; if(p==2, 1, p==3, -2*(-1)^e, p%4==1, e+1, 1-e%2)))};
    
  • PARI
    {a(n) = if( n<0, 0, n = 3*n+1; sumdiv(n, d, (d%4==1) - (d%4==3)))}; /* Michael Somos, Apr 19 2007 */

Formula

Expansion of chi(x) * f(x^3) * f(-x^6) in powers of x where chi(), f() are Ramanujan theta functions. - Michael Somos, Sep 02 2015
Expansion of q^(-1/3) * eta(q^2)^2 * et(q^6)^4 / (eta(q) * eta(q^3) * eta(q^4) * eta(q^12)) in powers of q.
Euler transform of period 12 sequence [1, -1, 2, 0, 1, -4, 1, 0, 2, -1, 1, -2, ...]. - Michael Somos, Apr 19 2007
G.f. is a period 1 Fourier series which satisfies f(-1 / (36 t)) = 2 (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A258228. - Michael Somos, Sep 02 2015
a(n) = A002654(3*n + 1) = A035154(3*n + 1) = A113446(3*n + 1) = A122864(3*n + 1) = A163746(3*n + 1).
a(n) = (-1)^n * A258277(n). a(2*n + 1) = A122856(n). - Michael Somos, Sep 02 2015
a(4*n) = A002175(n). a(4*n + 2) = 0. - Michael Somos, Jan 19 2017

A002175 Excess of number of divisors of 12n+1 of form 4k+1 over those of form 4k+3.

Original entry on oeis.org

1, 2, 3, 2, 1, 2, 2, 4, 2, 2, 1, 0, 4, 2, 3, 2, 2, 4, 0, 2, 2, 0, 4, 2, 3, 0, 2, 6, 2, 2, 1, 2, 0, 2, 2, 2, 2, 4, 2, 0, 4, 4, 4, 0, 1, 2, 0, 4, 2, 0, 2, 2, 5, 2, 0, 2, 2, 4, 4, 2, 0, 2, 4, 2, 2, 0, 4, 0, 0, 2, 3, 2, 4, 2, 0, 4, 0, 6, 2, 4, 1, 0, 4, 2, 2, 2, 2, 0, 0, 2, 0, 2, 8, 2, 2, 0, 2, 4, 0, 4, 2, 2, 3, 2, 2
Offset: 0

Views

Author

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Number of ways to write n as an ordered sum of 2 generalized pentagonal numbers. - Ilya Gutkovskiy, Aug 14 2017

Examples

			G.f. = 1 + 2*x + 3*x^2 + 2*x^3 + x^4 + 2*x^5 + 2*x^6 + 4*x^7 + 2*x^8 + 2*x^9 + ...
G.f. = q + 2*q^13 + 3*q^25 + 2*q^37 + q^49 + 2*q^61 + 2*q^73 + 4*q^85 + 2*q^97 + ...
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Maple
    series(mul( ( (1 + q^n)*(1 - q^(3*n))/(1 + q^(3*n)) )^2, n = 1..100), q, 101):
    seq(coeftayl(%, q = 0, n), n = 0..100); # Peter Bala, Jan 05 2025
  • Mathematica
    ed[n_]:=Module[{divs=Divisors[12n+1]},Count[divs,?(Mod[#,4] == 1&)]- Count[divs,?(Mod[#,4]==3&)]]; Array[ed,110,0] (* Harvey P. Dale, Jul 01 2012 *)
    a[ n_] := If[ n < 0, 0, With[ {m = 12 n + 1}, Sum[ KroneckerSymbol[ 4, d], {d, Divisors[m]}]]]; (* Michael Somos, Apr 23 2014 *)
    a[ n_] := SeriesCoefficient[ (QPochhammer[ x^2] QPochhammer[ x^3]^2 / (QPochhammer[ x] QPochhammer[ x^6]))^2, {x, 0, n}]; (* Michael Somos, Apr 23 2014 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 4, 0, x^3] / QPochhammer[ x, x^2])^2, {x, 0, n}]; (* Michael Somos, May 25 2015 *)
  • PARI
    {a(n) = if( n<0, 0, n = 12*n + 1; sumdiv( n, d, (d%4==1) - (d%4==3)))}; /* Michael Somos, Sep 19 2005 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^2 + A) * eta(x^3 + A)^2 / (eta(x + A) * eta(x^6 + A)))^2, n))}; /* Michael Somos, Jun 02 2012 */

Formula

Expansion of (phi(-x^3) / chi(-x))^2 in powers of x where phi(), chi() are Ramanujan theta functions.
Expansion of q^(-1/12) * (eta(q^2) * eta(q^3)^2 / (eta(q) * eta(q^6)))^2 in powers of q. - Michael Somos, Sep 19 2005
Euler transform of period 6 sequence [ 2, 0, -2, 0, 2, -2, ...]. - Michael Somos, Sep 19 2005
G.f. is a period 1 Fourier series which satisfies f(-1 / (72 t)) = 2 (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A258279. - Michael Somos, May 25 2015
From Michael Somos, Jun 02 2012: (Start)
a(n) = A008441(3*n) = A121363(3*n) = A122865(4*n) = A122856(8*n).
a(n) = A116604(6*n) = A125079(6*n) = A129447(6*n) = A138741(6*n).
From Michael Somos, May 25 2015: (Start)
a(n) = A258277(4*n) = A258278(8*n) = A258291(3*n).
a(n) = - A258210(12*n + 1) = A258228(12*n + 1) = A258256(12*n + 1).
2*a(n) = A258279(12*n + 1) = - A258292(12*n + 1). (End)
G.f.: (Sum_{k = -oo..oo} x^(k*(3*k-1)/2))^2. - Ilya Gutkovskiy, Aug 14 2017
G.f.: ( Product_{n >= 1} (1 + q^n)*(1 - q^(3*n))/(1 + q^(3*n)) )^2. - Peter Bala, Jan 05 2025

A258210 Expansion of f(-q) * f(-q^2) * chi(-q^3) in powers of q where chi(), f() are Ramanujan theta functions.

Original entry on oeis.org

1, -1, -2, 0, 1, 4, 0, 0, -2, -4, 2, 0, 0, -2, 0, 0, 1, 4, 4, 0, -4, 0, 0, 0, 0, -3, -4, 0, 0, 4, 0, 0, -2, 0, 2, 0, 4, -2, 0, 0, 2, 4, 0, 0, 0, -8, 0, 0, 0, -1, -6, 0, 2, 4, 0, 0, 0, 0, 2, 0, 0, -2, 0, 0, 1, 8, 0, 0, -4, 0, 0, 0, 4, -2, -4, 0, 0, 0, 0, 0, -4
Offset: 0

Views

Author

Michael Somos, May 23 2015

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Denoted by a_6(n) in Kassel and Reutenauer 2015. - Michael Somos, Jun 04 2015

Examples

			G.f. = 1 - q - 2*q^2 + q^4 + 4*q^5 - 2*q^8 - 4*q^9 + 2*q^10 - 2*q^13 + ...
		

Crossrefs

For the square of this series see A252650.

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ QPochhammer[ q]^2 / (QPochhammer[ q, q^6] QPochhammer[ q^5, q^6]), {q, 0, n}];
    a[ n_] := SeriesCoefficient[ (1/2) EllipticThetaPrime[ 1, 0, q^(1/2)] / EllipticTheta[ 1, Pi/6, q^(1/2)], {q, 0, n}];
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A) * eta(x^2 + A) * eta(x^3 + A) / eta(x^6 + A), n))};
    
  • PARI
    {a(n) = if( n<1, n==0, (-1)^n * (1 - (n%3==2)*3) * sumdiv(n, d, [0, 1, 2, -1][d%4 + 1] * if(d%9, 1, 4) * (-1)^((d%8==6) + n+d)))}; /* Michael Somos, Jun 04 2015 */

Formula

Expansion of f(-q)^2 * f(-q^6) / f(-q, -q^5) in powers of q where f(,) is Ramanujan's general theta function.
Expansion of eta(q) * eta(q^2) * eta(q^3) / eta(q^6) in powers of q.
Euler transform of period 6 sequence [ -1, -2, -2, -2, -1, -2, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (72 t)) = 12 (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A121444.
G.f.: Product_{k>0} (1 - x^k) * (1 - x^2*k) / (1 + x^(3*k)).
a(n) = (-1)^n * A258228(n). Convolution inverse of A077285.
a(4*n + 3) = 0. a(6*n + 2) = -2 * A122865(n). a(6*n + 4) = A122856(n). a(12*n + 1) = -1 * A002175(n).
a(9*n + 3) = a(9*n + 6) = 0. a(9*n) = A104794(n). a(3*n + 1) = -A258277(n). a(3*n + 2) = -2*A258278(n). - Michael Somos, May 01 2016
G.f.: Product_{i>0} 1/(1 + Sum_{j>0} j*x^(j*i)). - Seiichi Manyama, Oct 08 2017

A113652 Expansion of (1 - theta_4(q)^2) / 4 in powers of q.

Original entry on oeis.org

1, -1, 0, -1, 2, 0, 0, -1, 1, -2, 0, 0, 2, 0, 0, -1, 2, -1, 0, -2, 0, 0, 0, 0, 3, -2, 0, 0, 2, 0, 0, -1, 0, -2, 0, -1, 2, 0, 0, -2, 2, 0, 0, 0, 2, 0, 0, 0, 1, -3, 0, -2, 2, 0, 0, 0, 0, -2, 0, 0, 2, 0, 0, -1, 4, 0, 0, -2, 0, 0, 0, -1, 2, -2, 0, 0, 0, 0, 0, -2, 1, -2, 0, 0, 4, 0, 0, 0, 2, -2, 0, 0, 0, 0, 0, 0, 2, -1, 0, -3, 2, 0, 0, -2, 0
Offset: 1

Views

Author

Michael Somos, Nov 03 2005

Keywords

Examples

			G.f. = x - x^2 - x^4 + 2*x^5 - x^8 + x^9 - 2*x^10 + 2*x^13 - x^16 + 2*x^17 + ...
		

References

  • B. C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag, see p. 114 Entry 8(v).
  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 576.
  • J. M. Borwein and P. B. Borwein, Pi and the AGM, Wiley, 1987.
  • P. A. MacMahon, Combinatory Analysis, Cambridge Univ. Press, London and New York, Vol. 1, 1915 and Vol. 2, 1916; see vol. 2, p 28, Article 269.

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 1, 0, -(-1)^n DivisorSum[ n, KroneckerSymbol[ -4, #] &]]; (* Michael Somos, Jun 06 2015 *)
    a[ n_] := SeriesCoefficient[ (1 - EllipticTheta[4, 0, q]^2) / 4, {q, 0, n}]; (* Michael Somos, Jun 06 2015 *)
    a[ n_] := With[ {m = InverseEllipticNomeQ @ -q}, SeriesCoefficient[(1 - EllipticK[m] / (Pi/2)) / 4, {q, 0, n}]]; (* Michael Somos, Jun 06 2015 *)
  • PARI
    {a(n) = if( n<1, 0, -(-1)^n * sumdiv(n, d, kronecker( -4, d)))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, 0, A = factor(n); prod(k=1, matsize(A)[1], [p, e] = A[k,]; if( p==2, -1, p%4==1, e+1, !(e%2))))};
    
  • PARI
    {a(n) = if( n<1, 0, direuler(p=2, n, if( p==2, 1 - X/(1 - X), 1 / ((1 - X) * (1 - kronecker( -4, p)*X))) )[n])};
    
  • PARI
    {a(n) = my(A); if( n<1, 0, A = x * O(x^n); polcoeff( (1 - eta(x + A)^4 / eta(x^2 + A)^2) / 4, n))};

Formula

a(n) is multiplicative with a(2^e) = -1 if e>0, a(p^e) = e+1 if p == 1 (mod 4), (1 + (-1)^e)/2 if p == 3 (mod 4).
Expansion of (1 - eta(q)^4 / eta(q^2)^2) / 4 in powers of q.
Moebius transform is period 8 sequence [ 1, -2, -1, 0, 1, 2, -1, 0, ...].
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = u2 - 2*u3 + u6 - u1^2 + 3*u3^2 + 2*u1*u3 - 4*u2*u6.
G.f.: Sum_{k>0} -(-1)^k * x^((k^2 + k)/2) / (1 + x^k).
G.f.: Sum_{k>0} -(-1)^k * x^k / (1 + x^(2*k)).
G.f.: Sum_{k>0} -(-1)^k * x^(2*k - 1) / (1 + x^(2*k - 1)).
a(n) = -(-1)^n * A002654(n). a(n) = - A104794(n) / 4 unless n = 0.
a(2*n) = - A002654(n). a(3*n + 1) = A258277(n). a(3*n + 2) = - A258278(n). a(4*n + 1) = A008441(n). a(4*n + 3) = 0. a(6*n + 2) = - A122856(n). a(6*n + 4 ) = - A122856(n). - Michael Somos, Jun 06 2015
a(8*n + 1) = A113407(n). a(8*n + 5) = 2 * A053692(n). a(9*n + 3) = a(9*n + 6) = 0. - Michael Somos, Jun 06 2015

A258279 Expansion of psi(q)^2 * chi(-q^3)^2 in powers of q where psi(), chi() are Ramanujan theta functions.

Original entry on oeis.org

1, 2, 1, 0, -2, -2, 0, 0, 1, -4, -4, 0, 0, 4, 0, 0, -2, -2, 4, 0, 2, 0, 0, 0, 0, 6, 2, 0, 0, -2, 0, 0, 1, 0, -4, 0, 4, 4, 0, 0, -4, -2, 0, 0, 0, -8, 0, 0, 0, 2, 3, 0, -4, -2, 0, 0, 0, 0, -4, 0, 0, 4, 0, 0, -2, -4, 0, 0, 2, 0, 0, 0, 4, 4, 2, 0, 0, 0, 0, 0, 2
Offset: 0

Views

Author

Michael Somos, May 25 2015

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = 1 + 2*q + q^2 - 2*q^4 - 2*q^5 + q^8 - 4*q^9 - 4*q^10 + 4*q^13 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, Pi/6, q]^2, {q, 0, n}];
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^2 + A)^2 * eta(x^3 + A) / (eta(x + A) * eta(x^6 + A)))^2, n))};

Formula

Expansion of eta(q^2)^4 * eta(q^3)^2 / (eta(q)^2 * eta(q^6)^2) in powers of q.
Euler transform of period 6 sequence [ 2, -2, 0, -2, 2, -2, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (72 t)) = 36 (t/i) g(t) where q = exp(2 Pi i t) and g(t) is the g.f. for A002175.
G.f.: Product_{k>0} (1 - x^(2*k))^2 / (1 - x^k + x^(2*k))^2.
Convolution square of A089810.
a(2*n) = A258228(n). a(3*n + 1) = 2 * A258277(n). a(3*n + 2) = A258278(n). a(4*n + 3) = 0. a(6*n + 2) = A122865(n). a(6*n + 4) = -2 * A122856(n). a(12*n + 1) = 2 * A002175(n). a(12*n + 5) = -2 * A121444(n).
a(18*n) = A004018(n). a(18*n + 3) = a(18*n + 6) = a(18*n + 12) = 0.

A132003 Expansion of (phi(q^3) / phi(q)) * phi(-q^2) * phi(-q^6) in powers of q where phi() is a Ramanujan theta function.

Original entry on oeis.org

1, -2, 2, -2, 2, -4, 2, 0, 2, -2, 4, 0, 2, -4, 0, -4, 2, -4, 2, 0, 4, 0, 0, 0, 2, -6, 4, -2, 0, -4, 4, 0, 2, 0, 4, 0, 2, -4, 0, -4, 4, -4, 0, 0, 0, -4, 0, 0, 2, -2, 6, -4, 4, -4, 2, 0, 0, 0, 4, 0, 4, -4, 0, 0, 2, -8, 0, 0, 4, 0, 0, 0, 2, -4, 4, -6, 0, 0, 4, 0
Offset: 0

Views

Author

Michael Somos, Aug 06 2007

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = 1 - 2*x + 2*x^2 - 2*x^3 + 2*x^4 - 4*x^5 + 2*x^6 + 2*x^8 - 2*x^9 + 4*x^10 + ...
		

References

  • N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 85, Equation (32.72).

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 1, Boole[n == 0], -2 DivisorSum[ n, (-1)^(n + #) KroneckerSymbol[ -36, #] &]]; (* Michael Somos, Nov 01 2015 *)
    a[ n_] := If[ n < 1, Boole[n == 0], -2 Times @@ (Which[ # < 5, -(-1)^#, Mod[#, 4] == 3, 1 - Mod[#2, 2], True, #2 + 1] & @@@ FactorInteger @ n)]; (* Michael Somos, Nov 01 2015 *)
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q^3] EllipticTheta[ 4, 0, q^2] EllipticTheta[ 4, 0, q^6] / EllipticTheta[ 3, 0, q], {q, 0, n}]; (* Michael Somos, Nov 01 2015 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 4, 0, q]^2 + EllipticTheta[ 4, 0, q^3]^2) / 2, {q, 0, n}]; (* Michael Somos, Mar 05 2023 *)
  • PARI
    {a(n) = if( n<1, n==0, -2 * sumdiv(n, d, (-1)^(n+d) * kronecker(-36, d)))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 * eta(x^4 + A) * eta(x^6 + A)^7 / (eta(x^2 + A)^3 * eta(x^3 + A)^2 * eta(x^12 + A)^3), n))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, n==0, A = factor(n); -2 * prod(k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==3, 1, p==2, -1, p%4==1, e+1, 1-e%2)))};

Formula

Expansion of eta(q)^2 * eta(q^4) * eta(q^6)^7 / (eta(q^2)^3 * eta(q^3)^2 * eta(q^12)^3) in powers of q.
a(n) = -2*b(n) where b() is multiplicative with b(2^e) = 2*0^e - 1, b(3^e) = 1, b(p^e) = e+1 if p == 1 (mod 4), b(p^e) = (1 + (-1)^e)/2 if p == 3 (mod 4).
Euler transform of period 12 sequence [-2, 1, 0, 0, -2, -4, -2, 0, 0, 1, -2, -2, ...].
G.f.: 1 - 2 * Sum_{k>0} Kronecker(-36, k) * x^k / (1 + x^k).
a(n) = - A132004(n) unless n=0.
a(2*n) = A122857(n). a(2*n + 1) = -2 * A125079(n). a(3*n) = a(n). a(3*n + 1) = -2 * A258277(n). a(3*n + 2) = 2 * A258278(n). - Michael Somos, Nov 01 2015
a(12*n + 7) = a(12*n + 11) = 0. a(4*n + 1) = -2 * A008441(n).
a(n) = (-1)^n * A122857(n). Expansion of (phi(-q)^2 + phi(-q^3)^2) / 2 in powers of q where phi() is a Ramanujan theta function. - Michael Somos, Mar 05 2023

A132004 Expansion of (1 - phi(q^3) / phi(q) * phi(-q^2) * phi(-q^6)) / 2 in powers of q where phi() is a Ramanujan theta function.

Original entry on oeis.org

1, -1, 1, -1, 2, -1, 0, -1, 1, -2, 0, -1, 2, 0, 2, -1, 2, -1, 0, -2, 0, 0, 0, -1, 3, -2, 1, 0, 2, -2, 0, -1, 0, -2, 0, -1, 2, 0, 2, -2, 2, 0, 0, 0, 2, 0, 0, -1, 1, -3, 2, -2, 2, -1, 0, 0, 0, -2, 0, -2, 2, 0, 0, -1, 4, 0, 0, -2, 0, 0, 0, -1, 2, -2, 3, 0, 0, -2
Offset: 1

Views

Author

Michael Somos, Aug 06 2007

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = x - x^2 + x^3 - x^4 + 2*x^5 - x^6 - x^8 + x^9 - 2*x^10 - x^12 + 2*x^13 + ...
		

References

  • N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 85, Equation (32.72).

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 1, 0, DivisorSum[ n, (-1)^(n + #) KroneckerSymbol[ -36, #] &]]; (* Michael Somos, Nov 01 2015 *)
    a[ n_] := If[ n < 1, 0, Times @@ (Which[ # < 5, -(-1)^#, Mod[#, 4] == 3, 1 - Mod[#2, 2], True, #2 + 1] & @@@ FactorInteger @ n)]; (* Michael Somos, Nov 01 2015 *)
  • PARI
    {a(n) = if( n<1, 0, sumdiv( n, d, (-1)^(n+d) * kronecker( -36, d)))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (1 - eta(x + A)^2 * eta(x^4 + A) * eta(x^6 + A)^7 / (eta(x^2 + A)^3 * eta(x^3 + A)^2 * eta(x^12 + A)^3)) / 2, n))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, 0, A = factor(n); prod(k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==3, 1, p==2, -1, p%4==1, e+1, 1-e%2)))};

Formula

Expansion of (1 - eta(q)^2 * eta(q^4) * eta(q^6)^7 / (eta(q^2)^3 * eta(q^3)^2 * eta(q^12)^3)) / 2 in powers of q.
a(n) is multiplicative with a(2^e) = 2*0^e - 1, a(3^e) = 1, a(p^e) = e + 1 if p == 1 (mod 4), a(p^e) = (1 + (-1)^e) / 2 if p == 3 (mod 4).
G.f.: Sum_{k>0} x^k / (1 + x^k) * Kronecker(-36, k).
a(3*n) = a(n). -2 * a(n) = A132003(n) unless n = 0. a(2*n) = - A035154(n). a(2*n + 1) = A125079(n).
a(n) = (-1)^n * A035154(n). a(12*n + 7) = a(12*n + 11) = 0. - Michael Somos, Nov 01 2015
a(3*n + 1) = A258277(n). a(3*n + 2) = - A258278(n). a(4*n + 1) = A008441(n). a(4*n + 2) = - A125079(n). - Michael Somos, Nov 01 2015
a(6*n) = - A035154(n). a(6*n + 2) = - A122865(n). a(6*n + 4) = - A122856(n). - Michael Somos, Nov 01 2015
a(8*n + 1) = A113407(n). a(8*n + 5) = 2 * A053692(n). - Michael Somos, Nov 01 2015

A138952 Expansion of (eta(q^2)^7 * eta(q^3)^2 * eta(q^12) / (eta(q)^2 * eta(q^4)^3 * eta(q^6)^3) - 1) / 2 in powers of q.

Original entry on oeis.org

1, -1, -3, -1, 2, 3, 0, -1, 1, -2, 0, 3, 2, 0, -6, -1, 2, -1, 0, -2, 0, 0, 0, 3, 3, -2, -3, 0, 2, 6, 0, -1, 0, -2, 0, -1, 2, 0, -6, -2, 2, 0, 0, 0, 2, 0, 0, 3, 1, -3, -6, -2, 2, 3, 0, 0, 0, -2, 0, 6, 2, 0, 0, -1, 4, 0, 0, -2, 0, 0, 0, -1, 2, -2, -9, 0, 0, 6, 0
Offset: 1

Views

Author

Michael Somos, Apr 03 2008

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = q - q^2 - 3*q^3 - q^4 + 2*q^5 + 3*q^6 - q^8 + q^9 - 2*q^10 + 3*q^12 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 1, 0, -(-1)^n DivisorSum[ n, KroneckerSymbol[ -4, n/#] {1, 1, -2}[[Mod[#, 3, 1]]] &]]; (* Michael Somos, Sep 07 2015 *)
    a[ n_] := If[ n < 1, 0, Times @@ (Which[ # == 1, 1, # == 2, -1, # == 3, -1 + 2 (-1)^#2, Mod[#, 12] < 6, #2 + 1, True, 1 - Mod[#2, 2]] & @@@ FactorInteger@n)]; (* Michael Somos, Sep 07 2015 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 3, 0, q] EllipticTheta[ 4, 0, q^2] QPochhammer[ q^3] / QPochhammer[ -q^3] - 1) / 2, {q, 0, n}]; (* Michael Somos, Sep 07 2015 *)
  • PARI
    {a(n) = if( n<1, 0, -(-1)^n * sumdiv(n, d, kronecker(-4, n/d) * [-2, 1, 1][d%3 + 1]))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, 0, A = factor(n); prod( k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==2, -1, p==3, -1 + 2 * (-1)^e, p%12 < 6, e+1, 1-e%2 )))};

Formula

Expansion of (phi(q) * phi(-q^2) * chi(-q^3) / chi(q^3) - 1) / 2 in powers of q where phi(), chi() are Ramanujan theta functions.
Moebius transform is period 24 sequence [1, -2, -4, 0, 1, 8, -1, 0, 4, -2, -1, 0, 1, 2, -4, 0, 1, -8, -1, 0, 4, 2, -1, 0, ...].
a(n) is multiplicative with a(2^e) = -1 if e>0, a(3^e) = -1 + 2 * (-1)^e, a(p^e) = e+1 if p == 1, 5 (mod 12), a(p^e) = (1 + (-1)^e) / 2 if p == 7, 11 (mod 12).
a(12*n + 7) = a(12*n + 11) = 0.
a(n) = -(-1)^n * A138950(n). 2 * a(n) = A138951(n).
a(2*n) = - A138950(n). a(2*n + 1) = A116604(n). - Michael Somos, Sep 07 2015
a(3*n + 1) = A258277(n). a(3*n + 2) = - A258278(n). - Michael Somos, Sep 07 2015

A258291 Expansion of q^(-1/4) * eta(q) * eta(q^2) * eta(q^6) / eta(q^3) in powers of q.

Original entry on oeis.org

1, -1, -2, 2, -1, 0, 3, -1, 0, 2, -1, -4, 1, -1, 0, 2, -2, 0, 2, 0, -2, 4, -1, 0, 2, -1, 0, 2, -1, -4, 1, -2, 0, 0, -1, 0, 4, -1, -4, 2, 0, 0, 3, -1, 0, 2, -2, 0, 2, -1, 0, 4, 0, 0, 0, -2, -6, 2, -1, 0, 2, -1, 0, 0, -1, -4, 4, -1, 0, 2, -1, 0, 3, -1, 0, 0, -2
Offset: 0

Views

Author

Michael Somos, May 25 2015

Keywords

Examples

			G.f. = 1 - x - 2*x^2 + 2*x^3 - x^4 + 3*x^6 - x^7 + 2*x^9 - x^10 - 4*x^11 + ...
G.f. = q - q^5 - 2*q^9 + 2*q^13 - q^17 + 3*q^25 - q^29 + 2*q^37 - q^41 + ...
		

Crossrefs

Programs

  • Mathematica
    QP := QPochhammer; CoefficientList[Series[QP[q]*QP[q^2]*QP[q^6]/QP[q^3], {q, 0, 50}], q] (* G. C. Greubel, Aug 04 2018 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A) * eta(x^2 + A) * eta(x^6 + A) / eta(x^3 + A), n))};

Formula

Euler transform of period 6 sequence [ -1, -2, 0, -2, -1, -2, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (72 t)) = 9 (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A258277.
G.f.: Product_{k>0} (1 - x^k) * (1 - x^(2*k)) * (1 + x^(3*k)).
a(3*n) = A002175(n). a(3*n + 1) = - A121444(n). a(9*n + 2) = -2 * A008441(n). a(9*n + 5) = a(9*n + 8) = 0.

A258322 Expansion of phi(-q) * phi(-q^9) in powers of q where phi() is a Ramanujan theta function.

Original entry on oeis.org

1, -2, 0, 0, 2, 0, 0, 0, 0, -4, 4, 0, 0, -4, 0, 0, 2, 0, 4, 0, 0, 0, 0, 0, 0, -6, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 4, -4, 0, 0, 4, 0, 0, 0, 0, -8, 0, 0, 0, -2, 0, 0, 4, 0, 0, 0, 0, 0, 4, 0, 0, -4, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 4, -4, 0, 0, 0, 0, 0, 0, 0, -4, 4
Offset: 0

Views

Author

Michael Somos, May 26 2015

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = 1 - 2*q + 2*q^4 - 4*q^9 + 4*q^10 - 4*q^13 + 2*q^16 + 4*q^18 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ QPochhammer[ q]^2 QPochhammer[ q^9]^2 / (QPochhammer[ q^2] QPochhammer[ q^18]), {q, 0, n}];
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 * eta(x^9 + A)^2 / (eta(x^2 + A) * eta(x^18 + A)), n))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, n==0, if( n%3 == 2, 0, A = factor(n); -2 * prod( k=1, matsize(A)[1], [p, e] = A[k,]; if( p==2, -1, p%4==3, if(p==3, 2, 1)*(1-e%2), e+1))))};

Formula

Expansion of eta(q)^2 * eta(q^9)^2 / (eta(q^2) * eta(q^18)) in powers of q.
Euler transform of period 18 sequence [ -2, -1, -2, -1, -2, -1, -2, -1, -4, -1, -2, -1, -2, -1, -2, -1, -2, -2, ...]. - Michael Somos, May 26 2015
a(n) = -2*b(n) where b() is multiplicative with a(0) = 1, b(2^e) = -1 if e>0, b(3^e) = 1 + (-1)^e if e>0, b(p^e) = (1 + (-1)^e) / 2 if p == 3 (mod 4), b(p^e) = e+1 if p == 1 (mod 4).
a(n) = -2 * A257900(n) unless n = 0 or n == 2 (mod 3).
a(3*n + 1) = -2 * A258277(n). a(3*n + 2) = 0. a(9*n) = -4 * A113652(n). a(9*n + 3) = a(9*n + 6) = 0.
Sum_{k=1..n} abs(a(k)) ~ (Pi/3) * n. - Amiram Eldar, Jan 29 2024
Showing 1-10 of 14 results. Next