OFFSET
1,1
COMMENTS
a(2(p-1)) is divisible by p^2 for p=7,13,19,31,37,43,61,67.. A002476 Primes of form 6n + 1.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..200
FORMULA
a(n) = Sum[Sum[(2k)!/(k!)^2,{k,1,m}],{m,1,n}].
a(n) = 2 * Sum[ A079309[k], {k,1,n} ] = Sum[ A066796[k], {k,1,n} ]. - Alexander Adamchuk, Sep 01 2006
G.f.: x*(1/Sqrt[1-4*x]-1)/(x(x-1)^2) [From Harvey P. Dale, May 24 2011]
Recurrence: n*a(n) = 2*(3*n-1)*a(n-1) - (9*n-4)*a(n-2) + 2*(2*n-1)*a(n-3). - Vaclav Kotesovec, Oct 19 2012
a(n) ~ 2^(2*n+4)/(9*sqrt(Pi*n)). - Vaclav Kotesovec, Oct 19 2012
MATHEMATICA
Table[Sum[Sum[(2k)!/(k!)^2, {k, 1, m}], {m, 1, n}], {n, 1, 50}]
CoefficientList[Series[(1/Sqrt[1-4 x]-1)/((x-1)^2 x), {x, 0, 50}], x] (* Harvey P. Dale, May 24 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Jul 04 2006
STATUS
approved