OFFSET
0,2
LINKS
Robert Israel, Table of n, a(n) for n = 0..400
FORMULA
It appears that for n > 1, a(n) / a(n-1) = 2 * (4*e(n)-1) / (2*e(n)-1), where e(n) = A006519(n).
MAPLE
f:= proc(n) uses LinearAlgebra, MmaTranslator[Mma];
sqrt(Determinant(Matrix(2*n, 2*n, (i, j) -> signum(i-j)*BitXor(i, j))))
end proc:
f(0):= 1:
map(f, [$0..30]); # Robert Israel, Mar 31 2026
MATHEMATICA
a[0] = 1;
a[n_] := Sqrt@Det@Table[Sign[i - j] BitXor[i, j], {i, 2 n}, {j, 2 n}];
Table[a[n], {n, 0, 20}]
PROG
(PARI) a(n) = sqrtint(matdet(matrix(2*n, 2*n, i, j, sign(i-j)*bitxor(i, j)))); \\ Michel Marcus, Dec 28 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrei Zabolotskii, Dec 26 2022
STATUS
approved
