%I #6 May 10 2021 07:37:49
%S 0,1,19,713,45963,4571521,651249603,125978555961,31797923989563,
%T 10154867346496881,4003950222788879475,1910709271283079616425,
%U 1085491754899149563498475,724022706189621081117571425,560305448143863386421257597475,497969168730434344111574554745625
%N a(n+1) = (8*n^2+8*n+3)*a(n) - 16*n^4*a(n-1), with a(0)=0, a(1)=1.
%H Karl Dilcher and Lin Jiu, <a href="https://arxiv.org/abs/2105.01880">Hankel Determinants of shifted sequences of Bernoulli and Euler numbers</a>, arXiv:2105.01880 [math.NT], 2021. See Proposition 3.5. p. 7.
%o (PARI) a(n) = if (n<=1, n, n--; (8*n^2+8*n+3)*a(n) - 16*n^4*a(n-1));
%K nonn
%O 0,3
%A _Michel Marcus_, May 10 2021
|