OFFSET
0,2
MAPLE
a := n -> add(add(tan(j*Pi/(1 + 2*k))^(2*(n - k)), j = 0..k), k = 0..n):
seq(round(evalf(a(n), 64)), n = 0..20);
PROG
(PARI)
a(n) = { sum(k=0, n, trace(matcompanion(sum(m=0, k, x^m*binomial(2*k+1, 2*(k-m))*(-1)^(m+1)))^(n-k))+((n-k)==0)) } \\ Thomas Scheuerle, Nov 11 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Nov 11 2024
STATUS
approved