OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 151 is a term because A335581(151) = det(20897, 20963, 20261, 21191; 21673, 21419, 22157, 20431; 22961, 20947, 22229, 21559; 21977, 21499, 21149, 22271) = 13509856772096 = 2^14 * 824576219 and 824576219 is prime.
MAPLE
N:= 100: # for a(1)..a(N)
count:= 0:
R:= NULL:
L:= [seq(2*i-33, i=1..16)]:
for k from 1 while count < N do
for i from 1 to 16 do
for x from L[i]+32 by 32 do until isprime(x);
L[i]:= x;
od;
v:= abs(LinearAlgebra:-Determinant(Matrix(4, 4, L)))/2^14;
if isprime(v) then count:= count+1; R:= R, k fi;
od:
R;
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Jan 26 2021
STATUS
approved