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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A333561 a(n) = Sum_{j = 0..2*n} binomial(n+j-1,j)*2^j. 2
1, 7, 129, 2815, 65537, 1579007, 38862849, 970522623, 24494735361, 623210135551, 15956734640129, 410649406472191, 10612705274626049, 275241225206890495, 7159857331658817537, 186731505521384226815, 4880983719142471237633, 127836403093194475044863 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

Column 2 of the square array A333560. Compare with A119259(n) = Sum_{j = 0..n} binomial(n+j-1,j)*2^j.

We conjecture that this sequence satisfies the supercongruences a(n*p^k) == a(n*p^(k-1)) ( mod p^(3*k) ) for prime p >= 5 and positive integers n and k. Some examples are given below.

LINKS

Table of n, a(n) for n=0..17.

FORMULA

Conjectural o.g.f.: 1/(1 + x) + 8*x*f'(4*x)/(2*f(4*x) - 1), where f(x) = 1 + x + 3*x^2 + 12*x^3 + 55*x^4 + ... is the o.g.f. of A001764.

exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + 7*x + 89*x^2 + 1447*x^3 + ... appears to be the o.g.f. of A062747.

Conjectural recurrence: n*(n - 1)*(2*n - 1)*(3098*n - 6455)*a(n) = (n - 1)*(172988*n^3 - 585840*n^2 + 550321*n - 169824)*a(n-1) - 12*(11825*n^4 - 168518*n^3 + 627675*n^2 - 853766*n + 350744)*a(n-2) - 36*(n - 3)*(3*n - 7)*(3*n - 8)*(991*n - 724)*a(n-3) with a(1) = 7, a(2) = 129, a(3) = 2815.

From Vaclav Kotesovec, Mar 28 2020: (Start)

a(n) ~ 3^(3*n + 1/2) / (4*sqrt(Pi*n)).

Recurrence: n*(2*n - 1)*(7*n^2 - 20*n + 14)*a(n) = (364*n^4 - 1411*n^3 + 1818*n^2 - 868*n + 120)*a(n-1) + 6*(3*n - 5)*(3*n - 4)*(7*n^2 - 6*n + 1)*a(n-2). (End)

EXAMPLE

Examples of supercongruences:

a(11) - a(1) = 410649406472191 - 7 = (2^3)*3*(11^3)*12855290711 == 0 ( mod 11^3 ).

a(3*7) - a(3) = 61103847305642669128888090623 - 2815 = (2^8)*(7^5)* 87326419*162627033103121 == 0 ( mod 7^3 ).

a(5^2) - a(5) = 29754989698128108780761000609579007 - 1579007 = (2^11)*(5^6)*179*751*10267*673710468794491483 == 0 ( mod 5^6 ).

MAPLE

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

MATHEMATICA

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

PROG

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

CROSSREFS

Cf. A001764, A062747, A119259, A333560, A333562.

Sequence in context: A309130 A263686 A002614 * A187237 A095885 A332050

Adjacent sequences:  A333558 A333559 A333560 * A333562 A333563 A333564

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 October 1 13:04 EDT 2020. Contains 337443 sequences. (Running on oeis4.)