A172013 a(n) = 6*A142459(2*n, n)/(n+1).

3, 118, 20343, 8530698, 6711481694, 8575821262764, 16243345162977759, 42826533033277249154, 150138953276380791799098, 675925071086215282939520628, 3802445616812067139270851537718, 26147695687370407271086390933321188, 215852465255521412471161891166554453788

Offset

1,1

Links

G. C. Greubel, Table of n, a(n) for n = 1..180

Formula

a(n) = 6*A142459(2*n, n)/(n+1).

Mathematica

T[n_, k_, m_]:= T[n, k, m]= If[k==1 || k==n, 1, (m*n-m*k+1)*T[n-1, k-1, m] + (m*k-m+1)*T[n-1, k, m]];
A142459[n_, k_]:= A142459[n, k]= T[n, k, 4];
A172013[n_]:= A172013[n]= 6*A142459[2*n, n]/(n+1);
Table[A172013[n], {n, 30}] (* modified by G. C. Greubel, Mar 18 2022 *)

Prog

(Sage)
@CachedFunction
def T(n, k, m):
    if (k==1 or k==n): return 1
    else: return (m*(n-k)+1)*T(n-1, k-1, m) + (m*k-m+1)*T(n-1, k, m)
def A142459(n, k): return T(n, k, 4)
def A172013(n): return 6*A142459(2*n, n)/(n+1)
[A172013(n) for n in (1..30)] # G. C. Greubel, Mar 18 2022

Crossrefs

Cf. A142458, A142459, A172010.

Sequence in context: A283883 A173053 A180393 * A143781 A114077 A176996

Adjacent sequences: A172010 A172011 A172012 * A172014 A172015 A172016

Keyword

nonn

Author

Roger L. Bagula, Nov 19 2010

Extensions

Offset and formula corrected by G. C. Greubel, Mar 18 2022

Status

approved