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.

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

Original entry on oeis.org

1, -5, 6, 8, -23, 12, 14, -30, 18, 20, -40, 24, 31, -77, 30, 32, -60, 48, 38, -70, 42, 44, -138, 48, 57, -90, 54, 72, -100, 60, 62, -184, 84, 68, -120, 72, 74, -155, 96, 80, -239, 84, 108, -150, 90, 112, -160, 120, 98, -276, 102, 104, -240, 108, 110, -190, 114
Offset: 0

Views

Author

Michael Somos, Oct 06 2007

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).

Examples

			G.f. = 1 - 5*x + 6*x^2 + 8*x^3 - 23*x^4 + 12*x^5 + 14*x^6 - 30*x^7 + 18*x^8 + ...
G.f. = q - 5*q^3 + 6*q^5 + 8*q^7 - 23*q^9 + 12*q^11 + 14*q^13 - 30*q^15 + ...
		

Crossrefs

Programs

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

Formula

Expansion of psi(x)^4 - 9 * x * psi(x^3)^4 in powers of x where psi() is a Ramanujan theta function.
Expansion of x^(-1/2) * (b(x)^3 * c(x^2)^2 / (3 * c(x)))^(1/2) in powers of x where b(), c() are cubic AGM functions.
Expansion of q^(-1/2) * eta(q)^5 * eta(q^6)^3 / (eta(q^2) * eta(q^3)^3) in powers of q.
Euler transform of period 6 sequence [-5, -4, -2, -4, -5, -4, ...].
a(n) = b(2*n+1) where b() is multiplicative with b(2^e) = 0^e, b(3^e) = 4 - 3^(e+1), b(p^e) = (p^(e+1) - 1)/(p - 1) if p>5.
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = 18 (t/i)^2 g(t) where q = exp(2 Pi i t) and g(t) is the g.f. for A124449
G.f.: Product_{k>0} (1 - x^k)^2 * (1 - x^(2*k))^2 * (1 - x^k + x^(2*k))^3.
G.f.: Sum_{k>0} k * f(x^k) - 9 * k * f(x^(3*k)) where f(x) = x * (1 - x) / ((1 + x) * (1 + x^2)).
G.f.: f(x) - 3 * f(x^2) - 9 * f(x^3) + 2 * f(x^4) + 27 * f(x^6) - 18 * f(x^12) where f() is the g.f. of A000203.
a(n) = A131944(2*n + 1) = A118271(2*n + 1). a(3*n + 2) = 6 * A098098(n).

A185717 Expansion of q^(-1) * c(q^2) * (c(q) - c(q^4)) / 9 in powers of q^2 where c() is a cubic AGM theta function.

Original entry on oeis.org

1, 3, 6, 8, 9, 12, 14, 18, 18, 20, 24, 24, 31, 27, 30, 32, 36, 48, 38, 42, 42, 44, 54, 48, 57, 54, 54, 72, 60, 60, 62, 72, 84, 68, 72, 72, 74, 93, 96, 80, 81, 84, 108, 90, 90, 112, 96, 120, 98, 108, 102, 104, 144, 108, 110, 114, 114, 144, 126, 144, 133, 126, 156, 128
Offset: 0

Views

Author

Michael Somos, Feb 10 2011

Keywords

Comments

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

Examples

			1 + 3*x + 6*x^2 + 8*x^3 + 9*x^4 + 12*x^5 + 14*x^6 + 18*x^7 + 18*x^8 + ...
q + 3*q^3 + 6*q^5 + 8*q^7 + 9*q^9 + 12*q^11 + 14*q^13 + 18*q^15 + ...
		

Crossrefs

Programs

  • Mathematica
    A185717[n_] := SeriesCoefficient[(QPochhammer[q^3, q^3]/QPochhammer[-q^3, q^3])^4*(1/(QPochhammer[q, q^2]*QPochhammer[q^3, q^6])^3), {q, 0, n}];
    Table[A185717[n], {n, 0, 50}] (* G. C. Greubel, Jul 10 2017 *)
  • PARI
    {a(n) = if( n<0, 0, n = 2*n + 1; sumdiv( n, d, if( (n/d) % 3, 1, 0) * d))}
    
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^3 * eta(x^3 + A)^5 / (eta(x + A)^3 * eta(x^6 + A)), n))}

Formula

Expansion of phi(-x^3)^4 / (chi(-x) * chi(-x^3))^3 in powers of x where phi(), chi() are Ramanujan theta functions.
Euler transform of period 6 sequence [ 3, 0, -2, 0, 3, -4, ...].
a(n) = b(2*n + 1) where b(n) is multiplicative with b(2^e) = 0^e, b(3^e) = 3^e, b(p^e) = (p^(e+1) - 1) / (p - 1) if p>3.
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = 2 (t/i)^2 g(t) where q = exp(2 Pi i t) and g() is g.f. for A118271.
a(3*n + 1) = 3 * a(n). A078708(2*n + 1) = A121443(2*n + 1) = A124449(2*n + 1).
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^2/9 = 1.0966227... (A100044). - Amiram Eldar, Dec 28 2023

A214361 Expansion of c(q^2) * (c(q) + 2 * c(q^4)) / 9 in powers of q where c() is a cubic AGM theta function.

Original entry on oeis.org

1, 3, 3, 3, 6, 9, 8, 3, 9, 18, 12, 9, 14, 24, 18, 3, 18, 27, 20, 18, 24, 36, 24, 9, 31, 42, 27, 24, 30, 54, 32, 3, 36, 54, 48, 27, 38, 60, 42, 18, 42, 72, 44, 36, 54, 72, 48, 9, 57, 93, 54, 42, 54, 81, 72, 24, 60, 90, 60, 54, 62, 96, 72, 3, 84, 108, 68, 54
Offset: 1

Views

Author

Michael Somos, Jul 18 2012

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
b(n) = 6*a(n) is the number of solutions in integers to n = x^2 + y^2 + z^2 + w^2 where x + y + z is not divisible by 3. - Michael Somos, Jun 23 2018

Examples

			G.f. = q + 3*q^2 + 3*q^3 + 3*q^4 + 6*q^5 + 9*q^6 + 8*q^7 + 3*q^8 + 9*q^9 + 18*q^10 + ...
a(1) = 1, b(1) = 6 with solutions (w, x, y, z) = {(0, 0, 1, 0), (0, 1, 0, 0), (1, 0, 0, 0)} and their negatives. - _Michael Somos_, Jun 23 2018
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 3, 0, q]^4 - EllipticTheta[ 3, 0, q^3]^4) / 8, {q, 0, n}];
    a[ n_] := SeriesCoefficient[ (QPochhammer[ -q, q^2] QPochhammer[ -q^3, q^6])^3 EllipticTheta[ 2, 0, (-q)^(3/2)]^4 / (-16 (-q)^(1/2)), {q, 0, n}];
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x^2 + A)^6 * eta(x^3 + A) * eta(x^6 + A)^2 * eta(x^12 + A) / ( eta(x + A) * eta(x^4 + A))^3, n))};
    
  • PARI
    {a(n) = if( n<1, 0, sigma(n) + if( n%3==0, -1 * sigma(n/3)) + if( n%4==0, -4 * sigma(n/4)) + if( n%12==0, +4 * sigma(n/12)))};
    
  • 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, 3, p==3, 3^e, (p^(e+1) - 1) / (p - 1))))};

Formula

Expansion of (phi(q)^4 - phi(q^3)^4) / 8 = q * phi(q^3) * (chi(q) * psi(-q^3))^3 = q * psi(-q^3)^4 * (chi(q) * chi(q^3))^3 in powers of q where phi(), psi(), chi() are Ramanujan theta functions.
Expansion of eta(q^2)^6 * eta(q^3) * eta(q^6)^2 * eta(q^12) / ( eta(q) * eta(q^4))^3 in powers of q.
Euler transform of period 12 sequence [3, -3, 2, 0, 3, -6, 3, 0, 2, -3, 3, -4, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = 4 (t/i)^2 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A214456.
a(n) is multiplicative with a(2^e) = 3 if e>0, a(3^e) = 3^e, a(p^e) = (p^(e+1) - 1) / (p - 1) if p>3.
G.f.: x * Product_{k>0} (1 + (-x)^(3*k)) * (1 - x^(6*k))^4 / ( 1 + (-x)^k)^3.
a(n) = -(-1)^n * A124449(n). a(3*n) = 3*a(n). a(2*n) = a(4*n) = 3 * A121443(n). a(2*n + 1) = A185717(n).
From Amiram Eldar, Sep 12 2023: (Start)
Dirichlet g.f.: (1 - 1/2^(2*s-2)) * (1 - 1/3^s) * zeta(s-1) * zeta(s).
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^2/18 = 0.548311... (A086463). (End)
Showing 1-3 of 3 results.