OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000
MAPLE
comps:= remove(isprime, [$4 .. 1000]):
f:= proc(n) local M;
M:= Matrix(n, n, (i, j) -> `if`(i=j, 1, comps[abs(i-j)]));
LinearAlgebra:-Determinant(M)
end proc:
map(f, [$1..25]); # Robert Israel, Dec 03 2024
MATHEMATICA
Composite[n_Integer] := FixedPoint[n + PrimePi[ # ] + 1 &, n + PrimePi[n] + 1]; f[n_] := Det[ Table[ If[i == j, 1, Composite[ Abs[i - j]]], {i, 1, n}, {j, 1, n}]]; Table[ f[n], {n, 1, 20}]
CROSSREFS
KEYWORD
sign
AUTHOR
Robert G. Wilson v, May 26 2002
EXTENSIONS
a(21)-a(23) from Stefano Spezia, Jun 27 2024
STATUS
approved