A127606 a(n) = 2^(2*n*n) * Product_{1<=i,j<=n} (cos(i*Pi/(2*n+1))^2 + sin(j*Pi/(2*n+1))^2).

1, 4, 176, 79808, 372713728, 17931360207872, 8887976555024756736, 45390122553039546330628096, 2388340820825093234015277927170048, 1294826675280341699389150405743029631844352

Offset

0,2

Comments

a(n)/4^n is an integer. - Seiichi Manyama, Dec 31 2020

Links

Seiichi Manyama, Table of n, a(n) for n = 0..40

Wikipedia, Chebyshev polynomials

Wikipedia, Resultant

Formula

a(n) = 2^n * sqrt(Resultant(T_{2*n+1}(i*x/2), U_{2*n}(x/2))), where T_n(x) is a Chebyshev polynomial of the first kind, U_n(x) is a Chebyshev polynomial of the second kind and i = sqrt(-1). - Seiichi Manyama, Jan 09 2021

a(n) ~ 2^(1/8) * exp(G*(2*n + 1)^2/Pi) / (1 + sqrt(2))^(n + 1/2), where G is Catalan's constant A006752. - Vaclav Kotesovec, Mar 18 2023

Maple

for n from 0 to 12 do a[n]:=2^(2*n*n)*product(product(cos(i*Pi/(2*n+1))^2+ sin(j*Pi/(2*n+1))^2, j=1..n), i=1..n) od: seq(round(evalf(a[n], 300)), n=0..12);

Mathematica

Table[2^(2*n^2) * Product[Product[Cos[i*Pi/(2*n + 1)]^2 + Sin[j*Pi/(2*n + 1)]^2, {i, 1, n}], {j, 1, n}], {n, 0, 15}] // Round (* Vaclav Kotesovec, Mar 18 2023 *)

Prog

(PARI) default(realprecision, 120);
{a(n) = round(prod(i=1, n, prod(j=1, n, 4*cos(i*Pi/(2*n+1))^2+4*sin(j*Pi/(2*n+1))^2)))} \\ Seiichi Manyama, Dec 31 2020
(PARI) {a(n) = sqrtint(4^n*polresultant(polchebyshev(2*n+1, 1, I*x/2), polchebyshev(2*n, 2, x/2)))} \\ Seiichi Manyama, Jan 09 2021
(Python)
from math import isqrt
from sympy.abc import x
from sympy import resultant, chebyshevt, chebyshevu, I
def A127606(n): return isqrt(resultant(chebyshevt((n<<1)+1, I*x/2), chebyshevu(n<<1, x/2)))<<n if n else 1 # Chai Wah Wu, Nov 07 2023

Crossrefs

Cf. A004003, A127605.

Sequence in context: A343610 A113254 A210780 * A041945 A300387 A082393

Adjacent sequences: A127603 A127604 A127605 * A127607 A127608 A127609

Keyword

nonn

Author

Miklos Kristof, Apr 03 2007

Status

approved