[go: up one dir, main page]

login
A335591
Primes abs(A335581(k))/2^14 for k in A335582.
3
236821687, 1042169341, 824576219, 3936576629, 358825573, 14488357127, 85083966119, 31862900551, 18816350447, 50819969281, 38773854101, 52655522537, 1040940269, 379261901201, 309046612187, 428435406197, 2040640920173, 1956306006947, 1467786584689, 19138237171, 154667461493, 28186607929
OFFSET
1,1
LINKS
FORMULA
a(n) = abs(A335581(A335582(n)))/2^14.
EXAMPLE
A335582(3) = 151, A335581(151) = det(20897, 20963, 20261, 21191; 21673, 21419, 22157, 20431; 22961, 20947, 22229, 21559; 21977, 21499, 21149, 22271) = 13509856772096 = 2^14 * 824576219, so a(3) =
824576219.
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, v; fi;
od:
R;
CROSSREFS
Sequence in context: A029831 A289981 A203885 * A234058 A353554 A064869
KEYWORD
nonn,look
AUTHOR
J. M. Bergot and Robert Israel, Jan 26 2021
STATUS
approved