A081085 -id:A081085 - cp's OEIS frontend

A063007 T(n,k) = binomial(n,k)*binomial(n+k,k), 0 <= k <= n, triangle read by rows.

1, 1, 2, 1, 6, 6, 1, 12, 30, 20, 1, 20, 90, 140, 70, 1, 30, 210, 560, 630, 252, 1, 42, 420, 1680, 3150, 2772, 924, 1, 56, 756, 4200, 11550, 16632, 12012, 3432, 1, 72, 1260, 9240, 34650, 72072, 84084, 51480, 12870, 1, 90, 1980, 18480, 90090, 252252, 420420, 411840, 218790, 48620

Offset: 0

Henry Bottomley, Jul 02 2001

T(n,k) is the number of compatible k-sets of cluster variables in Fomin and Zelevinsky's Cluster algebra of finite type B_n. Take a row of this triangle regarded as a polynomial in x and rewrite as a polynomial in y := x+1. The coefficients of the polynomial in y give a row of triangle A008459 (squares of binomial coefficients). For example, x^2+6*x+6 = y^2+4*y+1. - Paul Boddington, Mar 07 2003
T(n,k) is the number of lattice paths from (0,0) to (n,n) using steps E=(1,0), N=(0,1) and D=(1,1) (i.e., bilateral Schroeder paths), having k N=(0,1) steps. E.g. T(2,0)=1 because we have DD; T(2,1) = 6 because we have NED, NDE, EDN, END, DEN and DNE; T(2,2)=6 because we have NNEE, NENE, NEEN, EENN, ENEN and ENNE. - Emeric Deutsch, Apr 20 2004
Another version of [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...] DELTA [0, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...] = 1; 1, 0; 1, 2, 0; 1, 6, 6, 0; 1, 12, 30, 20, 0; ..., where DELTA is the operator defined in A084938. - Philippe Deléham Apr 15 2005
Terms in row n are the coefficients of the Legendre polynomial P(n,2x+1) with increasing powers of x.
From Peter Bala, Oct 28 2008: (Start)
Row n of this triangle is the f-vector of the simplicial complex dual to an associahedron of type B_n (a cyclohedron) [Fomin & Reading, p.60]. See A008459 for the corresponding h-vectors for associahedra of type B_n and A001263 and A033282 respectively for the h-vectors and f-vectors for associahedra of type A_n.
An alternative description of this triangle in terms of f-vectors is as follows. Let A_n be the root lattice generated as a monoid by {e_i - e_j: 0 <= i,j <= n+1}. Let P(A_n) be the polytope formed by the convex hull of this generating set. Then the rows of this array are the f-vectors of a unimodular triangulation of P(A_n) [Ardila et al.]. A008459 is the corresponding array of h-vectors for these type A_n polytopes. See A127674 (without the signs) for the array of f-vectors for type C_n polytopes and A108556 for the array of f-vectors associated with type D_n polytopes.
The S-transform on the ring of polynomials is the linear transformation of polynomials that is defined on the basis monomials x^k by S(x^k) = binomial(x,k) = x(x-1)...(x-k+1)/k!. Let P_n(x) denote the S-transform of the n-th row polynomial of this array. In the notation of [Hetyei] these are the Stirling polynomials of the type B associahedra. The first few values are P_1(x) = 2*x + 1, P_2(x) = 3*x^2 + 3*x + 1 and P_3(x) = (10*x^3 + 15*x^2 + 11*x + 3)/3. These polynomials have their zeros on the vertical line Re x = -1/2 in the complex plane, that is, the polynomials P_n(-x) satisfy a Riemann hypothesis. See A142995 for further details. The sequence of values P_n(k) for k = 0,1,2,3, ... produces the n-th row of A108625. (End)
This is the row reversed version of triangle A104684. - Wolfdieter Lang, Sep 12 2016
T(n, k) is also the number of (n-k)-dimensional faces of a convex n-dimensional Lipschitz polytope of real functions f defined on the set X = {1, 2, ..., n+1} which satisfy the condition f(n+1) = 0 (see Gordon and Petrov). - Stefano Spezia, Sep 25 2021
The rows seem to give (up to sign) the coefficients in the expansion of the integer-valued polynomial ((x+1)*(x+2)*(x+3)*...*(x+n) / n!)^2 in the basis made of the binomial(x+i,i). - F. Chapoton, Oct 09 2022
Chapoton's observation above is correct: the precise expansion is ((x+1)*(x+2)*(x+3)*...*(x+n)/ n!)^2 = Sum_{k = 0..n} (-1)^k*T(n,n-k)*binomial(x+2*n-k, 2*n-k), as can be verified using the WZ algorithm. For example, n = 3 gives ((x+1)*(x+2)*(x+3)/3!)^2 = 20*binomial(x+6,6) - 30*binomial(x+5,5) + 12*binomial(x+4,4) - binomial(x+3,3). - Peter Bala, Jun 24 2023

			The triangle T(n, k) starts:
  n\k 0  1    2     3     4      5      6      7      8     9
  0:  1
  1:  1  2
  2:  1  6    6
  3:  1 12   30    20
  4:  1 20   90   140    70
  5:  1 30  210   560   630    252
  6:  1 42  420  1680  3150   2772    924
  7:  1 56  756  4200 11550  16632  12012   3432
  8:  1 72 1260  9240 34650  72072  84084  51480  12870
  9:  1 90 1980 18480 90090 252252 420420 411840 218790 48620
... reformatted by _Wolfdieter Lang_, Sep 12 2016
From _Petros Hadjicostas_, Jul 11 2020: (Start)
Its inverse (from Table II, p. 92, in Ser's book) is
   1;
  -1/2,  1/2;
   1/3, -1/2,    1/6;
  -1/4,  9/20,  -1/4,   1/20;
   1/5, -2/5,    2/7,  -1/10,  1/70;
  -1/6,  5/14, -25/84,  5/36, -1/28,  1/252;
   1/7, -9/28,  25/84, -1/6,   9/154, -1/84, 1/924;
   ... (End)

J. M. Borwein and P. B. Borwein, Pi and the AGM, Wiley, 1987, p. 366.
J. Ser, Les Calculs Formels des Séries de Factorielles. Gauthier-Villars, Paris, 1933, Table I, p. 92.
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.

T. D. Noe, Rows n = 0..100 of triangle, flattened
F. Ardila, M. Beck, S. Hosten, J. Pfeifle and K. Seashore, Root polytopes and growth series of root lattices, arXiv:0809.5123 [math.CO], 2008.
Peter Bala, Deformations of the Hadamard product of power series.
Cyril Banderier, Combinatoire analytique des chemins et des cartes, Thesis (2001), page 49.
Paul Barry, Continued fractions and transformations of integer sequences, JIS 12 (2009), #09.7.6.
H. J. Brothers, Pascal's Prism: Supplementary Material.
David Callan, A bijection for Delannoy paths, arXiv:2202.04649 [math.CO], 2022.
F. Chapoton, Enumerative properties of generalized associahedra, Séminaire Lotharingien de Combinatoire, B51b (2004), 16 pp.
Johann Cigler, Some remarks and conjectures related to lattice paths in strips along the x-axis, arXiv:1501.04750 [math.CO], 2015-2016.
Mark Dukes and Chris D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.
Mark Dukes and Chris D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, Electronic Journal Of Combinatorics, 23(1) (2016), #P1.45.
S. Fomin and N. Reading, Root systems and generalized associahedra, Lecture notes for IAS/Park-City 2004; arXiv:math/0505518 [math.CO], 2005-2008.
S. Fomin and A. Zelevinsky, Cluster algebras I: Foundations, J. Amer. Math. Soc. 15(2) (2002), 497-529.
S. Fomin and A. Zelevinsky, Y-systems and generalized associahedra, Ann. of Math. (2) 158 (2003), no. 3, 977-1018.
J. Gordon and F. Petrov, Combinatorics of the Lipschitz Polytope, Arnold Mathematical Journal (2016).
G. Hetyei, Face enumeration using generalized binomial coefficients. This is the draft version of Hetyei's paper referenced below. [Archived version]
Gabor Hetyei, The Stirling polynomial of a simplicial complex Discrete and Computational Geometry 35(3) (2006), 437-455.
Hsien-Kuei Hwang and Satoshi Kuriki, Integrated empirical measures and generalizations of classical goodness-of-fit statistics, arXiv:2404.06040 [math.ST], 2024. See p. 11.
C. Lanczos, Applied Analysis (Annotated scans of selected pages). See page 514.
T. Manneville and V. Pilaud, Compatibility fans for graphical nested complexes, arXiv preprint arXiv:1501.07152 [math.CO], 2015.
Thomas Selig, Combinatorial aspects of sandpile models on wheel and fan graphs, arXiv:2202.06487 [math.CO], 2022.
J. Ser, Les Calculs Formels des Séries de Factorielles, Gauthier-Villars, Paris, 1933 [Local copy].
J. Ser, Les Calculs Formels des Séries de Factorielles (Annotated scans of some selected pages) See Table I, page 92.
V. Strehl, Recurrences and Legendre transform, Séminaire Lotharingien de Combinatoire, B29b (1992), 22 pp.
R. A. Sulanke, Objects counted by the central Delannoy numbers., J. Integer Seq. 6 (2003), no. 1, Article 03.1.5.
D. Zagier, Integral solutions of Apery-like recurrence equations.

See A331430 for an essentially identical triangle, except with signed entries.
Columns include A000012, A002378, A033487 on the left and A000984, A002457, A002544 on the right.
Main diagonal is A006480.
Row sums are A001850. Alternating row sums are A033999.
Cf. A007318, A008459, A009766, A046899, A104684, A109983, A127674.
Cf. A033282 (f-vectors type A associahedra), A108625, A080721 (f-vectors type D associahedra).
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.)

a063007 n k = a063007_tabl !! n !! k
a063007_row n = a063007_tabl !! n
a063007_tabl = zipWith (zipWith (*)) a007318_tabl a046899_tabl
-- Reinhard Zumkeller, Nov 18 2014

/* As triangle: */ [[Binomial(n,k)*Binomial(n+k,k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Sep 03 2015

p := (n,x) -> orthopoly[P](n,1+2*x): seq(seq(coeff(p(n,x),x,k), k=0..n), n=0..9);

Flatten[Table[Binomial[n, k]Binomial[n + k, k], {n, 0, 10}, {k, 0, n}]] (* Harvey P. Dale, Dec 24 2011 *)
Table[CoefficientList[Hypergeometric2F1[-n, n + 1, 1, -x], x], {n, 0, 9}] // Flatten
(* Peter Luschny, Mar 09 2018 *)

{T(n, k) = local(t); if( n<0, 0, t = (x + x^2)^n; for( k=1, n, t=t'); polcoeff(t, k) / n!)} /* Michael Somos, Dec 19 2002 */

{T(n, k) = binomial(n, k) * binomial(n+k, k)} /* Michael Somos, Sep 22 2013 */

{T(n, k) = if( k<0 || k>n, 0, (n+k)! / (k!^2 * (n-k)!))} /* Michael Somos, Sep 22 2013 */

T(n, k) = (n+k)!/(k!^2*(n-k)!) = T(n-1, k)*(n+k)/(n-k) = T(n, k-1)*(n+k)*(n-k+1)/k^2 = T(n-1, k-1)*(n+k)*(n+k-1)/k^2.
binomial(x, n)^2 = Sum_{k>=0} T(n,k) * binomial(x, n+k). - Michael Somos, May 11 2012
T(n, k) = A109983(n, k+n). - Michael Somos, Sep 22 2013
G.f.: G(t, z) = 1/sqrt(1-2*z-4*t*z+z^2). Row generating polynomials = P_n(1+2z), i.e., T(n, k) = [z^k] P_n(1+2*z), where P_n are the Legendre polynomials. - Emeric Deutsch, Apr 20 2004
Sum_{k>=0} T(n, k)*A000172(k) = Sum_{k>=0} T(n, k)^2 = A005259(n). - Philippe Deléham, Jun 08 2005
1 + z*d/dz(log(G(t,z))) = 1 + (1 + 2*t)*z + (1 + 8*t + 8*t^2)*z^2 + ... is the o.g.f. for a signed version of A127674. - Peter Bala, Sep 02 2015
If R(n,t) denotes the n-th row polynomial then x^3 * exp( Sum_{n >= 1} R(n,t)*x^n/n ) = x^3 + (1 + 2*t)*x^4 + (1 + 5*t + 5*t^2)*x^5 + (1 + 9*t + 21*t^2 + 14*t^3)*x^6 + ... is an o.g.f for A033282. - Peter Bala, Oct 19 2015
P(n,x) := 1/(1 + x)*Integral_{t = 0..x} R(n,t) dt are (modulo differences of offset) the row polynomials of A033282. - Peter Bala, Jun 23 2016
From Peter Bala, Mar 09 2018: (Start)
R(n,x) = Sum_{k = 0..n} binomial(2*k,k)*binomial(n+k,n-k)*x^k.
R(n,x) = Sum_{k = 0..n} binomial(n,k)^2*x^k*(1 + x)^(n-k).
n*R(n,x) = (1 + 2*x)*(2*n - 1)*R(n-1,x) - (n - 1)*R(n-2,x).
R(n,x) = (-1)^n*R(n,-1 - x).
R(n,x) = 1/n! * (d/dx)^n ((x^2 + x)^n). (End)
The row polynomials are R(n,x) = hypergeom([-n, n + 1], [1], -x). - Peter Luschny, Mar 09 2018
T(n,k) = C(n+1,k)*A009766(n,k). - Bob Selcoe, Jan 18 2020 (Connects this triangle with the Catalan triangle. - N. J. A. Sloane, Jan 18 2020)
If we let A(n,k) = (-1)^(n+k)*(2*k+1)*(n*(n-1)*...*(n-(k-1)))/((n+1)*...*(n+(k+1))) for n >= 0 and k = 0..n, and we consider both T(n,k) and A(n,k) as infinite lower triangular arrays, then they are inverses of one another. (Empty products are by definition 1.) See the example below. The rational numbers |A(n,k)| appear in Table II on p. 92 in Ser's (1933) book. - Petros Hadjicostas, Jul 11 2020
From Peter Bala, Nov 28 2021: (Start)
Row polynomial R(n,x) = Sum_{k >= n} binomial(k,n)^2 * x^(k-n)/(1+x)^(k+1) for x > -1/2.
R(n,x) = 1/(1 + x)^(n+1) * hypergeom([n+1, n+1], [1], x/(1 + x)).
R(n,x) = (1 + x)^n * hypergeom([-n, -n], [1], x/(1 + x)).
R(n,x) = hypergeom([(n+1)/2, -n/2], [1], -4*x*(1 + x)).
If we set R(-1,x) = 1, we can run the recurrence n*R(n,x) = (1 + 2*x)*(2*n - 1)*R(n-1,x) - (n - 1)*R(n-2,x) backwards to give R(-n,x) = R(n-1,x).
R(n,x) = [t^n] ( (1 + t)*(1 + x*(1 + t)) )^n. (End)
n*T(n,k) = (2*n-1)*T(n-1,k) + (4*n-2)*T(n-1,k-1) - (n-1)*T(n-2,k). - Fabián Pereyra, Jun 30 2022
From Peter Bala, Oct 07 2024: (Start)
n-th row polynomial R(n,x) = Sum_{k = 0..n} binomial(n, k) * x^k o (1 + x)^(n-k), where o denotes the black diamond product of power series as defined by Dukes and White (see Bala, Section 4.4, exercise 3).
Denote this triangle by T. Then T * transpose(T) = A143007, the square array of crystal ball sequences for the A_n X A_n lattices.
Let S denote the triangle ((-1)^(n+k)*T(n, k))n,k >= 0, a signed version of this triangle. Then S^(-1) * T = A007318, Pascal's triangle; it appears that T * S^(-1) = A110098.
T = A007318 * A115951. (End)

A143007 Square array, read by antidiagonals, where row n equals the crystal ball sequence for the 2*n-dimensional lattice A_n x A_n.

Original entry on oeis.org

1, 1, 1, 1, 5, 1, 1, 13, 13, 1, 1, 25, 73, 25, 1, 1, 41, 253, 253, 41, 1, 1, 61, 661, 1445, 661, 61, 1, 1, 85, 1441, 5741, 5741, 1441, 85, 1, 1, 113, 2773, 17861, 33001, 17861, 2773, 113, 1, 1, 145, 4873, 46705, 142001, 142001, 46705, 4873, 145, 1

Offset: 0

Peter Bala, Jul 22 2008

The A_n lattice consists of all vectors v = (x_1,...,x_(n+1)) in Z^(n+1) such that x_1 + ... + x_(n+1) = 0. The lattice is equipped with the norm ||v|| = 1/2*(|x_1| + ... + |x_(n+1)|). Pairs of lattice points (v,w) in the product lattice A_n x A_n have norm ||(v,w)|| = ||v|| + ||w||. Then the k-th term in the crystal ball sequence for the A_n x A_n lattice gives the number of such pairs (v,w) for which ||(v,w)|| is less than or equal to k.
This array has a remarkable relationship with Apery's constant zeta(3). The row (or column) and main diagonal entries of the array occur in series acceleration formulas for zeta(3). For row n entries there holds zeta(3) = (1+1/2^3+...+1/n^3) + Sum_{k >= 1} 1/(k^3*T(n,k-1)*T(n,k)). Also, as consequence of Apery's proof of the irrationality of zeta(3), we have a series acceleration formula along the main diagonal of the table: zeta(3) = 6 * sum {n >= 1} 1/(n^3*T(n-1,n-1)*T(n,n)). Apery's result appears to generalize to the other diagonals of the table. Calculation suggests the following result may hold: zeta(3) = 1 + 1/2^3 + ... + 1/k^3 + Sum_{n >= 1} (2*n+k)*(3*n^2 +3*n*k +k^2)/(n^3*(n+k)^3*T(n-1,n+k-1)*T(n,n+k)).
For the corresponding results for the constant zeta(2), related to the crystal ball sequences of the lattices A_n, see A108625. For corresponding results for log(2), coming from either the crystal ball sequences of the hypercubic lattices A_1 x ... x A_1 or the lattices of type C_n, see A008288 and A142992 respectively.

			The table begins
n\k|0...1.....2......3.......4.......5
======================================
0..|1...1.....1......1.......1.......1
1..|1...5....13.....25......41......61 A001844
2..|1..13....73....253.....661....1441 A143008
3..|1..25...253...1445....5741...17861 A143009
4..|1..41...661...5741...33001..142001 A143010
5..|1..61..1441..17861..142001..819005 A143011
........
Example row 1 [1,5,13,...]:
The lattice A_1 x A_1 is equivalent to the square lattice of all integer lattice points v = (x,y) in Z x Z equipped with the taxicab norm ||v|| = (|x| + |y|). There are 4 lattice points (marked with a 1 on the figure below) satisfying ||v|| = 1 and 8 lattice points (marked with a 2 on the figure) satisfying ||v|| = 2. Hence the crystal ball sequence for the A_1 x A_1 lattice begins 1, 1+4 = 5, 1+4+8 = 13, ... .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . 2 . . . . .
. . . . 2 1 2 . . . .
. . . 2 1 0 1 2 . . .
. . . . 2 1 2 . . . .
. . . . . 2 . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
Row 1 = [1,5,13,...] is the sequence of partial sums of A008574; row 2 = [1,13,73,...] is the sequence of partial sums of A008530, so row 2 is the crystal ball sequence for the lattice A_2 x A_2 (the 4-dimensional di-isohexagonal orthogonal lattice).
Read as a triangle the array begins
n\k|0...1....2....3...4...5
===========================
0..|1
1..|1...1
2..|1...5....1
3..|1..13...13....1
4..|1..25...73...25...1
5..|1..41..253..253..41...1

G. C. Greubel, Antidiagonals n = 0..50, flattened
R. Bacher, P. de la Harpe and B. Venkov, Séries de croissance et séries d'Ehrhart associées aux réseaux de racines, C. R. Acad. Sci. Paris, 325 (Series 1) (1997), 1137-1142.
J. H. Conway and N. J. A. Sloane, Low dimensional lattices VII Coordination sequences, Proc. R. Soc. Lond., Ser. A, 453 (1997), 2369-2389.
Armin Straub, Multivariate Apéry numbers and supercongruences of rational functions, Algebra & Number Theory, Vol. 8, No. 8 (2014), pp. 1985-2008; arXiv preprint, arXiv:1401.0854 [math.NT], 2014.

Cf. A001844 (row 1), A005259 (main diagonal), A008288, A008530 (first differences of row 2), A008574 (first differences of row 1), A085478, A108625, A142992, A143003, A143004, A143005, A143006, A143008 (row 2), A143009 (row 3), A142010 (row 4), A143011 (row 5).
Cf. A227845 (antidiagonal sums), A246464.
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.)
Cf. A246563, A352653.

A:= func< n,k | (&+[(Binomial(n,j)*Binomial(n+k-j,k-j))^2: j in [0..n]]) >; // Array
A143007:= func< n,k | A(n-k,k) >; // Antidiagonal triangle
[A143007(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 05 2023

with(combinat): T:= (n,k) -> add(binomial(n+j,2*j)*binomial(2*j,j)^2*binomial(k+j,2*j), j = 0..n): for n from 0 to 9 do seq(T(n,k),k = 0..9) end do;

T[n_, k_]:= HypergeometricPFQ[{-k, k+1, -n, n+1}, {1, 1, 1}, 1]; Table[T[n-k, k], {n,0,12}, {k,0,n}]//Flatten (* Jean-François Alcover, Mar 06 2013 *)

/* Print as a square array: */
{T(n, k)=sum(j=0, n, binomial(n+j, 2*j)*binomial(2*j, j)^2*binomial(k+j, 2*j))}
for(n=0, 10, for(k=0,10, print1(T(n,k), ", "));print(""))

/* (1) G.f. A(x,y) when read as a triangle: */
{T(n,k)=local(A=1+x); A=sum(m=0, n, x^m * y^m / (1-x +x*O(x^n))^(2*m+1) * sum(k=0, m, binomial(m, k)^2*x^k)^2 ); polcoeff(polcoeff(A, n,x), k,y)}
for(n=0, 10, for(k=0,n, print1(T(n,k), ", "));print(""))

/* (2) G.f. A(x,y) when read as a triangle: */
{T(n,k)=local(A=1+x); A=sum(m=0, n, x^m/(1-x*y +x*O(x^n))^(2*m+1) * sum(k=0, m, binomial(m, k)^2 * x^k * y^k)^2 ); polcoeff(polcoeff(A, n,x), k,y)}
for(n=0, 10, for(k=0,n, print1(T(n,k), ", "));print(""))

/* (3) G.f. A(x,y) when read as a triangle: */
{T(n,k)=local(A=1+x); A=sum(m=0, n, x^m*sum(k=0, m, binomial(m , k)^2 * y^k * sum(j=0, k, binomial(k, j)^2 * x^j)+x*O(x^n))); polcoeff(polcoeff(A, n,x), k,y)}
for(n=0, 10, for(k=0,n, print1(T(n,k), ", "));print(""))

/* (4) G.f. A(x,y) when read as a triangle: */
{T(n,k)=local(A=1+x); A=sum(m=0, n, x^m*sum(k=0, m, binomial(m, k)^2 * y^(m-k) * sum(j=0, k, binomial(k, j)^2 * x^j * y^j)+x*O(x^n))); polcoeff(polcoeff(A, n,x), k,y)}
for(n=0, 10, for(k=0,n, print1(T(n,k), ", "));print(""))
/* End */

def A(n,k): return sum((binomial(n,j)*binomial(n+k-j,k-j))^2 for j in range(n+1)) # array
def A143007(n,k): return A(n-k,k) # antidiagonal triangle
flatten([[A143007(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Oct 05 2023

T(n,k) = Sum_{j = 0..n} C(n+j,2*j)*C(2*j,j)^2*C(k+j,2*j).
The array is symmetric T(n,k) = T(k,n).
The main diagonal [1,5,73,1445,...] is the sequence of Apery numbers A005259.
The entries in the k-th column satisfy the Apery-like recursion n^3*T(n,k) + (n-1)^3*T(n-2,k) = (2*n-1)*(n^2-n+1+2*k^2+2*k)*T(n-1,k).
The LDU factorization of the square array is L * D * transpose(L), where L is the lower triangular array A085478 and D is the diagonal matrix diag(C(2n,n)^2). O.g.f. for row n: The generating function for the coordination sequence of the lattice A_n is [Sum_{k = 0..n} C(n,k)^2*x^k ]/(1-x)^n. Thus the generating function for the coordination sequence of the product lattice A_n x A_n is {[Sum_{k = 0..n} C(n,k)^2*x^k]/(1-x)^n}^2 and hence the generating function for row n of this array, the crystal ball sequence of the lattice A_n x A_n, equals [Sum_{k = 0..n} C(n,k)^2*x^k]^2/(1-x)^(2n+1) = 1/(1-x)*[Legendre_P(n,(1+x)/(1-x))]^2. See [Conway & Sloane].
Series acceleration formulas for zeta(3): Row n: zeta(3) = (1 + 1/2^3 + ... + 1/n^3) + Sum_{k >= 1} 1/(k^3*T(n,k-1)*T(n,k)), n = 0,1,2,... . For example, the fourth row of the table (n = 3) gives zeta(3) = (1 + 1/2^3 + 1/3^3) + 1/(1^3*1*25) + 1/(2^3*25*253) + 1/(3^3*253*1445) + ... . See A143003 for further details.
Main diagonal: zeta(3) = 6 * Sum_{n >= 1} 1/(n^3*T(n-1,n-1)*T(n,n)). Conjectural result for other diagonals: zeta(3) = 1 + 1/2^3 + ... + 1/k^3 + Sum_{n >= 1} (2*n+k)*(3*n^2+3*n*k+k^2)/(n^3*(n+k)^3*T(n-1,n+k-1)*T(n,n+k)).
Sum_{k=0..n} T(n-k,k) = A227845(n) (antidiagonal sums). - Paul D. Hanna, Aug 27 2014
The main superdiagonal numbers S(n) := T(n,n+1) appear to satisfy the supercongruences S(m*p^r - 1) == S(m*p^(r-1) - 1) (mod p^(3*r)) for prime p >= 5 and m, r in N (this is true: see A352653. - Peter Bala, Apr 16 2022).
From Paul D. Hanna, Aug 27 2014: (Start)
G.f. A(x,y) = Sum_{n>=0, k=0..n} T(n,k)*x^n*y^k can be expressed by:
(1) Sum_{n>=0} x^n * y^n / (1-x)^(2*n+1) * [Sum_{k=0..n} C(n,k)^2 * x^k]^2,
(2) Sum_{n>=0} x^n / (1 - x*y)^(2*n+1) * [Sum_{k=0..n} C(n,k)^2 * x^k * y^k]^2,
(3) Sum_{n>=0} x^n * Sum_{k=0..n} C(n,k)^2 * y^k * Sum_{j=0..k} C(k,j)^2 * x^j,
(4) Sum_{n>=0} x^n * Sum_{k=0..n} C(n,k)^2 * y^(n-k) * Sum_{j=0..k} C(k,j)^2 * x^j * y^j. (End)
From Peter Bala, Jun 23 2023: (Start)
T(n,k) = Sum_{j = 0..n} C(n,j)^2 * C(n+k-j,k-j)^2.
T(n,k) = binomial(n+k,k)^2 * hypergeom([-n, -n, -k, -k],[-n - k, -n - k, 1], 1).
T(n,k) = hypergeom([n+1, -n, k+1, -k], [1, 1, 1], 1). (End)
From Peter Bala, Jun 28 2023: (Start)
T(n,k) = the coefficient of (x*z)^n*(y*t)^k in the expansion of 1/( (1 - x - y)*(1 - z - t) - x*y*z*t ).
T(n,k) = A(n, k, n, k) in the notation of Straub, equation 7.
The supercongruences T(n*p^r, k*p^r) == T(n*p^(r-1), k*p^(r-1)) (mod p^(3*r)) hold for all primes p >= 5 and positive integers n and k.
The formula T(n,k) = hypergeom([n+1, -n, k+1, -k], [1, 1, 1], 1) allows the table indexing to be extended to negative values of n and k; we have T(-n,k) = T(n-1,k) and T(n,-k) = T(n,k-1) leading to T(-n,-k) = T(n-1, k-1). (End)
From G. C. Greubel, Oct 05 2023: (Start)
Let t(n, k) = T(n-k, k) be the antidiagonal triangle, then:
t(n, k) = t(n, n-k).
Sum_{k=0..floor(n/2)} t(n-k,k) = A246563(n).
t(2*n+1, n+1) = A352653(n+1). (End)
From Peter Bala, Sep 27 2024: (Start)
The square array = A063007 * transpose(A063007) (LU factorization).
Let L denote the lower triangular array (l(n,k))n,k >= 0, where l(n, k) = (-1)^(n+k) * binomial(n, k)*binomial(n+k, k). (L is a signed version of A063007 and L = A063007 * A007318 ^(-1).)
Then the square array = L * transpose(A108625).
L^2 * transpose(A108625) = the Hadamard product of A108625 with itself (both identities can be verified using the MulZeil procedure in Doron Zeilberger's MultiZeilberger package to find recurrences for the double sums involved). (End)

Spelling/notation corrections by Charles R Greathouse IV, Mar 18 2010

A006077 (n+1)^2a(n+1) = (9n^2+9n+3)a(n) - 27n^2a(n-1), with a(0) = 1 and a(1) = 3.

Original entry on oeis.org

1, 3, 9, 21, 9, -297, -2421, -12933, -52407, -145293, -35091, 2954097, 25228971, 142080669, 602217261, 1724917221, 283305033, -38852066421, -337425235479, -1938308236731, -8364863310291, -24286959061533, -3011589296289, 574023003011199, 5028616107443691

Offset: 0

N. J. A. Sloane

sign

This is the Taylor expansion of a special point on a curve described by Beauville. - Matthijs Coster, Apr 28 2004
Conjecture: Let W(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. If n == 1 (mod 3) then W(n) = 0. When n == 0 or 2 (mod 3), W(n)*(-1)^(floor((n+1)/3))/6^n is always a positive odd integer. - Zhi-Wei Sun, Aug 21 2013
Conjecture: Let p == 1 (mod 3) be a prime, and write 4*p = x^2 + 27*y^2 with x, y integers and x == 1 (mod 3). Then W(p-1) == (-1)^{(p+1)/2}*(x-p/x) (mod p^2), where W(n) is defined as the above. - Zhi-Wei Sun, Aug 23 2013
This is one of the Apery-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017
Diagonal of rational functions 1/(1 - (x^2*y + y^2*z - z^2*x + 3*x*y*z)), 1/(1 - (x^3 + y^3 - z^3 + 3*x*y*z)), 1/(1 + x^3 + y^3 + z^3 - 3*x*y*z). - Gheorghe Coserea, Aug 04 2018

			G.f. = 1 + 3*x + 9*x^2 + 21*x^3 + 9*x^4 - 297*x^5 - 2421*x^6 - 12933*x^7 - ...

Matthijs Coster, Over 6 families van krommen [On 6 families of curves], Master's Thesis (unpublished), Aug 26 1983.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
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.

Seiichi Manyama, Table of n, a(n) for n = 0..1400 (terms 0..200 from T. D. Noe)
Arnaud Beauville, Les familles stables de courbes sur P_1 admettant quatre fibres singulières, Comptes Rendus, Académie Sciences Paris, no. 294, May 24 1982.
A. Bostan, S. Boukraa, J.-M. Maillard, and J.-A. Weil, Diagonals of rational functions and selected differential Galois groups, arXiv preprint arXiv:1507.03227 [math-ph], 2015.
M. Coster, Email, Nov 1990
Ofir Gorodetsky, New representations for all sporadic Apéry-like sequences, with applications to congruences, arXiv:2102.11839 [math.NT], 2021. See B p. 2.
S. Herfurtner, Elliptic surfaces with four singular fibres, Mathematische Annalen, 1991. Preprint.
Bradley Klee, Checking Weierstrass data, 2023.
Amita Malik and Armin Straub, Divisibility properties of sporadic Apéry-like numbers, Research in Number Theory, 2016, 2:5.
Stéphane Ouvry and Alexios Polychronakos, Lattice walk area combinatorics, some remarkable trigonometric sums and Apéry-like numbers, arXiv:2006.06445 [math-ph], 2020.
Zhi-Wei Sun, Three mysterious conjectures on Hankel-type determinants, a message to Number Theory List, August 22, 2013.
Zhi-Wei Sun, Connections between p = x^2+3*y^2 and Franel numbers, J. Number Theory 133(2013), 2914-2928.

Related to diagonal of rational functions: A268545-A268555.
Cf. A091401.
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.

a := n -> 3^n*hypergeom([-n/3, (1-n)/3, (2-n)/3], [1, 1], 1):
seq(simplify(a(n)), n=0..24); # Peter Luschny, Nov 01 2017

Table[Sum[(-1)^k*3^(n - 3*k)*Binomial[n, 3*k]*Binomial[2*k, k]* Binomial[3*k, k], {k, 0, Floor[n/3]}], {n, 0, 50}] (* G. C. Greubel, Oct 24 2017 *)
a[ n_] := SeriesCoefficient[ HypergeometricPFQ[ {1/3, 2/3}, {1}, x^3 / (x - 1/3)^3 ] / (1 - 3 x), {x, 0, n}]; (* Michael Somos, Nov 01 2017 *)

subst(eta(q)^3/eta(q^3), q, serreverse(eta(q^9)^3/eta(q)^3*q)) \\ (generating function) Helena Verrill (verrill(AT)math.lsu.edu), Apr 20 2009 [for (-1)^n*a(n)]

diag(expr, N=22, var=variables(expr)) = {
  my(a = vector(N));
  for (k = 1, #var, expr = taylor(expr, var[#var - k + 1], N));
  for (n = 1, N, a[n] = expr;
    for (k = 1, #var, a[n] = polcoeff(a[n], n-1)));
  return(a);
};
diag(1/(1 + x^3 + y^3 + z^3 - 3*x*y*z), 25)

seq(N) = {
  my(a = vector(N)); a[1] = 3; a[2] = 9;
  for (n = 2, N-1, a[n+1] = ((9*n^2+9*n+3)*a[n] - 27*n^2*a[n-1])/(n+1)^2);
  concat(1,a);
};
seq(24)
\\ test: y=subst(Ser(seq(202)), 'x, -'x/27); 0 == x*(x^2+9*x+27)*y'' + (3*x^2+18*x+27)*y' + (x+3)*y
\\ Gheorghe Coserea, Nov 09 2017

{a(n) = my(A); if( n<0, 0, A = x * O(x^n); (-1)^n * polcoeff( subst(eta(x + A)^3 / eta(x^3 + A), x, serreverse( x * eta(x^9 + A)^3 / eta(x + A)^3)), n))}; /* Michael Somos, Nov 01 2017 */

G.f.: hypergeom([1/3, 2/3], [1], x^3/(x-1/3)^3) / (1-3*x). - Mark van Hoeij, Oct 25 2011
a(n) = Sum_{k=0..floor(n/3)}(-1)^k*3^(n-3k)*C(n,3k)*C(2k,k)*C(3k,k). - Zhi-Wei Sun, Aug 21 2013
0 = x*(x^2+9*x+27)*y'' + (3*x^2 + 18*x + 27)*y' + (x + 3)*y, where y(x) = A(x/-27). - Gheorghe Coserea, Aug 26 2016
a(n) = 3^n*hypergeom([-n/3, (1-n)/3, (2-n)/3], [1, 1], 1). - Peter Luschny, Nov 01 2017
From Bradley Klee, Jun 05 2023: (Start)
The g.f. T(x) obeys a period-annihilating ODE:
0=3*(-1 + 9*x)*T(x) + (-1 + 9*x)^2*T'(x) +  x*(1 - 9*x + 27*x^2)*T''(x).
The periods ODE can be derived from the following Weierstrass data:
g2 = 3*(-8 + 9*(1 - 9*x)^3)*(1 - 9*x);
g3 = 8 - 36*(1 - 9*x)^3 + 27*(1 - 9*x)^6;
which determine an elliptic surface with four singular fibers. (End)

More terms from Kok Seng Chua (chuaks(AT)ihpc.nus.edu.sg), Jun 20 2000

A093388 (n+1)^2a(n+1) = (17n^2+17n+6)a(n) - 72n^2a(n-1).

Original entry on oeis.org

1, 6, 42, 312, 2394, 18756, 149136, 1199232, 9729882, 79527084, 654089292, 5408896752, 44941609584, 375002110944, 3141107339328, 26402533581312, 222635989516122, 1882882811380284, 15967419789558804, 135752058036988848, 1156869080242393644

Offset: 0

Matthijs Coster, Apr 29 2004

nonn

This is the Taylor expansion of a special point on a curve described by Beauville.

This is one of the Apery-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017

			A(x) = 1 + 6*x + 42*x^2 + 312*x^3 + 2394*x^4 + 18756*x^5 + ... is the g.f.

Matthijs Coster, Over 6 families van krommen [On 6 families of curves], Master's Thesis (unpublished), Aug 26 1983.

Seiichi Manyama, Table of n, a(n) for n = 0..1050 (terms 0..200 from Vincenzo Librandi)
Arnaud Beauville, Les familles stables de courbes sur P_1 admettant quatre fibres singulières, Comptes Rendus, Académie Sciences Paris, no. 294, May 24 1982.
Jonathan M. Borwein, Dirk Nuyens, Armin Straub and James Wan, Some Arithmetic Properties of Short Random Walk Integrals, May 2011.
Shaun Cooper, Apéry-like sequences defined by four-term recurrence relations, arXiv:2302.00757 [math.NT], 2023.
Matthijs Coster, Sequences
Ofir Gorodetsky, New representations for all sporadic Apéry-like sequences, with applications to congruences, arXiv:2102.11839 [math.NT], 2021. See F p. 2.
S. Herfurtner, Elliptic surfaces with four singular fibres, Mathematische Annalen, 1991. Preprint.
Bradley Klee, Checking Weierstrass data, 2023.
Robert S. Maier, On Rationally Parametrized Modular Equations, arXiv:math/0611041 [math.NT], 2006.
Amita Malik and Armin Straub, Divisibility properties of sporadic Apéry-like numbers, Research in Number Theory, 2016, 2:5
Armin Straub, Arithmetic aspects of random walks and methods in definite integration, Ph. D. Dissertation, School Of Science And Engineering, Tulane University, 2012. - From _N. J. A. Sloane_, Dec 16 2012
Zhi-Hong Sun, Congruences for Apéry-like numbers, arXiv:1803.10051 [math.NT], 2018.
Zhi-Hong Sun, New congruences involving Apéry-like numbers, arXiv:2004.07172 [math.NT], 2020.
H. Verrill, Some congruences related to modular forms, Section 2.2.
D. Zagier, Integral solutions of Apery-like recurrence equations. See line F in sporadic solutions table of page 5.

This is the seventh sequence in the family beginning A002894, A006077, A081085, A005258, A000172, A002893.
Cf. A091401.
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.

f:=proc(n) option remember; local m; if n=0 then RETURN(1); fi; if n=1 then RETURN(6); fi; m:=n-1; ((17*m^2+17*m+6)*f(n-1)-72*m^2*f(n-2))/n^2; end;

Table[(-1)^n*Sum[Binomial[n,k]*(-8)^k*Sum[Binomial[n-k,j]^3,{j,0,n-k}],{k,0,n}],{n,0,20}] (* Vaclav Kotesovec, Oct 14 2012 *)

a(n)=(-1)^n*sum(k=0,n,binomial(n,k)*(-8)^k*sum(j=0,n-k,binomial(n-k,j)^3));

seq(N) = {
  my(a = vector(N)); a[1] = 6; a[2] = 42;
  for (n=3, N, a[n] = ((17*n^2 - 17*n + 6)*a[n-1] - 72*(n-1)^2*a[n-2])/n^2);
  concat(1,a);
};
seq(20)  \\ Gheorghe Coserea, Aug 26 2016

a(n) = (-1)^n * Sum_{k=0..n} binomial(n, k) * (-8)^k * Sum_{j=0..n-k} binomial(n-k, j)^3. - Helena Verrill (verrill(AT)math.lsu.edu), Aug 09 2004
G.f.: hypergeom([1/3, 2/3],[1],x^2*(8*x-1)/(2*x-1/3)^3)/(1-6*x). - Mark van Hoeij, Oct 25 2011
a(n) ~ 3^(2*n+3/2)/(Pi*n). - Vaclav Kotesovec, Oct 14 2012
G.f. A(x) satisfies: 0 = x*(x+8)*(x+9)*y'' + (3*x^2 + 34*x + 72)*y' + (x+6)*y, where y(x) = A(-x/72). - Gheorghe Coserea, Aug 26 2016
From Bradley Klee, Jun 05 2023: (Start)
The g.f. T(x) obeys a period-annihilating ODE:
0=6*(-1 + 12*x)*T(x) + (1 - 34*x + 216*x^2)*T'(x) + x*(-1 + 8*x)*(-1 + 9*x)*T''(x).
The periods ODE can be derived from the following Weierstrass data:
g2 = 12*(-1 + 6*x)*(-1 + 18*x - 84*x^2 + 24*x^3);
g3 = -8*(1 - 12*x + 24*x^2)*(-1 + 24*x - 192*x^2 + 504*x^3 + 72*x^4);
which determine an elliptic surface with four singular fibers.  (End)

A125143 Almkvist-Zudilin numbers: Sum_{k=0..n} (-1)^(n-k) * ((3^(n-3k) (3k)!) / (k!)^3) binomial(n,3k) binomial(n+k,k).

Original entry on oeis.org

1, -3, 9, -3, -279, 2997, -19431, 65853, 292329, -7202523, 69363009, -407637387, 702049401, 17222388453, -261933431751, 2181064727997, -10299472204311, -15361051476987, 900537860383569, -10586290198314843, 74892552149042721, -235054958584593843

Offset: 0

R. K. Guy, Jan 11 2007

Apart from signs, this is one of the Apery-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017
Diagonal of rational function 1/(1 - (x + y + z + w - 27*x*y*z*w)). - Gheorghe Coserea, Oct 14 2018
Named after the Swedish mathematician Gert Einar Torsten Almkvist (1934-2018) and the Russian mathematician Wadim Walentinowitsch Zudilin (b. 1970). - Amiram Eldar, Jun 23 2021

G. Almkvist and W. Zudilin, Differential equations, mirror maps and zeta values. In Mirror Symmetry V, N. Yui, S.-T. Yau, and J. D. Lewis (eds.), AMS/IP Studies in Advanced Mathematics 38 (2007), International Press and Amer. Math. Soc., pp. 481-515. Cited in Chan & Verrill.
Helena Verrill, in a talk at the annual meeting of the Amer. Math. Soc., New Orleans, LA, Jan 2007 on "Series for 1/pi".

Seiichi Manyama, Table of n, a(n) for n = 0..1052 (terms 0..200 from Arkadiusz Wesolowski)
Gert Almkvist, Christian Krattenthaler, and Joakim Petersson, Some new formulas for pi, Experiment. Math., Vol. 12 (2003), pp. 441-456. (Math Rev MR2043994 by W. Zudilin)
G. Almkvist and W. Zudilin, Differential equations, mirror maps and zeta values, arXiv:math/0402386 [math.NT], 2004.
Tewodros Amdeberhan, and Roberto Tauraso, Supercongruences for the Almkvist-Zudilin numbers, arXiv:1506.08437 [math.NT], 2015.
Yuliy Baryshnikov, Stephen Melczer, Robin Pemantle and Armin Straub, Diagonal asymptotics for symmetric rational functions via ACSV, LIPIcs Proceedings of Analysis of Algorithms 2018, arXiv:1804.10929 [math.CO], 2018.
Heng Huat Chan and Helena Verrill, The Apery numbers, the Almkvist-Zudilin numbers and new series for 1/Pi, Math. Res. Lett., Vol. 16, No. 3 (2009), pp. 405-420.
Shaun Cooper, Apéry-like sequences defined by four-term recurrence relations, arXiv:2302.00757 [math.NT], 2023. See Table 2 p. 7.
Ofir Gorodetsky, New representations for all sporadic Apéry-like sequences, with applications to congruences, arXiv:2102.11839 [math.NT], 2021. See delta p. 3.
Ji-Cai Liu, A p-adic analogue of Chan and Verrill's formula for 1/Pi, arXiv:2008.06675 [math.NT], 2020.
Ji-Cai Liu, A supercongruence related to Ramanujan-type formula for 1/Pi, Bull. Math. Soc. Sci. Math. Roumanie Tome 67 (115), No. 4, 2024, 483-492. See p. 483.
Amita Malik and Armin Straub, Divisibility properties of sporadic Apéry-like numbers, Research in Number Theory, Vol. 2, (2016), Article 5.
Zhi-Hong Sun, Congruences for Apéry-like numbers, arXiv:1803.10051 [math.NT], 2018.
Zhi-Hong Sun, Super congruences concerning binomial coefficients and Apéry-like numbers, arXiv:2002.12072 [math.NT], 2020.
Zhi-Hong Sun, New congruences involving Apéry-like numbers, arXiv:2004.07172 [math.NT], 2020.

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.

Table[Sum[(-1)^(n-k)*((3^(n-3*k)*(3*k)!)/(k!)^3) *Binomial[n,3*k] *Binomial[n+k,k],{k,0,n}],{n,0,20}] (* Vaclav Kotesovec, Sep 11 2013 *)

a(n) = sum(k=0,n, (-1)^(n-k)*((3^(n-3*k)*(3*k)!)/(k!)^3)*binomial(n,3*k)*binomial(n+k,k) );

a(n) = Sum_{k=0..n} (-1)^(n-k) * ((3^(n-3*k) * (3*k)!) / (k!)^3) * binomial(n,3*k) * binomial(n+k,k) . - Arkadiusz Wesolowski, Jul 13 2011
Recurrence: n^3*a(n) = -(2*n-1)*(7*n^2 - 7*n + 3)*a(n-1) - 81*(n-1)^3*a(n-2). - Vaclav Kotesovec, Sep 11 2013
Lim sup n->infinity |a(n)|^(1/n) = 9. - Vaclav Kotesovec, Sep 11 2013
G.f. y=A(x) satisfies: 0 = x^2*(81*x^2 + 14*x + 1)*y''' + 3*x*(162*x^2 + 21*x + 1)*y'' + (21*x + 1)*(27*x + 1)*y' + 3*(27*x + 1)*y. - Gheorghe Coserea, Oct 15 2018
G.f.: hypergeom([1/8, 5/8], [1], -256*x^3/((81*x^2 + 14*x + 1)*(-x + 1)^2))^2/((81*x^2 + 14*x + 1)^(1/4)*sqrt(-x + 1)). - Sergey Yurkevich, Aug 31 2020

Edited and more terms added by Arkadiusz Wesolowski, Jul 13 2011

A229111 Expansion of the g.f. of A053723 in powers of the g.f. of A121591.

Original entry on oeis.org

1, -5, 35, -275, 2275, -19255, 163925, -1385725, 11483875, -91781375, 688658785, -4581861025, 22550427925, 8852899375, -2431720493125, 47471706909725, -699843878180125, 9141002535744625, -111232778205154375, 1288777160650004375, -14372445132730778975

Offset: 1

Michael Somos, Sep 30 2013

sign

In Verrill (1999) section 2.1, t = (eta(q^5) / eta(q))^6 the g.f. of A121591 and f = eta(q^5)^5 / eta(q) the g.f. of A053723.

Apart from signs, this is one of the Apery-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017

			G.f. = x - 5*x^2 + 35*x^3 - 275*x^4 + 2275*x^5 - 19255*x^6 + 163925*x^7 + ...

Seiichi Manyama, Table of n, a(n) for n = 1..958
Shaun Cooper, Apéry-like sequences defined by four-term recurrence relations, arXiv:2302.00757 [math.NT], 2023. See Table 2 p. 7.
Amita Malik and Armin Straub, Divisibility properties of sporadic Apéry-like numbers, Research in Number Theory, 2016, 2:5.
Ofir Gorodetsky, New representations for all sporadic Apéry-like sequences, with applications to congruences, arXiv:2102.11839 [math.NT], 2021. See eta p. 3.
L. O'Brien, Modular forms and two new integer sequences at level 7, Massey University, 2016.
H. Verrill, Some Congruences related to modular forms, Max Planck Institute, 1999.

Cf. A053723, A109064, A121591.
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.

a[n_] := a[n] = Switch[n, 1, 1, 2, -5, _, (1/(n-1)^3) ((1-2(n-1)) (11(n-2) (n-1)+5) a[n-1] - 125 (n-2)^3 a[n-2])];
a /@ Range[21] (* Jean-François Alcover, Jan 13 2020 *)

{a(n) = my(m = n-1); if( n<1, 0, if( n<3, [1, -5][n], -( (5*(m - 1))^3*a(n-2) + (2*m - 1)*(11*(m^2 - m) +5)*a(n-1) )/ m^3))};

{a(n) = sum(k=0, n-1, (-1)^k*binomial(n-1, k)^3*binomial(5*k-(n-1), 3*(n-1)))} \\ Seiichi Manyama, Sep 02 2020

n^3 * a(n+1) = -(2*n - 1)*(11*n*(n - 1) + 5) * a(n) - 125 * (n - 1)^3 * a(n-1).
a(n*p^k) == (p^3 + Kronecker(p, 5)) * a(n*p^(k-1)) - Kronecker(p, 5) * p^3*a(n*p^(-2)) (mod p^k). [Verrill, 1999]
a(n) = Sum_{k=0..n-1} (-1)^k * binomial(n-1,k)^3 * binomial(5*k-(n-1),3*(n-1)). - Seiichi Manyama, Sep 02 2020

A290575 Apéry-like numbers Sum_{k=0..n} (C(n,k) * C(2*k,n))^2.

Original entry on oeis.org

1, 4, 40, 544, 8536, 145504, 2618176, 48943360, 941244376, 18502137184, 370091343040, 7508629231360, 154145664817600, 3196100636757760, 66834662101834240, 1407913577733228544, 29849617614785770456, 636440695668355742560, 13638210075999240396736, 293565508750164008207104, 6344596821114216520841536

Offset: 0

Hugo Pfoertner, Aug 06 2017

Sequence epsilon in Almkvist, Straten, Zudilin article.

Seiichi Manyama, Table of n, a(n) for n = 0..500
G. Almkvist, D. van Straten, and W. Zudilin, Generalizations of Clausen’s formula and algebraic transformations of Calabi-Yau differential equations, Proc. Edinburgh Math. Soc.54 (2) (2011), 273-295.
Shaun Cooper, Apéry-like sequences defined by four-term recurrence relations, arXiv:2302.00757 [math.NT], 2023.
Ofir Gorodetsky, New representations for all sporadic Apéry-like sequences, with applications to congruences, arXiv:2102.11839 [math.NT], 2021. See epsilon p. 3.
Amita Malik and Armin Straub, Divisibility properties of sporadic Apéry-like numbers, Research in Number Theory, 2016, 2:5
Armin Straub and Wadim Zudilin, Positivity of rational functions and their diagonals, J. Approx. Theory, 195:57-69, 2015. arXiv:1312.3732 [math.NT].
Zhi-Hong Sun, Congruences for Apéry-like numbers, arXiv:1803.10051 [math.NT], 2018.
Zhi-Hong Sun, New congruences involving Apéry-like numbers, arXiv:2004.07172 [math.NT], 2020.
Jeremy Tan, Simplifying a binomial sum for bridge deals with specific voids, Mathematics Stack Exchange, 2024.

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.)

Table[Sum[(Binomial[n, k]*Binomial[2*k, n])^2, {k, 0, n}], {n, 0, 25}] (* G. C. Greubel, Oct 23 2017 *)
a[n_] := Binomial[2 n, n]^2 HypergeometricPFQ[{1/2 - n/2, 1/2 - n/2, -n/2, -n/2}, {1, 1/2 - n, 1/2 - n}, 1];
Table[a[n], {n, 0, 20}] (* Peter Luschny, Apr 10 2022 *)

C=binomial; a(n) = sum (k=0, n, C(n,k)^2 * C(k+k,n)^2);

a(-1)=0, a(0)=1, a(n+1) = ((2*n+1)*(12*n^2+12*n+4)*a(n)-16*n^3*a(n-1))/(n+1)^3.
a(n) = Sum_{k=ceiling(n/2)..n} binomial(n,k)^2*binomial(2*k,n)^2. [Gorodetsky] - Michel Marcus, Feb 25 2021
a(n) ~ 2^(2*n - 3/4) * (1 + sqrt(2))^(2*n+1) / (Pi*n)^(3/2). - Vaclav Kotesovec, Jul 10 2021
From Peter Bala, Apr 10 2022: (Start)
The g.f. is the diagonal of the rational function 1/(1 - (x + y + z + t) + 2*(x*y*z + x*y*t + x*z*t + y*z*t) + 4*x*y*z*t) (Straub and Zudilin)
The g.f. appears to be the diagonal of the rational function 1/(1 - x - y + z - t - 2*(x*z + y*z + z*t) + 4*(x*y*t + x*z*t) + 8*x*y*z*t).
If true, then a(n) = [(x*y*z)^n] ( (x + y + z + 1)*(x + y + z - 1)*(x + y - z - 1)*(x - y - z + 1) )^n . (End)
a(n) = binomial(2*n, n)^2 * hypergeom([1/2-n/2, 1/2-n/2, -n/2, -n/2], [1, 1/2-n, 1/2-n], 1). - Peter Luschny, Apr 10 2022
G.f.: hypergeom([1/8, 3/8],[1], 256*x^2 / (1 - 4*x)^4)^2 / (1 - 4*x). - Mark van Hoeij, Nov 12 2022
a(n) = [(w*x*y*z)^n] ((w+z)*(x+z)*(y+z)*(w+x+y+z))^n = Sum_{0 <= j <= i <= n} binomial(n,i)^2*binomial(i,j)^2*binomial(n+j,i). - Jeremy Tan, Mar 28 2024

A290576 Apéry-like numbers Sum_{k=0..n} Sum_{l=0..n} (C(n,k)^2C(n,l)C(k,l)*C(k+l,n)).

Original entry on oeis.org

1, 3, 27, 309, 4059, 57753, 866349, 13492251, 216077787, 3536145057, 58875891777, 994150929951, 16984143140589, 293036113226223, 5098773125244483, 89368239352074309, 1576424378494272987, 27964450505226314673, 498550055166916502121

Offset: 0

Hugo Pfoertner, Aug 06 2017

Sequence zeta (formula 4.12) in Almkvist, Straten, Zudilin article.

Robert Israel, Table of n, a(n) for n = 0..779
G. Almkvist, D. van Straten, and W. Zudilin, Generalizations of Clausen’s formula and algebraic transformations of Calabi-Yau differential equations, Proc. Edinburgh Math. Soc.54 (2) (2011), 273-295.
Ofir Gorodetsky, New representations for all sporadic Apéry-like sequences, with applications to congruences, arXiv:2102.11839 [math.NT], 2021. See zeta p. 3.
Timothy Huber, Daniel Schultz, and Dongxi Ye, Ramanujan-Sato series for 1/pi, Acta Arith. (2023) Vol. 207, 121-160. See p. 11.
Amita Malik and Armin Straub, Divisibility properties of sporadic Apéry-like numbers, Research in Number Theory, 2016, 2:5.

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.)

f:= gfun:-rectoproc({a(0)=1, a(1)=3, a(n+1) = ( (2*n+1)*(9*n^2+9*n+3)*a(n) + 27*n^3*a(n-1) ) / (n+1)^3}, a(n), remember):
map(f, [$0..30]); # Robert Israel, Aug 07 2017

Table[Sum[Sum[(Binomial[n, k]^2*Binomial[n, j] Binomial[k, j] Binomial[k + j, n]), {j, 0, n} ], {k, 0, n}], {n, 0, 18}] (* Michael De Vlieger, Aug 07 2017 *)

C=binomial;
a(n) = sum(k=0,n, sum(l=0,n, C(n,k)^2 * C(n,l) * C(k,l) * C(k+l,n) ));

a(0) = 1, a(1) = 3,
a(n+1) = ( (2*n+1)*(9*n^2+9*n+3)*a(n) + 27*n^3*a(n-1) ) / (n+1)^3.
a(n) ~ 3^(3*n/2 + 1) * (1+sqrt(3))^(2*n+1) / (2^(n + 5/2) * (Pi*n)^(3/2)). - Vaclav Kotesovec, Jul 10 2021
G.f.: hypergeom([1/12,5/12],[1],(12*x/(1-6*x-27*x^2))^3)^2/(1-6*x-27*x^2)^(1/2). - Mark van Hoeij, Nov 11 2022

A143003 a(0) = 0, a(1) = 1, a(n+1) = (2n+1)(n^2+n+5)a(n) - n^6a(n-1).

Original entry on oeis.org

0, 1, 21, 1091, 114520, 21298264, 6410456640, 2923097201856, 1920450126458880, 1747596822651334656, 2133806329230225408000, 3405545462439659704320000, 6950705677729940374290432000, 17807686090745585163974737920000

Offset: 0

Peter Bala, Jul 19 2008

This is the case m = 1 of the general recurrence a(0) = 0, a(1) = 1, a(n+1) = (2*n+1)*(n^2+n+2*m^2+2*m+1)*a(n) - n^6*a(n-1) (we suppress the dependence of a(n) on m), which arises when accelerating the convergence of the series Sum_{k>=1} 1/k^3 for Apery's constant zeta(3). For other cases see A066989 (m=0), A143004 (m=2), A143005 (m=3) and A143006 (m=4).
The solution to the general recurrence may be expressed as a sum: a(n) = n!^3*p_m(n)*Sum_{k = 1..n} 1/(k^3*p_m(k-1)*p_m(k)), where p_m(x) = Sum_{k = 0..n} C(2*k,k)^2*C(n+k,2*k)*C(x+k,2*k) is a polynomial in x of degree 2*m.
The first few are p_0(x) = 1, p_1(x) = 2*x^2 + 2*x + 1, p_2(x) = (3*x^4 + 6*x^3 + 9*x^2 + 6*x + 2)/2 and p_3(x) = (10*x^6 + 30*x^5 + 85*x^4 + 120*x^3 + 121*x^2 + 66*x + 18)/18. For fixed n, the sequence [p_n(k)]k>=0 is the crystal ball sequence for the product lattice A_n x A_n. See A143007 for the table of values [p_n(k)] n,k >= 0. Observe that [p_n(n)] n >= 0 is the sequence of Apery numbers A005259.
The reciprocity law p_m(n) = p_n(m) holds for nonnegative integers m and n. In particular we have p_m(1) = 2*m^2 + 2*m + 1 and p_m(2) = (3*m^4 + 6*m^3 + 9*m^2 + 6*m + 2)/2.
The polynomial p_m(x) is the unique polynomial solution of the difference equation (x+1)^3*f(x+1) + x^3*f(x-1) = (2*x+1)*(x^2+x+2*m^2+2*m+1)*f(x), normalized so that f(0) = 1. The reciprocity law now yields the Apery-like recursion m^3*p_m(x) + (m-1)^3*p_(m-2)(x) = (2*m-1)*(m^2-m+1+2*x^2+2*x)*p_(m-1)(x).
The polynomial functions p_m(x) have their zeros on the vertical line Re x = -1/2 in the complex plane; that is, the polynomials p_m(x-1), m = 1,2,3,..., satisfy a Riemann hypothesis (adapt the proof of the lemma on p. 4 of [BUMP et al.]).
The general recurrence in the first paragraph above has a second solution b(n) = n!^3*p_m(n) with initial conditions b(0) = 1, b(1) = 2*m^2+2*m+1. Hence the behavior of a(n) for large n is given by lim_{n -> infinity} a(n)/b(n) = Sum_{k>=1} 1/(k^3*p_m(k-1)*p_m(k)) = 1/((2*m^2+2*m+1) - 1^6/(3*(2*m^2+2*m+3) - 2^6/(5*(2*m^2+2*m+7) - 3^6/(7*(2*m^2+2*m+13) - ...)))) = Sum_{k>=1} 1/(m+k)^3. The final equality follows from a result of Ramanujan; see [Berndt, Chapter 12, Entry 32(iii)].
For the corresponding results for the constant zeta(2) see A142995. For corresponding results for the constant log(2) see A142979 and A142992.

Bruce C. Berndt, Ramanujan's Notebooks Part II, Springer-Verlag.

Seiichi Manyama, Table of n, a(n) for n = 0..181
D. Bump, K. Choi, P. Kurlberg and J. Vaaler, A local Riemann hypothesis, I, Math. Zeit. 233, (2000), 1-19.

Cf. A066989, A008288, A108625, A142979, A142992, A142995, A143004, A143005, A143006, A143007.

p := n -> 2*n^2+2*n+1: a := n -> n!^3*p(n)*sum (1/(k^3*p(k-1)*p(k)), k = 1..n): seq(a(n), n = 0..14)

RecurrenceTable[{a[0]==0,a[1]==1,a[n+1]==(2n+1)(n^2+n+5)a[n]- n^6 a[n-1]}, a[n],{n,15}] (* Harvey P. Dale, Jun 20 2011 *)

a(n) = n!^3*p(n)*Sum_{k = 1..n} 1/(k^3*p(k-1)*p(k)), where p(n) = 2*n^2 + 2*n + 1 = A001844(n).
Recurrence: a(0) = 0, a(1) = 1, a(n+1) = (2*n+1)*(n^2+n+5)*a(n) - n^6*a(n-1).
The sequence b(n):= n!^3*p(n) satisfies the same recurrence with the initial conditions b(0) = 1, b(1) = 5. Hence we obtain the finite continued fraction expansion a(n)/b(n) = 1/(5 - 1^6/(21 - 2^6/(55 - 3^6/(119 - ... - (n-1)^6/((2*n-1)*(n^2-n+5)))))), for n >= 2. The behavior of a(n) for large n is given by lim_{n -> infinity} a(n)/b(n) = Sum_{k>=1} 1/(k^3*(4*k^4 + 1)) = 1/(5 - 1^6/(21 - 2^6/(55 - 3^6/(119 - ... - n^6/((2*n+1)*(n^2+n+5) - ...))))) = zeta(3) - 1, where the final equality follows from a result of Ramanujan; see [Berndt, Chapter 12, Entry 32(iii) at x = 1].

A183204 Central terms of triangle A181544.

Original entry on oeis.org

1, 4, 48, 760, 13840, 273504, 5703096, 123519792, 2751843600, 62659854400, 1451780950048, 34116354472512, 811208174862904, 19481055861877120, 471822589361293680, 11511531876280913760, 282665135367572129040

Offset: 0

Paul D. Hanna, Dec 30 2010

The g.f. for row n of triangle A181544 is (1-x)^(3n+1)*Sum_{k>=0}C(n+k-1,k)^3*x^k.
This sequence is s_7 in Cooper's paper. - Jason Kimberley, Nov 06 2012
Diagonal of the rational function R(x,y,z,w) = 1/(1 - (w*x*y + w*x*z + w*y*z + x*y + x*z + y + z)). - Gheorghe Coserea, Jul 14 2016
This is one of the Apery-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017
Every prime eventually divides some term of this sequence. - Amita Malik, Aug 20 2017

			Triangle A181544 begins:
(1);
1, (4), 1;
1, 20, (48), 20, 1;
1, 54, 405, (760), 405, 54, 1;
1, 112, 1828, 8464, (13840), 8464, 1828, 112, 1; ...

Seiichi Manyama, Table of n, a(n) for n = 0..702 (terms 0..499 from Jason Kimberley)
S. Cooper, Sporadic sequences, modular forms and new series for 1/pi, Ramanujan J., December 2012, Volume 29, Issue 1, pp 163-183.
Shaun Cooper, Jesús Guillera, Armin Straub, and Wadim Zudilin, Crouching AGM, Hidden Modularity, arXiv:1604.01106 [math.NT], 5-April-2016.
Ofir Gorodetsky, New representations for all sporadic Apéry-like sequences, with applications to congruences, arXiv:2102.11839 [math.NT], 2021. See s7 p. 3.
Timothy Huber, Daniel Schultz, and Dongxi Ye, Ramanujan-Sato series for 1/pi, Acta Arith. (2023) Vol. 207, 121-160. See p. 11.
Amita Malik and Armin Straub, Divisibility properties of sporadic Apéry-like numbers, Research in Number Theory, 2016, 2:5.
Robert Osburn, Armin Straub, and Wadim Zudilin, A modular supercongruence for 6F5: an Apéry-like story, arXiv:1701.04098 [math.NT], 2017.
Wadim Zudilin, A generating function of the squares of Legendre polynomials, preprint arXiv:1210.2493 [math.CA], 2012.

Cf. A108625, A181544, A183205, A245086.
Related to diagonal of rational functions: A268545-A268555.
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.)

P:=PolynomialRing(Integers()); C:=Binomial;
A183204:=func; // or directly:
A183204:=func;
[A183204(n):n in[0..16]]; // Jason Kimberley, Oct 29 2012

Table[Sum[Binomial[n,j]^2 * Binomial[2*j,n] * Binomial[j+n,j],{j,0,n}],{n,0,20}] (* Vaclav Kotesovec, Apr 05 2015 *)

{a(n)=polcoeff((1-x)^(3*n+1)*sum(j=0, 2*n, binomial(n+j, j)^3*x^j), n)}

a(n) = [x^n] (1-x)^(3n+1) * Sum_{k>=0} C(n+k-1,k)^3*x^k.
a(n) = Sum_{j = 0..n} C(n,j)^2 * C(2*j,n) * C(j+n,j). [Formula of Wadim Zudilin provided by Jason Kimberley, Nov 06 2012]
1/Pi = sqrt(7) Sum_{n>=0} (-1)^n a(n) (11895n + 1286)/22^(3n+3). [Cooper, equation (41)] - Jason Kimberley, Nov 06 2012
G.f.: sqrt((1-13*x+(1-26*x-27*x^2)^(1/2))/(1-21*x+8*x^2+(1-8*x)*(1-26*x-27*x^2)^(1/2)))*hypergeom([1/12,5/12],[1],13824*x^7/(1-21*x+8*x^2+(1-8*x)*(1-26*x-27*x^2)^(1/2))^3)^2. - Mark van Hoeij, May 07 2013
a(n) ~ 3^(3*n+3/2) / (4 * (Pi*n)^(3/2)). - Vaclav Kotesovec, Apr 05 2015
G.f. A(x) satisfies 1/(1+4*x)^2 * A( x/(1+4*x)^3 ) = 1/(1+2*x)^2 * A( x^2/(1+2*x)^3 ) [see Cooper, Guillera, Straub, Zudilin]. - Joerg Arndt, Apr 08 2016
a(n) = (-1)^n*binomial(3n+1,n)* 4F3({-n,n+1,n+1,n+1};{1,1,2(n+1)}; 1). - M. Lawrence Glasser, May 15 2016
Conjecture D-finite with recurrence: n^3*a(n) - (2*n-1)*(13*n^2-13*n+4)*a(n-1) - 3*(n-1)*(3*n-4)*(3*n-2)*a(n-2) = 0. - R. J. Mathar, May 15 2016
0 = (-x^2+26*x^3+27*x^4)*y''' + (-3*x+117*x^2+162*x^3)*y'' + (-1+86*x+186*x^2)*y' + (4+24*x)*y, where y is g.f. - Gheorghe Coserea, Jul 14 2016
From Jeremy Tan, Mar 14 2024: (Start)
The conjectured D-finite recurrence can be proved by Zeilberger's algorithm.
a(n) = Sum_{k=0..n} binomial(n,k)^2 * binomial(n+k,n) * binomial(2*n-k,n) = [(w*x*y*z)^n] ((w+y)*(x+z)*(y+z)*(w+x+y+z))^n. (End)
a(n) = Sum_{0 <= j, k <= n} binomial(n, k)^2 * binomial(n, j)^2 * binomial(k+j, n) = Sum_{k = 0..n} binomial(n, k)^2 * A108625(n, k). - Peter Bala, Jul 08 2024
From Peter Bala, Sep 18 2024: (Start)
a(n) = Sum_{k = 0..n} (-1)^(n-k)*binomial(n+k, k)^3*binomial(3*n+1, n-k). Cf A245086.
a(n) = Sum_{k = 0..n} (-1)^(n+k)*binomial(n, k)*A143007(n, k) (verified using the MulZeil procedure in Doron Zeilberger's MultiZeilberger package). (End)

cp's OEIS Frontend

A063007 T(n,k) = binomial(n,k)*binomial(n+k,k), 0 <= k <= n, triangle read by rows.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

Formula

A143007 Square array, read by antidiagonals, where row n equals the crystal ball sequence for the 2*n-dimensional lattice A_n x A_n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

SageMath

Formula

Extensions

A006077 (n+1)^2*a(n+1) = (9n^2+9n+3)*a(n) - 27*n^2*a(n-1), with a(0) = 1 and a(1) = 3.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

PARI

PARI

Formula

Extensions

A093388 (n+1)^2*a(n+1) = (17n^2+17n+6)*a(n) - 72*n^2*a(n-1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

A125143 Almkvist-Zudilin numbers: Sum_{k=0..n} (-1)^(n-k) * ((3^(n-3*k) * (3*k)!) / (k!)^3) * binomial(n,3*k) * binomial(n+k,k).

Views

Author

Keywords

Comments

References

Links

Crossrefs

A006077 (n+1)^2a(n+1) = (9n^2+9n+3)a(n) - 27n^2a(n-1), with a(0) = 1 and a(1) = 3.

A093388 (n+1)^2a(n+1) = (17n^2+17n+6)a(n) - 72n^2a(n-1).

A125143 Almkvist-Zudilin numbers: Sum_{k=0..n} (-1)^(n-k) * ((3^(n-3k) (3k)!) / (k!)^3) binomial(n,3k) binomial(n+k,k).

A290576 Apéry-like numbers Sum_{k=0..n} Sum_{l=0..n} (C(n,k)^2C(n,l)C(k,l)*C(k+l,n)).

A143003 a(0) = 0, a(1) = 1, a(n+1) = (2n+1)(n^2+n+5)a(n) - n^6a(n-1).