A103212 a(n) = (1/n) * Sum_{i=0..n-1} C(n,i)*C(n,i+1)*(n-1)^i*n^(n-i) for n>=1, a(0)=1.

1, 1, 6, 93, 2380, 85405, 3956106, 224939113, 15175702200, 1185580310121, 105302043709390, 10482085765658661, 1156062800841590148, 139945327558704629221, 18449221488652046992914, 2631255715262150125502865, 403689862107153669227378416, 66297391981691913179574751633

Offset

0,3

Links

Andrew Howroyd, Table of n, a(n) for n = 0..200

Formula

a(n) = A103209(n, n-1). [corrected by Vaclav Kotesovec, Sep 24 2017]

a(n) ~ 2^(2*n) * n^(n-3/2) / (sqrt(Pi) * exp(1/2)). - Vaclav Kotesovec, Sep 24 2017

Mathematica

Table[HypergeometricPFQ[{-n, n+1}, {2}, -n+1], {n, 0, 20}] (* Vaclav Kotesovec, Sep 24 2017 *)
Flatten[{1, 1, Table[Sum[Binomial[n, k]*Binomial[n, k+1]*(n-1)^k*n^(n-k), {k, 0, n-1}]/n, {n, 2, 20}]}] (* Vaclav Kotesovec, Sep 24 2017 *)

Prog

(PARI) a(n) = {if(n==0, 1, sum(i=0, n-1, binomial(n, i)*binomial(n, i+1)*(n-1)^i*n^(n-i))/n)} \\ Andrew Howroyd, Apr 14 2021

Crossrefs

Cf. A103209, A292798.

Sequence in context: A331623 A158121 A328427 * A359928 A033935 A218682

Adjacent sequences: A103209 A103210 A103211 * A103213 A103214 A103215

Keyword

nonn

Author

Ralf Stephan, Jan 27 2005

Extensions

Prepended a(0)=1 from Vaclav Kotesovec, Sep 24 2017

Terms a(15) and beyond from Andrew Howroyd, Apr 14 2021

Status

approved