login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A333593 a(n) = Sum_{k = 0..n} (-1)^(n + k)*binomial(n + k - 1, k)^2. 2
1, 0, 6, 72, 910, 12000, 163086, 2266544, 32043726, 459167040, 6651400756, 97214919648, 1431514320886, 21213380196736, 316072831033350, 4731683468079072, 71128104013487310, 1073134004384407680, 16243463355081280080, 246585461357885877920 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

COMMENTS

It is known that Sum_{k = 0..2*n} (-1)^(n+k)*C(2*n,k)^2 = C(2*n,n). Here, we consider by analogy Sum_{k = 0..n} (-1)^(n+k)*C(-n,k)^2, where C(-n,k) = (-1)^k * C(n+k-1,k) for integer n and nonnegative integer k.

The sequence b(n) = C(2*n,n) of central binomial coefficients satisfies the supercongruences b(n*p^k) = b(n*p^(k-1)) ( mod p^(3*k) ) for all prime p >= 5 and any positive integers n and k -  see Mestrovic. We conjecture that the present sequence also satisfies these supercongruences. Some examples of the supercongruences are given below. For a proof of the particular case of these supercongruences, a(p) == 0 ( mod p^3 ) for prime p >= 5, see the Bala link. The sequence a(p)/(2*p^3) for prime p >= 5 begins [48, 304, 36519504, 4827806144, 109213719151680, 17975321574419440, ...]. Cf. A034602.

More generally, calculation suggests that for positive integer A and integer B, the sequence a(A,B;n) := Sum_{k = 0..A*n} (-1)^(n+k)* C(-B*n,k)^2 = Sum_{k = 0..A*n} (-1)^(n+k)*C(B*n+k-1,k)^2 may also satisfy the above supercongruences for all prime p >= 5.

LINKS

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

P. Bala, Notes on A333593

R. Mestrovic, Wolstenholme's theorem: Its Generalizations and Extensions in the last hundred and fifty years (1862-2011), arXiv:1111.3057 [math.NT], 2011.

FORMULA

a(n) ~ 2^(4*n) / (5*Pi*n). - Vaclav Kotesovec, Mar 28 2020

EXAMPLE

Examples of supercongruences:

a(11)  = 97214919648 = (2^5)*3*(7^2)*(11^3)*15527 == 0 ( mod 11^3 ).

a(2*7) - a(2) = 316072831033350 - 6 = (2^13)*3*(7^3)*11*691*4933 == 0 ( mod 7^3 ).

a(5^2) - a(5) = 3164395891098711251676512000 - 12000 = (2^5)*(5^6)*29* 124891891*1747384859327 == 0 ( mod 5^6 ).

MAPLE

seq( add( (-1)^(n+k)*binomial(n+k-1, k)^2, k = 0..n ), n = 0..25);

MATHEMATICA

Table[Binomial[2*n-1, n]^2 * HypergeometricPFQ[{1, -n, -n}, {1 - 2 n, 1 - 2 n}, -1], {n, 0, 20}] (* Vaclav Kotesovec, Mar 28 2020 *)

PROG

(PARI) a(n) = sum(k=0, n, (-1)^(n+k)*binomial(n+k-1, k)^2); \\ Michel Marcus, Mar 29 2020

CROSSREFS

Cf. A000984, A034602, A333592.

Sequence in context: A007031 A067419 A113331 * A214159 A303342 A332705

Adjacent sequences:  A333590 A333591 A333592 * A333594 A333595 A333596

KEYWORD

nonn,easy

AUTHOR

Peter Bala, Mar 27 2020

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 18 13:21 EST 2021. Contains 340254 sequences. (Running on oeis4.)