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.

Previous Showing 21-30 of 67 results. Next

A036917 G.f.: (4/Pi^2)*EllipticK(4*x^(1/2))^2.

Original entry on oeis.org

1, 8, 88, 1088, 14296, 195008, 2728384, 38879744, 561787864, 8206324928, 120929313088, 1794924383744, 26802975999424, 402298219288064, 6064992788397568, 91786654611673088, 1393772628452578264, 21227503080738294464, 324160111169327247424
Offset: 0

Views

Author

Keywords

Examples

			G.f. = 1 + 8*x + 88*x^2 +  1088*x^3 + 14296*x^5 + 195008*x^5 + ... - _Michael Somos_, May 29 2023
		

References

  • M. Petkovsek et al., "A=B", Peters, p. ix of second printing.

Crossrefs

The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692,A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Haskell
    a036917 n = sum $ map
       (\k -> (a007318 (2*n-2*k) (n-k))^2 * (a007318 (2*k) k)^2) [0..n]
    -- Reinhard Zumkeller, May 24 2012
    
  • Mathematica
    a[n_] := (16 (n - 1/2)(2*n^2 - 2*n + 1)a[n - 1] - 256(n - 1)^3 a[n - 2])/n^3; a[0] = 1; a[1] = 8; Array[a, 19, 0] (* Or *)
    f[n_] := Sum[(Binomial[2 (n - k), n - k] Binomial[2 k, k])^2, {k, 0, n}]; Array[f, 19, 0] (* Or *)
    lmt = 20; Take[ 4^Range[0, 2 lmt]*CoefficientList[ Series[(4/Pi^2) EllipticK[4 x^(1/2)]^2, {x, 0, lmt}], x^(1/2)], lmt] (* Robert G. Wilson v *)
    a[n_] := HypergeometricPFQ[{1/2, 1/2, -n, -n}, {1, 1/2-n, 1/2-n}, 1] * 4^n * (2n-1)!!^2 / n!^2 (* Vladimir Reshetnikov, Mar 08 2014 *)
    a[ n_] := SeriesCoefficient[ EllipticTheta[3, 0, EllipticNomeQ[16*x]]^4, {x, 0, n}]; (* Michael Somos, May 30 2023 *)
  • PARI
    for(n=0,25, print1(sum(k=0,n, (binomial(2*n-2*k,n-k) *binomial(2*k,k))^2), ", ")) \\ G. C. Greubel, Oct 24 2017
    
  • PARI
    a(n) = if(n<0, 0, polcoeff(agm(1, sqrt(1 - 16*x + x*O(x^n)))^-2, n)); /* Michael Somos, May 29 2023 */

Formula

a(n) = (16*(n-1/2)*(2*n^2-2*n+1)*a(n-1)-256*(n-1)^3*a(n-2))/n^3.
a(n) = Sum_{k=0..n} (C(2 * (n-k), n-k) * C(2 * k, k))^2. [corrected by Tito Piezas III, Oct 19 2010]
a(n) = hypergeom([1/2, 1/2, -n, -n], [1, 1/2-n, 1/2-n], 1) * 4^n * (2n-1)!!^2 / n!^2. - Vladimir Reshetnikov, Mar 08 2014
a(n) ~ 2^(4*n+1) * log(n) / (n*Pi^2) * (1 + (4*log(2) + gamma)/log(n)), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Nov 28 2015
G.f. y=A(x) satisfies: 0 = x^2*(16*x - 1)^2*y''' + 3*x*(16*x - 1)*(32*x - 1)*y'' + (1792*x^2 - 112*x + 1)*y' + 8*(32*x - 1)*y. - Gheorghe Coserea, Jul 03 2018
G.f.: 1 / AGM(1, sqrt(1 - 16*x))^2. - Vaclav Kotesovec, Oct 01 2019
It appears that a(n) is equal to the coefficient of (x*y*z*t)^n in the expansion of (1+x+y+z-t)^n * (1+x+y-z+t)^n * (1+x-y+z+t)^n * (1-x+y+z+t)^n. Cf. A000172. - Peter Bala, Sep 21 2021
G.f. y = A(x) satisfies 0 = x*(1 - 16*x)*(2*y''*y - y'*y') + 2*(1 - 32*x)*y*y' - 16*y*y. - Michael Somos, May 29 2023
Expansion of theta_3(0, q)^4 in powers of m/16 where the modulus m = k^2. - Michael Somos, May 30 2023
From Paul D. Hanna, Mar 25 2024: (Start)
G.f. ( Sum_{n>=0} binomial(2*n,n)^2 * x^n )^2.
G.f. Sum_{n>=0} binomial(2*n,n)^3 * x^n * (1 - 16*x)^n. (End)

Extensions

Replaced complicated definition via a formula with simple generating function provided by Vladeta Jovovic, Dec 01 2003. Thanks to Paul D. Hanna for suggesting this. - N. J. A. Sloane, Mar 25 2024

A143583 Apéry-like numbers: a(n) = (1/C(2n,n))*Sum_{k=0..n} C(2k,k)*C(4k,2k)*C(2n-2k,n-k)*C(4n-4k,2n-2k).

Original entry on oeis.org

1, 12, 164, 2352, 34596, 516912, 7806224, 118803648, 1818757924, 27972399792, 431824158864, 6686855325888, 103814819552016, 1615296581684928, 25180747436810304, 393189646497706752, 6148451986328464164, 96269310864931432368, 1509065592479205772304
Offset: 0

Views

Author

Peter Bala, Aug 25 2008

Keywords

Comments

These numbers bear some analogy to the Apéry numbers A005258. They appear in the evaluation of the spectral zeta function of the non-commutative harmonic oscillator zeta_Q(s) at s = 2 and satisfy a recurrence relation similar to the one satisfied by the Apéry numbers.

Examples

			G.f. = 1 + 12*x + 164*x^2 + 2352*x^3 + 34596*x^4 + 516912*x^5 + ...
		

Crossrefs

Cf. A005258.
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Maple
    a := n -> 1/binomial(2*n, n)*add(binomial(2*k, k)*binomial(4*k, 2*k)*binomial(2*n-2*k, n-k)*binomial(4*n-4*k, 2*n-2*k), k = 0..n): seq(a(n), n = 0..25);
    series( 2*EllipticK(4*x^(1/2))/(Pi*sqrt(1-16*x)), x=0, 20); # Mark van Hoeij, Apr 06 2013
    A143583 := n -> 16^n*hypergeom([1/2, 1/2, -n], [1, 1], 1):
    seq(simplify(A143583(n)), n = 0..18); # Peter Luschny, Nov 12 2022
  • Mathematica
    Table[1/Binomial[2*n,n]*Sum[Binomial[2*k,k]*Binomial[4*k,2*k]*Binomial[2*n-2*k,n-k]*Binomial[4*n-4*k,2*n-2*k],{k,0,n}],{n,0,20}] (* Vaclav Kotesovec, Oct 11 2013 *)

Formula

a(n) = (1/C(2n,n))*sum {k = 0..n} C(2k,k)*C(4k,2k)*C(2n-2k,n-k)*C(4n-4k,2n-2k).
Recurrence relation:
a(0) = 1, a(1) = 12, n^2*a(n) = 4*(8*n^2-8*n+3)*a(n-1) - 256*(n-1)^2*a(n-2).
Congruences:
For odd prime p, a(m*p^r) = a(m*p^(r-1)) (mod p^r) for any m,r in N.
a(n) ~ 16^n/(Pi*sqrt(Pi*n)) * (log(n) + gamma + 6*log(2)), where gamma is the Euler-Mascheroni constant (A001620). - Vaclav Kotesovec, Oct 11 2013
a(n) = sum {k = 0..n} 4^(n-k) C(2k,k)^2*C(2n-2k,n-k). - Tito Piezas III, Dec 12 2014
a(n) = hypergeom([1/2,1/2,n+1],[1,n+3/2],1)*2^(5*n+1)*n!/((2*n+1)!!*Pi) - G. A. Edgar, Dec 10 2016
a(n) = binomial(4*n,2*n)*hypergeom([1/4,3/4,-n,-n], [1,1/4-n,3/4-n], 1). - Peter Luschny, May 14 2020
From Peter Luschny, Nov 12 2022: (Start)
a(n) = 16^n*Sum_{k=0..n} (-1)^k*binomial(-1/2, k)^2*binomial(n, k).
a(n) = 16^n*hypergeom([1/2, 1/2, -n], [1, 1], 1). (End)

A226535 Expansion of b(-q) in powers of q where b() is a cubic AGM theta function.

Original entry on oeis.org

1, 3, 0, -6, -3, 0, 0, 6, 0, -6, 0, 0, 6, 6, 0, 0, -3, 0, 0, 6, 0, -12, 0, 0, 0, 3, 0, -6, -6, 0, 0, 6, 0, 0, 0, 0, 6, 6, 0, -12, 0, 0, 0, 6, 0, 0, 0, 0, 6, 9, 0, 0, -6, 0, 0, 0, 0, -12, 0, 0, 0, 6, 0, -12, -3, 0, 0, 6, 0, 0, 0, 0, 0, 6, 0, -6, -6, 0, 0, 6, 0
Offset: 0

Views

Author

Michael Somos, Sep 22 2013

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Zagier (2009) denotes the g.f. as f(z) in Case B which is associated with F(t) the g.f. of A006077.

Examples

			G.f. = 1 + 3*q - 6*q^3 - 3*q^4 + 6*q^7 - 6*q^9 + 6*q^12 + 6*q^13 - 3*q^16 + ...
		

References

  • D. Zagier, Integral solutions of Apery-like recurrence equations, in: Groups and Symmetries: from Neolithic Scots to John McKay, CRM Proc. Lecture Notes 47, Amer. Math. Soc., Providence, RI, 2009, pp. 349-366.

Crossrefs

The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

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

Formula

Expansion of f(q)^3 / f(q^3) in powers of q where f() is a Ramanujan theta function.
Expansion of 2*b(q^4) - b(q) = b(q^2)^3 / (b(q) * b(q^4)) in powers of q where b() is a cubic AGM theta function.
Expansion of eta(q^2)^9 * eta(q^3) * eta(q^12) / (eta(q) * eta(q^4) * eta(q^6))^3 in powers of q.
Euler transform of period 12 sequence [ 3, -6, 2, -3, 3, -4, 3, -3, 2, -6, 3, -2, ...].
Moebius transform is period 36 sequence [ 3, -3, -9, -3, -3, 9, 3, 3, 0, 3, -3, 9, 3, -3, 9, -3, -3, 0, 3, 3, -9, 3, -3, -9, 3, -3, 0, -3, -3, -9, 3, 3, 9, 3, -3, 0, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (36 t)) = 972^(1/2) (t / i) g(t) where q = exp(2 Pi i t) and g() is the g.f. of A227696.
G.f.: f(q) = F(t(q)) where F() is the g.f. of A006077 and t() is the g.f. of A227454.
G.f.: Product_{k>0} (1 - (-x)^k)^3 / (1 - (-x)^(3*k)).
a(3*n + 2) = a(4*n + 2) = 0.
a(n) = (-1)^n * A005928(n) = (-1)^(((n+1) mod 6 ) > 3) * A113062(n). A113062(n) = |a(n)|.
a(3*n) = A180318(n). a(2*n + 1) = 3 * A123530(n). a(4*n) = A005928(n).

A227216 Expansion of f(-q^2, -q^3)^5 / f(-q)^3 in powers of q where f() is a Ramanujan theta function.

Original entry on oeis.org

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

Views

Author

Michael Somos, Sep 21 2013

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Zagier (2009) refers to Case D corresponding to the Apery numbers (A005258).

Examples

			G.f. = 1 + 3*q + 4*q^2 + 2*q^3 + q^4 + 3*q^5 + 6*q^6 + 4*q^7 - q^9 + ...
		

References

  • D. Zagier, Integral solutions of Apery-like recurrence equations, in: Groups and Symmetries: from Neolithic Scots to John McKay, CRM Proc. Lecture Notes 47, Amer. Math. Soc., Providence, RI, 2009, pp. 349-366.

Crossrefs

The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(5), 1), 20); A[1] + 3*A[2]; /* Michael Somos, Jun 10 2014 */
  • Mathematica
    a[ n_] := If[ n < 1, Boole[ n == 0], Sum[ Re[(3 - I) {1, I, -I, -1, 0}[[ Mod[ d, 5, 1] ]] ], {d, Divisors @ n}]];
    a[ n_] := SeriesCoefficient[ QPochhammer[ q]^2 / (QPochhammer[ q, q^5] QPochhammer[ q^4, q^5])^5, {q, 0, n}]; (* Michael Somos, Jun 10 2014 *)
  • PARI
    {a(n) = if( n<1, n==0, sumdiv(n, d, real( (3 - I) * [ 0, 1, I, -I, -1][ d%5 + 1])))};
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( prod(k=1, n, (1 - x^k)^[ 2, -3, 2, 2, -3][k%5 + 1], 1 + x * O(x^n)), n))};
    
  • Sage
    A = ModularForms( Gamma1(5), 1, prec=20) . basis(); A[0] + 3*A[1]; # Michael Somos, Jun 10 2014
    

Formula

Expansion of f(-q)^2 * (f(-q^5) / f(-q, -q^4))^5 = f(-q^2, -q^3)^2 * (f(-q^5) / f(-q, -q^4))^3 in powers of q where f() is a Ramanujan theta function.
Euler transform of period 5 sequence [ 3, -2, -2, 3, -2, ...].
Moebius transform is period 5 sequence [ 3, 1, -1, -3, 0, ...]. - Michael Somos, Jun 10 2014
G.f. = g(t(q)) where g(), t() are the g.f. for A005258 and A078905.
G.f.: (Product_{k>0} (1 - x^k)^2) / (Product_{k>0} (1 - x^(5*k - 1)) * (1 - x^(5*k - 4)))^5.

A143413 Apéry-like numbers for the constant e: a(n) = 1/(n-1)!*Sum_{k = 0..n+1} (-1)^k*C(n+1,k)*(2*n-k)! for n >= 1.

Original entry on oeis.org

-1, 1, 11, 181, 3539, 81901, 2203319, 67741129, 2346167879, 90449857081, 3843107102339, 178468044946621, 8994348275804891, 488964835817842021, 28523735794360301039, 1777328098986754744081, 117817961601577138782479, 8279178465722546926265329
Offset: 0

Views

Author

Peter Bala, Aug 14 2008

Keywords

Comments

This sequence satisfies the recursion (n-1)^2*a(n) - n^2*a(n-2) = (2*n-1) *(2*n^2 - 2*n+1)*a(n-1), which leads to a rapidly converging series for Napier's constant: e = 2 * Sum_{n >= 1} (-1)^n * n^2/(a(n)* a(n-1)).
Notice the striking parallels with the theory of the Apéry numbers A(n) = A005258(n), which satisfy a similar recurrence relation n^2*A(n) - (n-1)^2*A(n-2) = (11*n^2-11*n+3)*A(n-1) and which appear in the series acceleration formula zeta(2) = 5*Sum {n >= 1} 1/(n^2*A(n)*A(n-1)) = 5*[1/(1*3) + 1/(2^2*3*19) + 1/(3^2*19*147) + ...].

Crossrefs

The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692,A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Maple
    a := n -> 1/(n-1)!*add((-1)^k*binomial(n+1,k)*(2*n-k)!, k = 0..n+1):
    seq(a(n), n = 1..19);
    # Alternative
    a := n -> `if`(n<2, 2*n-1, (2*n)!/(n-1)!*hypergeom([-n-1], [-2*n], -1)):
    seq(simplify(a(n)), n=0..17); # Peter Luschny, Nov 14 2018
  • Mathematica
    Join[{-1}, Table[(1/(n-1)!)*Sum[(-1)^k*Binomial[n+1,k]*(2*n-k)!, {k, 0, n+1}], {n, 1, 50}]] (* G. C. Greubel, Oct 24 2017 *)
  • PARI
    concat([-1], for(n=1,25, print1((1/(n-1)!)*sum(k=0,n+1, (-1)^k*binomial(n+1,k)*(2*n-k)!), ", "))) \\ G. C. Greubel, Oct 24 2017

Formula

a(0):= -1, a(n) = 1/(n-1)!*sum {k = 0..n+1} (-1)^k*C(n+1,k)*(2*n-k)! for n >= 1.
Apart from the initial term, this sequence is the second superdiagonal of the square array A060475; equivalently, the second subdiagonal of the square array A086764.
Recurrence relation: a(0) = -1, a(1) = 1, (n-1)^2*a(n) - n^2*a(n-2) = (2*n-1)*(2*n^2-2*n+1)*a(n-1), n >= 2.
Let b(n) denote the solution to this recurrence with initial conditions b(0) = 0, b(1) = 2. Then b(n) = A143414(n) = 1/(n-1)!*sum {k = 0..n-1} C(n-1,k)*(2*n-k)!. The rational number b(n)/a(n) is equal to the Padé approximation to exp(x) of degree (n-1,n+1) evaluated at x = 1 and b(n)/a(n) -> e very rapidly.
For example, b(100)/a(100) - e is approximately 1.934 * 10^(-436). The identity b(n)*a(n-1) - b(n-1)*a(n) = (-1)^n *2*n^2 leads to rapidly converging series for e and 1/e: e = 2 * Sum_{n >= 1} (-1)^n * n^2/(a(n)*a(n-1)) = 2*[1 + 2^2/(1*11) - 3^2/(11*181) + 4^2/(181*3539) - ...]; 1/e = 1/2 - 2*Sum_{n >= 2} (-1)^n * n^2/(b(n)*b(n-1)) = 1/2 - 2*[2^2/(2*30) - 3^2/(30*492) + 4^2/(492*9620) - ...].
Conjectural congruences: for r >= 0 and odd prime p, calculation suggests that a(p^r*(p+1)) + a(p^r) == 0 (mod p^(r+1)).
a(n) = ((2*n)!/(n-1)!)*hypergeom([-n-1], [-2*n], -1) for n >= 2. - Peter Luschny, Nov 14 2018
a(n) ~ 2^(2*n + 1/2) * n^(n+1) / exp(n + 1/2). - Vaclav Kotesovec, Jul 11 2021

A214262 Expansion of eta(q)^5 * eta(q^3) * eta(q^6)^4 / eta(q^2)^4 in powers of q.

Original entry on oeis.org

1, -5, 9, -11, 24, -45, 50, -53, 81, -120, 120, -99, 170, -250, 216, -203, 288, -405, 362, -264, 450, -600, 528, -477, 601, -850, 729, -550, 840, -1080, 962, -821, 1080, -1440, 1200, -891, 1370, -1810, 1530, -1272, 1680, -2250, 1850, -1320, 1944, -2640, 2208
Offset: 1

Views

Author

Michael Somos, Jul 09 2012

Keywords

Comments

Zagier (2009) writes "... associated to the weight 3 Eisenstein series g(z) = Sigma b(n)q^n = q - 5q^2 + 9q^3 - 11q^4 + ...".
Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).

Examples

			G.f. = q - 5*q^2 + 9*q^3 - 11*q^4 + 24*q^5 - 45*q^6 + 50*q^7 - 53*q^8 + 81*q^9 + ...
		

References

  • D. Zagier, Integral solutions of Apery-like recurrence equations, in: Groups and Symmetries: from Neolithic Scots to John McKay, CRM Proc. Lecture Notes 47, Amer. Math. Soc., Providence, RI, 2009, pp. 349-366.

Crossrefs

Cf. A111661.
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Mathematica
    a[ n_] := If[ n < 1, 0, DivisorSum[ n, -(-1)^# #^2 JacobiSymbol[ -3, n/#] &]]; (* Michael Somos, Oct 06 2013 *)
    a[ n_] := SeriesCoefficient[ q QPochhammer[ q]^5 QPochhammer[ q^3] QPochhammer[ q^6]^4 / QPochhammer[ q^2]^4, {q, 0, n}]; (* Michael Somos, Oct 06 2013 *)
  • PARI
    {a(n) = if( n<1, 0, sumdiv( n, d, -(-1)^d * d^2 * kronecker( -3, n/d)))}; /* Michael Somos, Oct 06 2013 */
    
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x + A)^5 * eta(x^3 + A) * eta(x^6 + A)^4 / eta(x^2 + A)^4, n))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<0, 0, A = factor( n); prod( k=1, matsize(A)[1], if(p = A[k,1], e = A[k,2]; if( p==3, 9^e, if( p==2, -(4^(e+1) + 9*(-1)^(e+1)) / 5, if( p%6==1, ((p^2)^(e+1) - 1) / (p^2 - 1), ((p^2)^(e+1) - (-1)^(e+1)) / (p^2 + 1)))))))};

Formula

Expansion of (1/9) * c(q) * b(q)^2 * c(q^2) / b(q^2) = (c(q)^3 - 8*c(q^2)^3) / 27 in powers of q where b(), c() are cubic AGM theta functions.
Euler transform of period 6 sequence [ -5, -1, -6, -1, -5, -6, ...]. - Michael Somos, Oct 06 2013
a(n) is multiplicative with a(3^e) = 9^e, a(2^e) = -(4^(e+1) + 9*(-1)^(e+1)) / 5, a(p^e) = ((p^2)^(e+1) - 1) / (p^2 - 1) if p == 1 (mod 6), a(p^e) = ((p^2)^(e+1) - (-1)^(e+1)) / (p^2 + 1) if p == 5 (mod 6).
G.f. is a period 1 Fourier series which satisfies f(-1 / (6 t)) = 192^(1/2) (t/i)^3 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A111661.
G.f.: Sum_{k>0} -(-1)^k * k^2 * x^k / (1 + x^k + x^(2*k)) = Sum_{k>0} Kronecker( -3, k) * (x^k - x^(2*k)) / (1 + x^k)^3.

A260832 a(n) = numerator(Jtilde2(n)).

Original entry on oeis.org

1, 3, 41, 147, 8649, 32307, 487889, 1856307, 454689481, 1748274987, 26989009929, 104482114467, 6488426222001, 25239009088827, 393449178700161, 1535897056631667, 1537112996582116041, 6016831929058214523, 94316599529950360769, 369994845516850143483, 23244865440911268112681
Offset: 0

Views

Author

Michel Marcus, Nov 17 2015

Keywords

Comments

Jtilde2(n) are Apéry-like rational numbers that arise in the calculation of zetaQ(2), the spectral zeta function for the non-commutative harmonic oscillator using a Gaussian hypergeometric function.

Crossrefs

Cf. A056982 (denominators), A013661 (zeta(2)), A264541 (Jtilde3).
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Maple
    a := n -> numer(simplify(hypergeom([1/2, 1/2, -n], [1, 1], 1))):
    seq(a(n), n = 0..20); # Peter Luschny, Dec 08 2022
  • Mathematica
    Numerator[Table[Sum[ (-1)^k*Binomial[-1/2, k]^2*Binomial[n, k], {k, 0, n}], {n,0,50}]] (* G. C. Greubel, Feb 15 2017 *)
  • PARI
    a(n) = numerator(sum(k=0, n, (-1)^k*binomial(-1/2,k)^2*binomial(n, k)));
    
  • PARI
    a(n) = numerator(sum(k=0, n, binomial(2*k, k)*binomial(4*k, 2*k)* binomial(2*(n-k),n-k)*binomial(4*(n-k),2*(n-k))) / (2^(4*n)* binomial(2*n,n)));

Formula

Jtilde2(n) = J2(n)/J2(0) with J2(0) = 3*zeta(2) (normalization).
And 4n^2*J2(n) - (8n^2-8n+3)*J2(n-1) + 4(n-1)^2*J2(n-2) = 0 with J2(0) = 3*zeta(2) and J2(1) = 9*zeta(2)/4.
Jtilde2(n) = Sum_{k=0..n} (-1)^k*binomial(-1/2,k)^2*binomial(n,k).
Jtilde2(n) = Sum_{k=0..n} binomial(2*k,k)*binomial(4*k,2*k)*binomial(2*(n-k),n-k)*binomial(4*(n-k),2*(n-k))/(2^(4*n)*binomial(2*n,n)).
From Andrey Zabolotskiy, Oct 04 2016 and Dec 08 2022: (Start)
Jtilde2(n) = Integral_{ x >= 0 } (L_n(x))^2*exp(-x)/sqrt(Pi*x) dx, where L_n(x) is the Laguerre polynomial (A021009).
G.f. of Jtilde2(n): 2F1(1/2,1/2;1;z/(z-1))/(1-z).
Jtilde2(n) = A143583(n) / 16^n. (End)
a(n) = numerator(hypergeom([1/2, 1/2, -n], [1, 1], 1)). - Peter Luschny, Dec 08 2022

A279619 Expansion of g.f. of A002652 in powers of the g.f. of A279618.

Original entry on oeis.org

1, 2, 22, 336, 6006, 117348, 2428272, 52303680, 1160427510, 26337699740, 608642155660, 14272471122560, 338764038330480, 8123136091556640, 196484811079765440, 4788469475873867520, 117465323079289162230, 2898183118626011393100
Offset: 1

Views

Author

Lynette O'Brien, Dec 15 2016

Keywords

Comments

G.f. is the square root of the g.f. for A183204.
This sequence is c_n in Theorem 6.1 in O'Brien's thesis.
Also see Conjecture 5.4 in Chan, Cooper and Sica's paper.

Examples

			G.f. = 1 + 2*x + 22*x^2 + 336*x^3 + 6006*x^4 + ....
		

References

  • L. O'Brien, Modular forms and two new integer sequences at level 7, Massey University, 2016.

Crossrefs

The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Magma
    I:=[2, 22]; [1] cat [n le 2 select I[n] else ((26*n^2-39*n+15)* Self(n-1) + 3*(3*n-4)*(3*n-5)*Self(n-2))/n^2 : n in [1..50]] // G. C. Greubel, Jul 04 2018
  • Mathematica
    RecurrenceTable[{a[n+1] == ((26*n^2+13*n+2)*a[n] + 3*(3*n-1)*(3*n-2)*a[n-1])/ (n + 1)^2, a[-1] == 0, a[0] == 1}, a, {n, 0, 50}] (* G. C. Greubel, Jul 04 2018 *)
    CoefficientList[Series[Sqrt[7]*(1/(25 - 80*x + 24*Sqrt[1 - 27*x]*Sqrt[1+x]))^(1/4) * Hypergeometric2F1[1/12, 5/12, 1, 13824*x^7/(1 - 21*x + 8*x^2 + Sqrt[1 - 27*x] * (1 - 8*x)*Sqrt[1+x])^3], {x, 0, 20}], x] (* Vaclav Kotesovec, Jul 04 2018 *)

Formula

(n+1)^2*a_7(n+1) = (26*n^2+13*n+2)*a_7(n) + 3*(3*n-1)*(3*n-2)*a_7(n-1), a(0)=1, a(-1)=0.
Conjecture: For any positive integer n and any prime p with p equiv. 0,1,2 or 4 modulo 7, a(n) equiv. a(n)=a(n_0)a(n_1)...a(n_r) modulo p, where n=n_0+n_1p+...n_rp^r is the base p representation of n.
Conjecture: a(n)~ C n^(-3/2) 27^n where C=0.0955223052681267146513079107870296256727946666510071798669948234917659...

A227454 Expansion of q * (f(q^9) / f(q))^3 in powers of q where f() is a Ramanujan theta function.

Original entry on oeis.org

1, -3, 9, -22, 51, -108, 221, -429, 810, -1476, 2631, -4572, 7802, -13056, 21519, -34918, 55935, -88452, 138332, -213990, 327852, -497592, 748833, -1117692, 1655719, -2434938, 3556791, -5161808, 7445631, -10677096, 15226658, -21599469, 30485268, -42817788
Offset: 1

Views

Author

Michael Somos, Sep 22 2013

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Zagier (2009) denotes the g.f. as t(z) in Case B which is associated with F(t) the g.f. of A006077.

Examples

			G.f. = q - 3*q^2 + 9*q^3 - 22*q^4 + 51*q^5 - 108*q^6 + 221*q^7 - 429*q^8 + ...
		

References

  • D. Zagier, Integral solutions of Apery-like recurrence equations, in: Groups and Symmetries: from Neolithic Scots to John McKay, CRM Proc. Lecture Notes 47, Amer. Math. Soc., Providence, RI, 2009, pp. 349-366.

Crossrefs

The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ q (QPochhammer[ -q^9] / QPochhammer[ -q])^3, {q, 0, n}]
  • PARI
    {a(n) = local(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( (eta(x + A) * eta(x^4 + A) * eta(x^18 + A)^3 / (eta(x^2 + A)^3 * eta(x^9 + A) * eta(x^36 + A)))^3, n))}

Formula

Expansion of c(-q^3) / (-3 * b(-q)) in powers of q where b(), c() are cubic AGM theta functions.
Expansion of (eta(q) * eta(q^4) * eta(q^18)^3 / (eta(q^2)^3 * eta(q^9) * eta(q^36)))^3 in powers of q.
Euler transform of period 36 sequence [ -3, 6, -3, 3, -3, 6, -3, 3, 0, 6, -3, 3, -3, 6, -3, 3, -3, 0, -3, 3, -3, 6, -3, 3, -3, 6, 0, 3, -3, 6, -3, 3, -3, 6, -3, 0, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (36 t)) = (1/27) g(t) where q = exp(2 Pi i t) and g() is the g.f. of A227498.
G.f. t(q) satisfies f(q) = F(t(q)) where F() is the g.f. of A006077 and f() is the g.f. of A226535
G.f.: x * (Product_{k>0} (1 - (-x)^(9*k)) / (1 - (-x)^k))^3.
a(n) = -(-1)^n * A121589(n).

A264541 a(n) = numerator(Jtilde3(n)).

Original entry on oeis.org

0, 1, 65, 13247, 704707, 660278641, 357852111131, 309349386395887, 240498440880062263, 148443546307725010253, 61760947097005048531, 13658972396318235617977, 723464275788899734058353751, 489812222050789870424202126629, 2614176630672654770175367214389, 204702102697072009862200307064701369
Offset: 0

Views

Author

Michel Marcus, Nov 17 2015

Keywords

Comments

Jtilde3(n) are Apéry-like rational numbers that arise in the calculation of zetaQ(3), the spectral zeta function for the non-commutative harmonic oscillator using a Gaussian hypergeometric function.

Crossrefs

Cf. A002117 (zeta(3)), A260832 (Jtilde2), A264542 (denominators).
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Mathematica
    Numerator[Table[-2*Sum[(-1)^k*Binomial[-1/2, k]^2*Binomial[n, k]*Sum[ 1/(Binomial[-1/2, j]^2*(2*j + 1)^3), {j, 0, k - 1}], {k, 0, n}], {n, 0, 50}]] (* G. C. Greubel, Oct 24 2017 *)
  • PARI
    a(n) = numerator(-2*sum(k=0, n, (-1)^k*binomial(-1/2, k)^2*binomial(n, k)*sum(j=0, k-1, 1/(binomial(-1/2,j)^2*(2*j+1)^3))));

Formula

Jtilde3(n) = J3(n) - J3(0)*Jtilde2(n) (normalization).
4n^2*J3(n) - (8n^2-8n+3)*J3(n-1) + 4(n-1)^2*J3(n-2) = 2^n*(n-1)!/(2n-1)!! with J3(0)=7*zeta(3) and J3(1)=21*zeta(3)/4 + 1/2.
Previous Showing 21-30 of 67 results. Next