OFFSET
0,2
COMMENTS
For any n > 0, n^3 divides Sum_{k=0..n-1} (2k+1)*a(k).
LINKS
Zhi-Wei Sun, Table of n, a(n) for n = 0..100
Zhi-Wei Sun, Two new kinds of numbers and their arithmetic properties, arXiv:1408.5381 [math.NT], 2014-2018.
FORMULA
Recurrence: n^3*(2*n-3)*(8*n^4 - 48*n^3 + 96*n^2 - 72*n + 13)*a(n) = (2*n-1)*(280*n^7 - 2096*n^6 + 5728*n^5 - 6536*n^4 + 1383*n^3 + 3160*n^2 - 2432*n + 552)*a(n-1) - (2*n-3)*(280*n^7 - 1824*n^6 + 4096*n^5 - 3384*n^4 - 345*n^3 + 2046*n^2 - 1100*n + 192)*a(n-2) + (n-2)^3*(2*n-1)*(8*n^4 - 16*n^3 + 8*n - 3)*a(n-3). - Vaclav Kotesovec, Aug 27 2014
a(n) ~ sqrt(24+17*sqrt(2)) * (17+12*sqrt(2))^n / (4*Pi^(3/2)*sqrt(n)). - Vaclav Kotesovec, Aug 27 2014
EXAMPLE
a(1) = 13 since Sum_{k=0..1} (2k+1)*C(1,k)^2*C(1+k,k)^2 = 1 + 3*2^2 = 13.
MAPLE
A246462:=n->add((2*k+1)*binomial(n, k)^2*binomial(n+k, k)^2, k=0..n): seq(A246462(n), n=0..20); # Wesley Ivan Hurt, Aug 27 2014
MATHEMATICA
a[n_]:=Sum[(2k+1)*Binomial[n, k]^2*Binomial[n+k, k]^2, {k, 0, n}]
Table[a[n], {n, 0, 15}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Zhi-Wei Sun, Aug 26 2014
STATUS
approved