A309010 -id:A309010 - cp's OEIS frontend

A005260 a(n) = Sum_{k = 0..n} binomial(n,k)^4.

1, 2, 18, 164, 1810, 21252, 263844, 3395016, 44916498, 607041380, 8345319268, 116335834056, 1640651321764, 23365271704712, 335556407724360, 4854133484555664, 70666388112940818, 1034529673001901732, 15220552520052960516, 224929755893153896200, 3337324864503769353060

Offset: 0

N. J. A. Sloane

This sequence is s_10 in Cooper's paper. - Jason Kimberley, Nov 25 2012
Diagonal of the rational function R(x,y,z,w) = 1/(1 - (w*x*y + w*x*z + w*y*z + x*y*z + w*x + y*z)). - Gheorghe Coserea, Jul 13 2016
This is one of the Apéry-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017
Every prime eventually divides some term of this sequence. - Amita Malik, Aug 20 2017
Two walkers, A and B, stand on the South-West and North-East corners of an n X n grid, respectively. A walks by either North or East steps while B walks by either South or West steps. Sequence values a(n) < binomial(2*n,n)^2 count the simultaneous walks where A and B meet after exactly n steps and change places after 2*n steps. - Bradley Klee, Apr 01 2019
a(n) is the constant term in the expansion of ((1 + x) * (1 + y) * (1 + z) + (1 + 1/x) * (1 + 1/y) * (1 + 1/z))^n. - Seiichi Manyama, Oct 27 2019

			G.f. = 1 + 2*x + 18*x^2 + 164*x^3 + 1810*x^4 + 21252*x^5 + 263844*x^6 + ...

H. W. Gould, Combinatorial Identities, Morgantown, 1972, (X.14), p. 79.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Seiichi Manyama, Table of n, a(n) for n = 0..834 (terms 0..250 from Jason Kimberley)
B. Adamczewski, J. P. Bell, and E. Delaygue, Algebraic independence of G-functions and congruences "a la Lucas", arXiv preprint arXiv:1603.04187 [math.NT], 2016.
Hacene Belbachir and Yassine Otmani, A Strehl Version of Fourth Franel Sequence, arXiv:2012.02563 [math.CO], 2020.
F. Beukers, Another congruence for the Apéry numbers, J. Number Theory 25 (1987), no. 2, 201-210.
W. Y. C. Chen, Q.-H. Hou, and Y-P. Mu, A telescoping method for double summations, J. Comp. Appl. Math. 196 (2006) 553-566, eq (5.5).
S. Cooper, Sporadic sequences, modular forms and new series for 1/pi, Ramanujan J. (2012).
M. Coster, Email, Nov 1990
E. Delaygue, Arithmetic properties of Apéry-like numbers, arXiv preprint arXiv:1310.4131 [math.NT], 2013.
C. Elsner, On recurrence formulas for sums involving binomial coefficients, Fib. Q., 43,1 (2005), 31-45.
Ofir Gorodetsky, New representations for all sporadic Apéry-like sequences, with applications to congruences, arXiv:2102.11839 [math.NT], 2021. See s10 p. 3.
Darij Grinberg, Introduction to Modern Algebra (UMN Spring 2019 Math 4281 Notes), University of Minnesota (2019).
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.
M. A. Perlstadt, Some Recurrences for Sums of Powers of Binomial Coefficients, Journal of Number Theory 27 (1987), pp. 304-309.
V. Strehl, Recurrences and Legendre transform, Séminaire Lotharingien de Combinatoire, B29b (1992), 22 pp.
Zhi-Wei Sun, Open conjectures on congruences, Nanjing Univ. J. Math. Biquarterly 36(2019), no.1, 1-99. (Cf. Conjectures 49-51.)
Eric Weisstein's World of Mathematics, Binomial Sums
Mark C. Wilson, Diagonal asymptotics for products of combinatorial classes, preprint of Combinatorics, Probability and Computing, 24(1), 2015, 354-372.
Jin Yuan, Zhi-Juan Lu, and Asmus L. Schmidt, On recurrences for sums of powers of binomial coefficients, J. Numb. Theory 128 (2008) 2784-2794.

Column k=4 of A309010.
Cf. A000172, A096192, A328725.
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.)
Sum_{k = 0..n} C(n,k)^m for m = 1..12: A000079, A000984, A000172, A005260, A005261, A069865, A182421, A182422, A182446, A182447, A342294, A342295.
Row sums of A202750.

A005260 := proc(n)
        add( (binomial(n,k))^4,k=0..n) ;
end proc:
seq(A005260(n),n=0..10) ; # R. J. Mathar, Nov 19 2012

Table[Sum[Binomial[n, k]^4, {k, 0, n}], {n, 0, 20}] (* Wesley Ivan Hurt, Mar 09 2014 *)

{a(n) = sum(k=0, n, binomial(n, k)^4)};

def A005260(n):
    m, g = 1, 0
    for k in range(n+1):
        g += m
        m = m*(n-k)**4//(k+1)**4
    return g # Chai Wah Wu, Oct 04 2022

a(n) ~ 2^(1/2)*Pi^(-3/2)*n^(-3/2)*2^(4*n). - Joe Keane (jgk(AT)jgk.org), Jun 21 2002
D-finite with recurrence: n^3*a(n) = 2*(2*n - 1)*(3*n^2 - 3*n + 1)*a(n-1) + (4*n - 3)*(4*n - 4)*(4*n - 5)*a(n-2). [Yuan]
G.f.: 5*hypergeom([1/8, 3/8],[1], (4/5)*((1-16*x)^(1/2)+(1+4*x)^(1/2))*(-(1-16*x)^(1/2)+(1+4*x)^(1/2))^5/(2*(1-16*x)^(1/2)+3*(1+4*x)^(1/2))^4)^2/(2*(1-16*x)^(1/2)+3*(1+4*x)^(1/2)). - Mark van Hoeij, Oct 29 2011
1/Pi = sqrt(15)/18 * Sum_{n >= 0} a(n)*(4*n + 1)/36^n (Cooper, equation (5)) = sqrt(15)/18 * Sum_{n >= 0} a(n)*A016813(n)/A009980(n). - Jason Kimberley, Nov 26 2012
0 = (-x^2 + 12*x^3 + 64*x^4)*y''' + (-3*x + 54*x^2 + 384*x^3)*y'' + (-1 + 40*x + 444*x^2)*y' + (2 + 60*x)*y, where y is g.f. - Gheorghe Coserea, Jul 13 2016
For r a nonnegative integer, Sum_{k = r..n} C(k,r)^4*C(n,k)^4 = C(n,r)^4*a(n-r), where we take a(n) = 0 for n < 0. - Peter Bala, Jul 27 2016
a(n) = hypergeom([-n, -n, -n, -n], [1, 1, 1], 1). - Peter Luschny, Jul 27 2016
Sum_{n>=0} a(n) * x^n / (n!)^4 = (Sum_{n>=0} x^n / (n!)^4)^2. - Ilya Gutkovskiy, Jul 17 2020
a(n) = Sum_{k=0..n} C(n,k)*C(n+k,k)*C(2k,k)*C(2n-2k,n-k)*(-1)^(n-k). This can be proved via the Zeilberger algorithm. - Zhi-Wei Sun, Aug 23 2020
a(n) = (-1)^n*binomial(2*n, n)*hypergeom([1/2, -n, -n, n + 1], [1, 1, 1/2 - n], 1). - Peter Luschny, Aug 24 2020
a(n) = Sum_{k=0..n} binomial(n,k)^2*binomial(2*k,n)*binomial(2*n-k,n) [Theorem 1 in Belbachir and Otmani]. - Michel Marcus, Dec 06 2020
a(n) = [x^n] (1 - x)^(2*n) P(n,(1 + x)/(1 - x))^2, where P(n,x) denotes the n-th Legendre polynomial. See Gould, p. 66. This formula is equivalent to the binomial sum identity of Zhi-Wei Sun given above. - Peter Bala, Mar 24 2022
From Peter Bala, Oct 31 2024: (Start)
For n >= 1, a(n) = 2 * Sum_{k = 0..n-1} binomial(n, k)^3 * binomial(n-1, k).
For n >= 1, a(n) = 2 * hypergeom([-n, -n, -n, -n + 1], [1, 1, 1], 1). (End)
G.f.: Sum_{k>=0} Sum_{l=0..p*k} Sum_{m=0..l} (-1)^m*binomial(p*k+1,m)*binomial(l+k-m,k)^p*x^(l+k)/(1-x)^(p*k+1), where p = 4. - Miles Wilson, Apr 12 2025

Edited by Michael Somos, Aug 09 2002

Minor edits by Vaclav Kotesovec, Aug 28 2014

A005261 a(n) = Sum_{k = 0..n} C(n,k)^5.

Original entry on oeis.org

1, 2, 34, 488, 9826, 206252, 4734304, 113245568, 2816649826, 72001228052, 1883210876284, 50168588906768, 1357245464138656, 37198352117916992, 1030920212982957184, 28847760730478655488, 814066783370083977826

Offset: 0

N. J. A. Sloane

a(n) is the constant term in the expansion of ((1 + w) * (1 + x) * (1 + y) * (1 + z) + (1 + 1/w) * (1 + 1/x) * (1 + 1/y) * (1 + 1/z))^n. - Seiichi Manyama, Oct 27 2019

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 0..200
C. Elsner, On recurrence formulas for sums involving binomial coefficients, Fib. Q., 43,1 (2005), 31-45.
M. A. Perlstadt, Some Recurrences for Sums of Powers of Binomial Coefficients, Journal of Number Theory 27 (1987), pp. 304-309.
Eric Weisstein's World of Mathematics, Binomial sums.
Mark C. Wilson, Diagonal asymptotics for products of combinatorial classes.
Jin Yuan, Zhi-Juan Lu, Asmus L. Schmidt, On recurrences for sums of powers of binomial coefficients, J. Numb. Theory 128 (2008) 2784-2794

Column k=5 of A309010.

Sum_{k = 0..n} C(n,k)^m for m = 1..12: A000079, A000984, A000172, A005260, A005261, A069865, A182421, A182422, A182446, A182447, A342294, A342295.

a := n -> hypergeom([seq(-n, i=1..5)],[seq(1, i=1..4)], -1):
seq(simplify(a(n)),n=0..16); # Peter Luschny, Jul 27 2016

RecurrenceTable[{32*(55n^2+33n+6)*(n-1)^4*a[n-2]-(19415n^6-27181n^5+7453n^4+3289n^3-956n^2-276n+96)*a[n-1]-(1155n^6+693n^5-732n^4-715n^3+45n^2+210n+56)*a[n]+(55n^2-77n+28)*(n+1)^4*a[n+1]==0,a[0]==1,a[1]==2,a[2]==34},a,{n,0,25}]
(* or directly *)
Table[Sum[Binomial[n,k]^5,{k,0,n}],{n,0,25}] (* Vaclav Kotesovec, Apr 27 2012 *)

a(n) = sum(k=0, n, binomial(n, k)^5); \\ Michel Marcus, Mar 09 2016

def A005261(n):
    m, g = 1, 0
    for k in range(n+1):
        g += m
        m = m*(n-k)**5//(k+1)**5
    return g # Chai Wah Wu, Oct 04 2022

a(n) ~ 4*5^(-1/2)*Pi^-2*n^-2*2^(5*n). - Joe Keane (jgk(AT)jgk.org), Jun 21 2002
Recurrence (M. A. Perlstadt, 1987): 32*(55*n^2 + 33*n + 6)*(n - 1)^4*a(n-2) - (19415*n^6 - 27181*n^5 + 7453*n^4 + 3289*n^3 - 956*n^2 - 276*n + 96)*a(n-1) - (1155*n^6 + 693*n^5 - 732*n^4 - 715*n^3 + 45*n^2 + 210*n + 56)*a(n) + (55*n^2 - 77*n + 28)*(n + 1)^4*a(n+1) = 0. [Vaclav Kotesovec, Apr 27 2012]
For r a nonnegative integer, Sum_{k = r..n} C(k,r)^5*C(n,k)^5 = C(n,r)^5*a(n-r), where we take a(n) = 0 for n < 0. - Peter Bala, Jul 27 2016
Sum_{n>=0} a(n) * x^n / (n!)^5 = (Sum_{n>=0} x^n / (n!)^5)^2. - Ilya Gutkovskiy, Jul 17 2020
From Peter Bala, Nov 01 2024: (Start)
For n >= 1, a(n) = 2 * Sum_{k = 0..n-1} binomial(n, k)^4 * binomial(n-1, k).
For n >= 1, a(n) = 2 * hypergeom([-n, -n, -n, -n, -n + 1], [1, 1, 1, 1], -1). (End)

More terms from Matthew Conroy, Mar 16 2006

A342294 a(n) = Sum_{k = 0..n} binomial(n,k)^11.

Original entry on oeis.org

1, 2, 2050, 354296, 371185666, 200097656252, 222100237312864, 193798873701831680, 231719476114879600642, 257097895846251291074612, 330463219813679264204224300, 419460465362069257397304825200, 573863850341313751827291703127200

Offset: 0

N. J. A. Sloane, Mar 27 2021

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..306
M. A. Perlstadt, Some Recurrences for Sums of Powers of Binomial Coefficients, Journal of Number Theory 27 (1987), pp. 304-309.

Column 11 of A309010.

Sum_{k = 0..n} C(n,k)^m for m = 1..12: A000079, A000984, A000172, A005260, A005261, A069865, A182421, A182422, A182446, A182447, A342294, A342295.

Table[Sum[Binomial[n,k]^11,{k,0,n}],{n,0,15}] (* Harvey P. Dale, May 08 2025 *)

a(n) = sum(k=0, n, binomial(n, k)^11); \\ Michel Marcus, Mar 27 2021

a(n) ~ 2^(p*n)/sqrt(p) * (2/(Pi*n))^((p-1)/2) * (1 - (p-1)^2/(4*p*n)), set p=11. - Vaclav Kotesovec, Aug 04 2022

A342295 a(n) = Sum_{k = 0..n} binomial(n,k)^12.

Original entry on oeis.org

1, 2, 4098, 1062884, 2210336770, 2000488281252, 4355497029345924, 6773152698818628936, 15744083665278445490178, 32270900877696351763796420, 80314784333143089874093429348, 192454957455454582636391397662856, 509571049488109525160616367158261124

Offset: 0

N. J. A. Sloane, Mar 27 2021

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..281
M. A. Perlstadt, Some Recurrences for Sums of Powers of Binomial Coefficients, Journal of Number Theory 27 (1987), pp. 304-309.

Column 12 of A309010.

Sum_{k = 0..n} C(n,k)^m for m = 1..12: A000079, A000984, A000172, A005260, A005261, A069865, A182421, A182422, A182446, A182447, A342294, A342295.

a(n) = sum(k=0, n, binomial(n, k)^12); \\ Michel Marcus, Mar 27 2021

a(n) ~ 2^(p*n)/sqrt(p) * (2/(Pi*n))^((p-1)/2) * (1 - (p-1)^2/(4*p*n)), set p=12. - Vaclav Kotesovec, Aug 04 2022

A328747 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) is Sum_{i=0..n} (-1)^(n-i)binomial(n,i)Sum_{j=0..i} binomial(i,j)^k.

Original entry on oeis.org

1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 3, 1, 0, 1, 1, 7, 7, 1, 0, 1, 1, 15, 31, 19, 1, 0, 1, 1, 31, 115, 175, 51, 1, 0, 1, 1, 63, 391, 1255, 991, 141, 1, 0, 1, 1, 127, 1267, 8071, 13671, 5881, 393, 1, 0, 1, 1, 255, 3991, 49399, 161671, 160461, 35617, 1107, 1, 0

Offset: 0

Seiichi Manyama, Oct 27 2019

T(n,k) is the constant term in the expansion of (-1 + Product_{j=1..k-1} (1 + x_j) + Product_{j=1..k-1} (1 + 1/x_j))^n for k > 0.

For fixed k > 0, T(n,k) ~ (2^k - 1)^(n + (k-1)/2) / (2^((k-1)^2/2) * sqrt(k) * (Pi*n)^((k-1)/2)). - Vaclav Kotesovec, Oct 28 2019

			Square array begins:
   1, 1,  1,   1,     1,      1, ...
   1, 1,  1,   1,     1,      1, ...
   0, 1,  3,   7,    15,     31, ...
   0, 1,  7,  31,   115,    391, ...
   0, 1, 19, 175,  1255,   8071, ...
   0, 1, 51, 991, 13671, 161671, ...

Seiichi Manyama, Antidiagonals n = 0..100, flattened

Columns k=0..5 give A019590(n+1), A000012, A002426, A172634, A328725, A328750.
Main diagonal gives A328811.
T(n,n+1) gives A328813.
Cf. A309010, A328748, A328807.

T[n_, k_] := Sum[(-1)^(n-i) * Binomial[n, i] * Sum[Binomial[i, j]^k, {j, 0, i}], {i, 0, n}]; Table[T[k, n - k], {n, 0, 10}, {k, 0, n}] // Flatten (* Amiram Eldar, May 06 2021 *)

A328748 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) is Sum_{i=0..n} (-2)^(n-i)binomial(n,i)Sum_{j=0..i} binomial(i,j)^k.

Original entry on oeis.org

1, 1, 0, 1, 0, -1, 1, 0, 0, 2, 1, 0, 2, 0, -3, 1, 0, 6, 0, 0, 4, 1, 0, 14, 12, 6, 0, -5, 1, 0, 30, 72, 90, 0, 0, 6, 1, 0, 62, 300, 882, 360, 20, 0, -7, 1, 0, 126, 1080, 6690, 8400, 2040, 0, 0, 8, 1, 0, 254, 3612, 44706, 124920, 95180, 10080, 70, 0, -9

Offset: 0

Seiichi Manyama, Oct 27 2019

T(n,k) is the constant term in the expansion of (-2 + Product_{j=1..k-1} (1 + x_j) + Product_{j=1..k-1} (1 + 1/x_j))^n for k > 0.

			Square array begins:
    1, 1, 1,   1,    1,      1, ...
    0, 0, 0,   0,    0,      0, ...
   -1, 0, 2,   6,   14,     30, ...
    2, 0, 0,  12,   72,    300, ...
   -3, 0, 6,  90,  882,   6690, ...
    4, 0, 0, 360, 8400, 124920, ...

Seiichi Manyama, Antidiagonals n = 0..100, flattened

Columns k=0..5 give A097141(n+1), A000007, A126869, A002898, A328735, A328751.
T(n,n+1) gives A328814.
Cf. A309010, A328747, A328807.

T[n_, k_] := Sum[(-2)^(n-i) * Binomial[n, i] * Sum[Binomial[i, j]^k, {j, 0, i}], {i, 0, n}]; Table[T[k, n - k], {n, 0, 10}, {k, 0, n}] // Flatten (* Amiram Eldar, May 06 2021 *)

A328807 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) is Sum_{i=0..n} binomial(n,i)*Sum_{j=0..i} binomial(i,j)^k.

Original entry on oeis.org

1, 1, 3, 1, 3, 8, 1, 3, 9, 20, 1, 3, 11, 27, 48, 1, 3, 15, 45, 81, 112, 1, 3, 23, 93, 195, 243, 256, 1, 3, 39, 225, 639, 873, 729, 576, 1, 3, 71, 597, 2583, 4653, 3989, 2187, 1280, 1, 3, 135, 1665, 11991, 32133, 35169, 18483, 6561, 2816

Offset: 0

Seiichi Manyama, Oct 28 2019

T(n,k) is the constant term in the expansion of (1 + Product_{j=1..k-1} (1 + x_j) + Product_{j=1..k-1} (1 + 1/x_j))^n for k > 0.

For fixed k > 0 is T(n,k) ~ (2^k + 1)^(n + (k-1)/2) / (2^((k-1)^2/2) * sqrt(k) * (Pi*n)^((k-1)/2)). - Vaclav Kotesovec, Oct 28 2019

			Square array begins:
     1,   1,   1,    1,     1,      1, ...
     3,   3,   3,    3,     3,      3, ...
     8,   9,  11,   15,    23,     39, ...
    20,  27,  45,   93,   225,    597, ...
    48,  81, 195,  639,  2583,  11991, ...
   112, 243, 873, 4653, 32133, 260613, ...

Seiichi Manyama, Antidiagonals n = 0..100, flattened

Columns k=0..5 give A001792, A000244, A026375, A002893, A328808, A328809.
Main diagonal gives A328810.
Cf. A309010, A328747, A328748.

T[n_, k_] := Sum[Binomial[n, i] * Sum[Binomial[i, j]^k, {j, 0, i}], {i, 0, n}]; Table[T[k, n - k], {n, 0, 9}, {k, 0, n}] // Flatten (* Amiram Eldar, May 06 2021 *)

cp's OEIS Frontend

A005260 a(n) = Sum_{k = 0..n} binomial(n,k)^4.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A005261 a(n) = Sum_{k = 0..n} C(n,k)^5.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A342294 a(n) = Sum_{k = 0..n} binomial(n,k)^11.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A342295 a(n) = Sum_{k = 0..n} binomial(n,k)^12.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A328747 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) is Sum_{i=0..n} (-1)^(n-i)*binomial(n,i)*Sum_{j=0..i} binomial(i,j)^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A328748 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) is Sum_{i=0..n} (-2)^(n-i)*binomial(n,i)*Sum_{j=0..i} binomial(i,j)^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A328807 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) is Sum_{i=0..n} binomial(n,i)*Sum_{j=0..i} binomial(i,j)^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

A328747 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) is Sum_{i=0..n} (-1)^(n-i)binomial(n,i)Sum_{j=0..i} binomial(i,j)^k.

A328748 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) is Sum_{i=0..n} (-2)^(n-i)binomial(n,i)Sum_{j=0..i} binomial(i,j)^k.

A350594 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) is Sum_{j=0..2n} (-1)^(n+j) binomial(2*n,j)^k.