|
|
A335582
|
|
Numbers k such that abs(A335581(k))/2^14 is prime.
|
|
3
|
|
|
65, 72, 151, 182, 207, 357, 416, 534, 537, 570, 574, 668, 829, 935, 945, 955, 1002, 1008, 1014, 1187, 1196, 1226, 1239, 1551, 1553, 1733, 1864, 1944, 1972, 1978, 2018, 2125, 2263, 2266, 2310, 2340, 2693, 2900, 3167, 3451, 3466, 3475, 3486, 3518, 3691, 3734, 3739, 3838, 3913, 3923, 3969, 4004
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
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
|
|
|
STATUS
|
approved
|
|
|
|