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 13 results. Next

A008441 Number of ways of writing n as the sum of 2 triangular numbers.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700). The present sequence gives the expansion coefficients of psi(q)^2.
Also the number of positive odd solutions to equation x^2 + y^2 = 8*n + 2. - Seiichi Manyama, May 28 2017

Examples

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

References

  • B. C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag. See p. 139 Example (iv).
  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 102.
  • R. W. Gosper, Strip Mining in the Abandoned Orefields of Nineteenth Century Mathematics, in Computers in Mathematics (Ed. D. V. Chudnovsky and R. D. Jenks). New York: Dekker, 1990. See p. 279.
  • R. W. Gosper, Experiments and discoveries in q-trigonometry, in Symbolic Computation, Number Theory, Special Functions, Physics and Combinatorics. Editors: F. G. Garvan and M. E. H. Ismail. Kluwer, Dordrecht, Netherlands, 2001, pp. 79-105. [See Pi_q.]
  • P. A. MacMahon, Combinatory Analysis, Cambridge Univ. Press, London and New York, Vol. 1, 1915 and Vol. 2, 1916. See vol. 2, p 31, Article 272.
  • Ivan Niven, Herbert S. Zuckerman and Hugh L. Montgomery, An Introduction to the Theory Of Numbers, Fifth Edition, John Wiley and Sons, Inc., NY 1991, p. 165.

Crossrefs

Cf. A004020, A005883, A104794, A052343, A199015 (partial sums).
Number of ways of writing n as a sum of k triangular numbers, for k=1,...: A010054, A008441, A008443, A008438, A008439, A008440, A226252, A007331, A226253, A226254, A226255, A014787, A014809.
Cf. A274621 (reciprocal series).

Programs

  • Haskell
    a052343 = (flip div 2) . (+ 1) . a008441
    -- Reinhard Zumkeller, Jul 25 2014
    
  • Magma
    A := Basis( ModularForms( Gamma1(8), 1), 420); A[2]; /* Michael Somos, Jan 31 2015 */
  • Maple
    sigmamr := proc(n,m,r) local a,d ; a := 0 ; for d in numtheory[divisors](n) do if modp(d,m) = r then a := a+1 ; end if; end do: a; end proc:
    A002654 := proc(n) sigmamr(n,4,1)-sigmamr(n,4,3) ; end proc:
    A008441 := proc(n) A002654(4*n+1) ; end proc:
    seq(A008441(n),n=0..90) ; # R. J. Mathar, Mar 23 2011
  • Mathematica
    Plus@@((-1)^(1/2 (Divisors[4#+1]-1)))& /@ Range[0, 104] (* Ant King, Dec 02 2010 *)
    a[ n_] := SeriesCoefficient[ (1/2) EllipticTheta[ 2, 0, q] EllipticTheta[ 3, 0, q], {q, 0, n + 1/4}]; (* Michael Somos, Jun 19 2012 *)
    a[ n_] := SeriesCoefficient[ (1/4) EllipticTheta[ 2, 0, q]^2, {q, 0, 2 n + 1/2}]; (* Michael Somos, Jun 19 2012 *)
    a[ n_] := If[ n < 0, 0, DivisorSum[ 4 n + 1, (-1)^Quotient[#, 2] &]];  (* Michael Somos, Jun 08 2014 *)
    QP = QPochhammer; s = QP[q^2]^4/QP[q]^2 + O[q]^100; CoefficientList[s, q] (* Jean-François Alcover, Nov 27 2015, adapted from PARI *)
    TriangleQ[n_] := IntegerQ@Sqrt[8n +1]; Table[Count[FrobeniusSolve[{1, 1}, n], {?TriangleQ}], {n, 0, 104}] (* Robert G. Wilson v, Apr 15 2017 *)
  • PARI
    {a(n) = if( n<1, n==0, polcoeff( sum(k=0, (sqrtint(8*n + 1) - 1)\2, x^(k * (k+1)/2), x * O(x^n))^2, n) )};
    
  • PARI
    {a(n) = if( n<0, 0, n = 4*n + 1; sumdiv(n, d, (-1)^(d\2)))}; /* Michael Somos, Sep 02 2005 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^4 / eta(x + A)^2, n))};
    
  • PARI
    {a(n) = if( n<0, 0, n = 4*n + 1; sumdiv( n, d, (d%4==1) - (d%4==3)))}; /* Michael Somos, Sep 14 2005 */
    
  • PARI
    { my(q='q+O('q^166)); Vec(eta(q^2)^4 / eta(q)^2) } \\ Joerg Arndt, Apr 16 2017
    
  • Sage
    ModularForms( Gamma1(8), 1, prec=420).1; # Michael Somos, Jun 08 2014
    

Formula

This sequence is the quadrisection of many sequences. Here are two examples:
a(n) = A002654(4n+1), the difference between the number of divisors of 4*n+1 of form 4*k+1 and the number of form 4*k-1. - David Broadhurst, Oct 20 2002
a(n) = b(4*n + 1), where b(n) is multiplicative and b(2^e) = 0^e, b(p^e) = (1 + (-1)^e) / 2 if p == 3 (mod 4), b(p^e) = e+1 if p == 1 (mod 4). - Michael Somos, Sep 14 2005
G.f.: (Sum_{k>=0} x^((k^2 + k)/2))^2 = (Sum_{k>=0} x^(k^2 + k)) * (Sum_{k in Z} x^(k^2)).
Expansion of Jacobi theta (theta_2(0, sqrt(q)))^2 / (4 * q^(1/4)).
Sum[d|(4n+1), (-1)^((d-1)/2) ].
Given g.f. A(x), then B(q) = q * A(q^4) satisfies 0 = f(B(q), B(q^2), B(q^4)) where f(u, v, w) = v^3 + 4 * v * w^2 - u^2 * w. - Michael Somos, Sep 14 2005
Given g.f. A(x), then B(q) = q * A(q^4) satisfies 0 = f(B(q), B(q^2), B(q^3), B(q^6)) where f(u1, u2, u3, u6) = u1 * u3 - (u2 - u6) * (u2 + 3*u6). - Michael Somos, Sep 14 2005
Expansion of Jacobi k/(4*q^(1/2)) * (2/Pi)* K(k) in powers of q^2. - Michael Somos, Sep 14 2005. Convolution of A001938 and A004018. This appears in the denominator of the Jacobi sn and cn formula given in the Abramowitz-Stegun reference, p. 575, 16.23.1 and 16.23.2, where m=k^2. - Wolfdieter Lang, Jul 05 2016
G.f.: Sum_{k>=0} a(k) * x^(2*k) = Sum_{k>=0} x^k / (1 + x^(2*k + 1)).
G.f.: Sum_{k in Z} x^k / (1 - x^(4*k + 1)). - Michael Somos, Nov 03 2005
Expansion of psi(x)^2 = phi(x) * psi(x^2) in powers of x where phi(), psi() are Ramanujan theta functions.
Moebius transform is period 8 sequence [ 1, -1, -1, 0, 1, 1, -1, 0, ...]. - Michael Somos, Jan 25 2008
G.f. is a period 1 Fourier series which satisfies f(-1 / (8 t)) = 1/2 (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A104794.
Euler transform of period 2 sequence [ 2, -2, ...].
G.f.: q^(-1/4) * eta(q^2)^4 / eta(q)^2. See also the Fine reference.
a(n) = Sum_{k=0..n} A010054(k)*A010054(n-k). - Reinhard Zumkeller, Nov 03 2009
A004020(n) = 2 * a(n). A005883(n) = 4 * a(n).
Convolution square of A010054.
G.f.: Product_{k>0} (1 - x^(2*k))^2 / (1 - x^(2*k-1))^2.
a(2*n) = A113407(n). a(2*n + 1) = A053692(n). a(3*n) = A002175(n). a(3*n + 1) = 2 * A121444(n). a(9*n + 2) = a(n). a(9*n + 5) = a(9*n + 8) = 0. - Michael Somos, Jun 08 2014
G.f.: exp( Sum_{n>=1} 2*(x^n/n) / (1 + x^n) ). - Paul D. Hanna, Mar 01 2016
a(n) = A001826(2+8*n) - A001842(2+8*n), the difference between the number of divisors 1 (mod 4) and 3 (mod 4) of 2+8*n. See the Ono et al. link, Corollary 1, or directly the Niven et al. reference, p. 165, Corollary (3.23). - Wolfdieter Lang, Jan 11 2017
Expansion of continued fraction 1 / (1 - x^1 + x^1*(1 - x^1)^2 / (1 - x^3 + x^2*(1 - x^2)^2 / (1 - x^5 + x^3*(1 - x^3)^2 / ...))) in powers of x^2. - Michael Somos, Apr 20 2017
Given g.f. A(x), and B(x) is the g.f. for A079006, then B(x) = A(x^2) / A(x) and B(x) * B(x^2) * B(x^4) * ... = 1 / A(x). - Michael Somos, Apr 20 2017
a(0) = 1, a(n) = (2/n)*Sum_{k=1..n} A002129(k)*a(n-k) for n > 0. - Seiichi Manyama, May 06 2017
From Paul D. Hanna, Aug 10 2019: (Start)
G.f.: Sum_{n>=0} x^n * Sum_{k=0..n} binomial(n,k) * (x^(2*n+1) - x^(2*k))^(n-k) = Sum_{n>=0} a(n)*x^(2*n).
G.f.: Sum_{n>=0} x^n * Sum_{k=0..n} binomial(n,k) * (x^(2*n+1) + x^(2*k))^(n-k) * (-1)^k = Sum_{n>=0} a(n)*x^(2*n). (End)
From Peter Bala, Jan 05 2021: (Start)
G.f.: Sum_{n = -oo..oo} x^(4*n^2+2*n) * (1 + x^(4*n+1))/(1 - x^(4*n+1)). See Agarwal, p. 285, equation 6.20 with i = j = 1 and mu = 4.
For prime p of the form 4*k + 3, a(n*p^2 + (p^2 - 1)/4) = a(n).
If n > 0 and p are coprime then a(n*p + (p^2 - 1)/4) = 0. The proofs are similar to those given for the corresponding results for A115110. Cf. A000729.
For prime p of the form 4*k + 1 and for n not congruent to (p - 1)/4 (mod p) we have a(n*p^2 + (p^2 - 1)/4) = 3*a(n) (since b(n), where b(4*n+1) = a(n), is multiplicative). (End)
From Peter Bala, Mar 22 2021: (Start)
G.f. A(q) satisfies:
A(q^2) = Sum_{n = -oo..oo} q^n/(1 - q^(4*n+2)) (set z = q, alpha = q^2, mu = 4 in Agarwal, equation 6.15).
A(q^2) = Sum_{n = -oo..oo} q^(2*n)/(1 - q^(4*n+1)) (set z = q^2, alpha = q, mu = 4 in Agarwal, equation 6.15).
A(q^2) = Sum_{n = -oo..oo} q^n/(1 + q^(2*n+1))^2 = Sum_{n = -oo..oo} q^(3*n+1)/(1 + q^(2*n+1))^2. (End)
G.f.: Sum_{k>=0} a(k) * q^k = Sum_{k>=0} (-1)^k * q^(k*(k+1)) + 2 * Sum_{n>=1, k>=0} (-1)^k * q^(k*(k+2*n+1)+n). - Mamuka Jibladze, May 17 2021
G.f.: Sum_{k>=0} a(k) * q^k = Sum_{k>=0} (-1)^k * q^(k*(k+1)) * (1 + q^(2*k+1))/(1 - q^(2*k+1)). - Mamuka Jibladze, Jun 06 2021
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/2 (A019669). - Amiram Eldar, Oct 15 2022

Extensions

More terms and information from Michael Somos, Mar 23 2003

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

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

Original entry on oeis.org

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

Views

Author

Michael Somos, Jan 15 2012

Keywords

Comments

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

Examples

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

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q] EllipticTheta[ 4, 0, q^4], {q, 0, n}]; (* Michael Somos, Jun 10 2015 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^5 / (eta(x + A)^2 * eta(x^8 + A)), n))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, n==0, A = factor(n); 2 * (-1)^(n%8==5) * prod( k=1, matsize(A)[1], [p, e] = A[k,]; if( p==2, 2 * (e>2) * (-1)^(e<4), p%4==1, e+1, !(e%2))))};

Formula

Expansion of eta(q^2)^5 / (eta(q)^2 * eta(q^8)) in powers of q.
Euler transform of period 8 sequence [ 2, -3, 2, -3, 2, -3, 2, -2, ...].
G.f.: Product_{k>0} (1 - x^(2*k))^5 / ((1 - x^k)^2 * (1 - x^(8*k))).
a(4*n + 2) = a(4*n + 3) = a(8*n + 4) = 0. a(8*n) = A104794(n). a(4*n + 1) = 2 * A134343(n).
a(n) = (-1)^n * A246950(n). a(8*n + 1) = 2 * A113407(n). a(8*n + 5) = -4 * A053692(n). - Michael Somos, Jun 10 2015

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

Original entry on oeis.org

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

Views

Author

Michael Somos, Sep 08 2014

Keywords

Comments

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

Examples

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

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(64), 1), 85); A[1] - 2*A[2] + 4*A[6] - 4*A[9] - 2*A[10] + 4*A[14] + 4*A[17] - 4*A[18] - 6*A[26] + 4*A[30] - 4*A[35] + 4*A[36]; /* Michael Somos, Jun 21 2015 */
  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 4, 0, q] EllipticTheta[ 4, 0, q^4], {q, 0, n}];
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 * eta(x^4 + A)^2 / (eta(x^2 + A) * eta(x^8 + A)), n))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, n==0, A = factor(n); 2 * (-1)^(n%8==1) * prod( k=1, matsize(A)[1], [p, e] = A[k,]; if( p==2, 2 * (e>2) * (-1)^(e<4), p%4==1, (e+1), !(e%2))))};
    

Formula

Expansion of f(-q, -q) * f(q, -q) in powers of q where f(,) is the Ramanujan general theta function. - Michael Somos, Jun 21 2015
Expansion of eta(q)^2 * eta(q^4)^2 / (eta(q^2) * eta(q^8)) in powers of q.
Euler transform of period 8 sequence [ -2, -1, -2, -3, -2, -1, -2, -2, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (64 t)) = 64 (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A053692.
G.f.: Product_{k>0} (1 - x^k)^2 * (1 + x^(2*k)) / (1 + x^(4*k)).
a(n) = (-1)^n * A204531(n). a(4*n + 2) = a(4*n + 3) = a(8*n + 4) = 0.
a(8*n) = A104794(n). a(4*n + 1) = - 2 * A134343(n).
a(8*n + 1) = -2 * A113407(n). a(8*n + 5) = 4 * A053692(n). - Michael Somos, Jun 10 2015

A319307 Expansion of theta_4(q)^16 in powers of q = exp(Pi i t).

Original entry on oeis.org

1, -32, 480, -4480, 29152, -140736, 525952, -1580800, 3994080, -8945824, 18626112, -36714624, 67978880, -118156480, 197120256, -321692928, 509145568, -772845120, 1143441760, -1681379200, 2428524096, -3392205824, 4658843520, -6411152640, 8705492608, -11488092896
Offset: 0

Views

Author

Seiichi Manyama, Sep 16 2018

Keywords

Crossrefs

theta_4(q)^b: A002448 (b=1), A104794 (b=2), A213384 (b=3), A096727 (b=4), A035016 (b=8), A286346 (b=12), this sequence (b=16), A319308 (b=20), A319309 (b=24), A319310 (b=28).
Cf. A000152.

Formula

Expansion of eta(q)^32 / eta(q^2)^16 in powers of q.

A319308 Expansion of theta_4(q)^20 in powers of q = exp(Pi i t).

Original entry on oeis.org

1, -40, 760, -9120, 77560, -497648, 2508000, -10232640, 34729720, -100906760, 259114704, -606957280, 1327461600, -2738111280, 5341699520, -9915552192, 17701924600, -30615844560, 51294999960, -83279292960, 131880275664, -204949382400, 312126610080, -464844224960, 680432137440
Offset: 0

Views

Author

Seiichi Manyama, Sep 16 2018

Keywords

Crossrefs

theta_4(q)^b: A002448 (b=1), A104794 (b=2), A213384 (b=3), A096727 (b=4), A035016 (b=8), A286346 (b=12), A319307 (b=16), this sequence (b=20), A319309 (b=24), A319310 (b=28).

Formula

Expansion of eta(q)^40 / eta(q^2)^20 in powers of q.

A319309 Expansion of theta_4(q)^24 in powers of q = exp(Pi i t).

Original entry on oeis.org

1, -48, 1104, -16192, 170064, -1362336, 8662720, -44981376, 195082320, -721175536, 2319457632, -6631997376, 17231109824, -41469483552, 93703589760, -200343312768, 407488018512, -793229226336, 1487286966928, -2697825744960, 4744779429216, -8110465650176
Offset: 0

Views

Author

Seiichi Manyama, Sep 16 2018

Keywords

Crossrefs

theta_4(q)^b: A002448 (b=1), A104794 (b=2), A213384 (b=3), A096727 (b=4), A035016 (b=8), A286346 (b=12), A319307 (b=16), A319308 (b=20), this sequence (b=24), A319310 (b=28).
Cf. A000156.

Formula

Expansion of eta(q)^48 / eta(q^2)^24 in powers of q.

A319310 Expansion of theta_4(q)^28 in powers of q = exp(Pi i t).

Original entry on oeis.org

1, -56, 1512, -26208, 327656, -3147984, 24189984, -152867520, 811401192, -3681079640, 14500933104, -50376047904, 156797510688, -444306558864, 1163495873088, -2851049839680, 6597606440936, -14512424533488, 30505974273096, -61591664700384, 119983597365744, -226303038736128
Offset: 0

Views

Author

Seiichi Manyama, Sep 16 2018

Keywords

Crossrefs

theta_4(q)^b: A002448 (b=1), A104794 (b=2), A213384 (b=3), A096727 (b=4), A035016 (b=8), A286346 (b=12), A319307 (b=16), A319308 (b=20), A319309 (b=24), this sequence (b=28).

Formula

Expansion of eta(q)^56 / eta(q^2)^28 in powers of q.

A227695 Expansion of psi(x)^2 * phi(-x)^6 in powers of x where phi(), psi() are Ramanujan theta functions.

Original entry on oeis.org

1, -10, 37, -50, -30, 128, -25, -34, -320, 310, 410, -370, -87, -410, 320, 30, 500, 384, -630, -640, -359, 300, -326, 2560, -110, -1098, -1280, -370, 1490, -1850, 269, 1500, 1216, 640, 570, -3328, 340, -2010, -1110, 1790, 768, 3200, 303, 750, -1600, -442
Offset: 0

Views

Author

Michael Somos, Sep 02 2013

Keywords

Comments

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

Examples

			G.f. = 1 - 10*x + 37*x^2 - 50*x^3 - 30*x^4 + 128*x^5 - 25*x^6 - 34*x^7 - 320*x^8 + ...
G.f. = q - 10*q^5 + 37*q^9 - 50*q^13 - 30*q^17 + 128*q^21 - 25*q^25 - 34*q^29 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ (QPochhammer[ x]^5 / QPochhammer[ x^2])^2, {x, 0, n}];
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x + A)^5 / eta(x^2 + A))^2, n))};

Formula

Expansion of q^(-1/4) * (eta(q)^5 / eta(q^2))^2 in powers of q.
Expansion of phi(-x)^5 * f(-x^2)^3 = phi(-x)^2 * f(-x)^6 in powers of x where phi(), f() are Ramanujan theta functions.
Euler transform of period 2 sequence [ -10, -8, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (32 t)) = 8192 (t / i)^4 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A227317.
G.f.: (Product_{k>0} (1 - x^k)^5 / (1 - x^(2*k)))^2.
Convolution of A000729 and A104794.
Showing 1-10 of 13 results. Next