|
%I #31 Apr 24 2018 02:10:50
%S 1,2,25,800,48841,4867200,719580625,147968000000,40399053800625,
%T 14140937699532800,6174655078400355625,3290389182409605120000,
%U 2101698235513021884765625,1585118602783467315200000000,1393789829051727854522489390625
%N a(n) = Product_{k=1..n} (k^2+(n-k+1)^2).
%H Seiichi Manyama, <a href="/A301616/b301616.txt">Table of n, a(n) for n = 0..200</a>
%F a(n) = A302661(n)^2 + A302662(n)^2.
%F a(n) ~ n^(2*n) / exp(2*n - Pi*(n + 1)/2). - _Vaclav Kotesovec_, Apr 11 2018
%p seq(mul(k^2+(n-k+1)^2,k=1..n),n=0..15); # _Muniru A Asiru_, Apr 11 2018
%t Table[Product[k^2 + (n - k + 1)^2, {k, 1, n}], {n, 0, 15}] (* _Vaclav Kotesovec_, Apr 11 2018 *)
%o (PARI) {a(n) = prod(k=1, n, k^2+(n-k+1)^2)}
%o (GAP) List([0..15],n->Product([1..n],k->k^2+(n-k+1)^2)); # _Muniru A Asiru_, Apr 11 2018
%Y Cf. A302661, A302662.
%K nonn
%O 0,2
%A _Seiichi Manyama_, Apr 11 2018
|
