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.

A002895 Domb numbers: number of 2n-step polygons on diamond lattice.

Original entry on oeis.org

1, 4, 28, 256, 2716, 31504, 387136, 4951552, 65218204, 878536624, 12046924528, 167595457792, 2359613230144, 33557651538688, 481365424895488, 6956365106016256, 101181938814289564, 1480129751586116848, 21761706991570726096, 321401321741959062016
Offset: 0

Views

Author

Keywords

Comments

a(n) is the (2n)th moment of the distance from the origin of a 4-step random walk in the plane. - Peter M.W. Gill (peter.gill(AT)nott.ac.uk), Mar 03 2004
Row sums of the cube of A008459. - Peter Bala, Mar 05 2013
Conjecture: Let D(n) be the (n+1) X (n+1) Hankel-type determinant with (i,j)-entry equal to a(i+j) for all i,j = 0..n. Then the number D(n)/12^n is always a positive odd integer. - Zhi-Wei Sun, Aug 14 2013
It appears that the expansions exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + 4*x + 22*x^2 + 152*x^3 + 1241*x^4 + ... and exp( Sum_{n >= 1} 1/4*a(n)*x^n/n ) = 1 + x + 4*x^2 + 25*x^3 + 199*x^4 + ... have integer coefficients. See A267219. - Peter Bala, Jan 12 2016
This is one of the Apéry-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017
Named after the British-Israeli theoretical physicist Cyril Domb (1920-2012). - Amiram Eldar, Mar 20 2021

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

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.)
For primes that do not divide the terms of the sequences A000172, A005258, A002893, A081085, A006077, A093388, A125143, A229111, A002895, A290575, A290576, A005259 see A260793, A291275-A291284 and A133370 respectively.

Programs

  • Maple
    A002895 := n -> add(binomial(n,k)^2*binomial(2*n-2*k,n-k)*binomial(2*k,k), k=0..n): seq(A002895(n), n=0..25); # Wesley Ivan Hurt, Dec 20 2015
    A002895 := n -> binomial(2*n,n)*hypergeom([1/2, -n, -n, -n],[1, 1, 1/2 - n], 1):
    seq(simplify(A002895(n)), n=0..19); # Peter Luschny, May 23 2017
  • Mathematica
    Table[Sum[Binomial[n,k]^2 Binomial[2n-2k,n-k]Binomial[2k,k],{k,0,n}], {n,0,30}] (* Harvey P. Dale, Aug 15 2011 *)
    a[n_] = Binomial[2*n, n]*HypergeometricPFQ[{1/2, -n, -n, -n}, {1, 1, 1/2-n}, 1]; (* or *) a[n_] := SeriesCoefficient[BesselI[0, 2*Sqrt[x]]^4, {x, 0, n}]*n!^2; Table[a[n], {n, 0, 19}] (* Jean-François Alcover, Dec 30 2013, after Vladeta Jovovic *)
    max = 19; Total /@ MatrixPower[Table[Binomial[n, k]^2, {n, 0, max}, {k, 0, max}], 3] (* Jean-François Alcover, Mar 24 2015, after Peter Bala *)
  • PARI
    C=binomial;
    a(n) = sum(k=0,n, C(n,k)^2 * C(2*n-2*k,n-k) * C(2*k,k) );
    /* Joerg Arndt, Apr 19 2013 */
    
  • Python
    from math import comb
    def A002895(n): return (sum(comb(n,k)**2*comb(n-k<<1,n-k)*comb(m:=k<<1,k) for k in range(n+1>>1))<<1) + (0 if n&1 else comb(n,n>>1)**4) # Chai Wah Wu, Jun 17 2025

Formula

a(n) = Sum_{k=0..n} binomial(n, k)^2 * binomial(2n-2k, n-k) * binomial(2k, k).
D-finite with recurrence: n^3*a(n) = 2*(2*n-1)*(5*n^2-5*n+2)*a(n-1) - 64*(n-1)^3*a(n-2). - Vladeta Jovovic, Jul 16 2004
Sum_{n>=0} a(n)*x^n/n!^2 = BesselI(0, 2*sqrt(x))^4. - Vladeta Jovovic, Aug 01 2006
G.f.: hypergeom([1/6, 1/3],[1],108*x^2/(1-4*x)^3)^2/(1-4*x). - Mark van Hoeij, Oct 29 2011
From Zhi-Wei Sun, Mar 20 2013: (Start)
Via the Zeilberger algorithm, Zhi-Wei Sun proved that:
(1) 4^n*a(n) = Sum_{k = 0..n} (binomial(2k,k)*binomial(2(n-k),n-k))^3/ binomial(n,k)^2,
(2) a(n) = Sum_{k = 0..n} (-1)^(n-k)*binomial(n,k)*binomial(2k,n)*binomial(2k,k)* binomial(2(n-k),n-k). (End)
a(n) ~ 2^(4*n+1)/((Pi*n)^(3/2)). - Vaclav Kotesovec, Aug 20 2013
G.f. y=A(x) satisfies: 0 = x^2*(4*x - 1)*(16*x - 1)*y''' + 3*x*(128*x^2 - 30*x + 1)*y'' + (448*x^2 - 68*x + 1)*y' + 4*(16*x - 1)*y. - Gheorghe Coserea, Jun 26 2018
a(n) = Sum_{p+q+r+s=n} (n!/(p!*q!*r!*s!))^2 with p,q,r,s >= 0. See Verrill, p. 5. - Peter Bala, Jan 06 2020
From Peter Bala, Jul 25 2024: (Start)
a(n) = 2*Sum_{k = 1..n} (k/n)*binomial(n, k)^2*binomial(2*n-2*k, n-k)* binomial(2*k, k) for n >= 1.
a(n-1) = (1/2)*Sum_{k = 1..n} (k/n)^3*binomial(n, k)^2*binomial(2*n-2*k, n-k)* binomial(2*k, k) for n >= 1. Cf. A081085. (End)

Extensions

More terms from Vladeta Jovovic, Mar 11 2003