A178791 The arithmetic mean of (2k+1)(-1)^k*A_k (k=0,...,n-1), where A_0, A_1,... are Apéry numbers given by A005259.
1, -7, 117, -2441, 57449, -1453635, 38609845, -1061792695, 29973352185, -863536596143, 25288254409373, -750531594051981, 22525211241191881, -682459907754004723, 20845409947239778533, -641211780685502724425
Offset: 1
Keywords
Examples
For n=3 we have a(3) = (A_0 - 3A_1 + 5A_2)/3 = (1 - 3*5 + 5*73)/3 = 117.
Links
- G. C. Greubel, Table of n, a(n) for n = 1..500
- Zhi-Wei Sun, Arithmetic properties of Apery numbers and central Delannoy numbers, arXiv:1006.2776 [math.NT], 2010-2011.
Programs
-
Maple
G := (-1/2)*(3*x-3+(x^2-34*x+1)^(1/2))*(x+1)^(-2)*hypergeom([1/3, 2/3], [1], (-1/2)*(x^2-7*x+1)*(x+1)^(-3)*(x^2-34*x+1)^(1/2)+(1/2)*(x^3+30*x^2 -24*x+1)*(x+1)^(-3))^2; ogf := 2*x*G/(x+1)+Int((x-1)*G/(x+1)^2,x); series(ogf, x=0, 25); series(-subs(x=-x,%), x=0, 25); # Mark van Hoeij, May 07 2013
-
Mathematica
Apery[n_]:= Sum[Binomial[n+k,k]^2Binomial[n,k]^2,{k,0,n}]; AA[n_]:= Sum[(2k+1)(-1)^k*Apery[k],{k,0,n-1}]/n; Table[AA[n],{n,25}]
-
PARI
A(n) = sum(k=0, n, (binomial(n, k)*binomial(n+k, k))^2); \\ A005259 a(n) = sum(k=0, n-1, (2*k+1)*(-1)^k*A(k))/n; \\ Michel Marcus, Jan 24 2019
Formula
G.f.: apart from the minus signs (just replace x by -x) the generating function is 2*x*G/(x+1) + Int((x-1)*G/(x+1)^2, x) where G is the generating function of A005259. - Mark van Hoeij, May 07 2013
a(n) ~ -(-1)^n * 2^(3/4) * (1 + sqrt(2))^(4*n) / (24 * (Pi*n)^(3/2)). - Vaclav Kotesovec, Jan 24 2019
Comments