OFFSET
0,2
COMMENTS
In general, for m>0, Product_{k=1..n} (m*k^2+1) is asymptotic to 2*m^(n+1/2) * n^(2*n+1) * sinh(Pi/sqrt(m)) / exp(2*n).
FORMULA
a(n) = (-1)^(n+1) * A101928(n+2).
a(n) ~ 2^(2*n+2) * n^(2*n+1) * sinh(Pi/2) / exp(2*n).
a(n) = 2^(2*n+1) * |Gamma(1 + i/2 + n)|^2 * sinh(Pi/2)/Pi. - Vladimir Reshetnikov, Oct 10 2016
MATHEMATICA
Table[Product[4*k^2+1, {k, 1, n}], {n, 0, 15}]
Round@Table[2^(2 n + 1) Abs[Gamma[1 + I/2 + n]]^2 Sinh[Pi/2]/Pi, {n, 0, 15}] (* Vladimir Reshetnikov, Oct 10 2016 *)
PROG
(PARI) a(n) = prod(k=1, n, (4*k^2+1)); \\ Michel Marcus, Oct 11 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Oct 10 2016
STATUS
approved