A348315 a(n) = Sum_{k=0..n} binomial(n^2 - k,n*k).

1, 1, 4, 64, 4382, 1357136, 1597653852, 8389021518585, 164828345435877580, 14256525628649472111712, 4602970880920727147946847283, 6484132480933772335644792339409450, 34112054985056318746734374876035089268523

Offset

0,3

Links

Seiichi Manyama, Table of n, a(n) for n = 0..57

Formula

a(n) = A306680(n^2,n) = [x^(n^2)] (1-x)^(n-1)/((1-x)^n - x^(n+1)).

a(n) ~ c * 2^(1/2 - n/2 + n^2) / (sqrt(Pi)*exp(1/8)*n), where c = Sum_{m = -infinity..+infinity} 1/(2^m * exp(m*(2*m+1)) = 1.77058122254033174512511... if n is even and c = Sum_{m = -infinity..+infinity} 1/(2^(m + 1/2) * exp((m+1)*(2*m+1))) = 1.81629595919505881855931... if n is odd. - Vaclav Kotesovec, Oct 12 2021

Mathematica

a[n_] := Sum[Binomial[n^2 - k, n*k], {k, 0, n}]; Array[a, 13, 0] (* Amiram Eldar, Oct 12 2021 *)

Prog

(PARI) a(n) = sum(k=0, n, binomial(n^2-k, n*k));
(PARI) a(n) = polcoef((1-x)^(n-1)/((1-x)^n-x^(n+1)+x*O(x^n^2)), n^2);

Crossrefs

Cf. A167009, A238696, A306680, A348322.

Sequence in context: A193611 A193755 A194501 * A053923 A359231 A326868

Adjacent sequences: A348312 A348313 A348314 * A348316 A348317 A348318

Keyword

nonn

Author

Seiichi Manyama, Oct 11 2021

Status

approved