%I #14 Jan 27 2021 11:39:02
%S 236821687,1042169341,824576219,3936576629,358825573,14488357127,
%T 85083966119,31862900551,18816350447,50819969281,38773854101,
%U 52655522537,1040940269,379261901201,309046612187,428435406197,2040640920173,1956306006947,1467786584689,19138237171,154667461493,28186607929
%N Primes abs(A335581(k))/2^14 for k in A335582.
%H Robert Israel, <a href="/A335591/b335591.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = abs(A335581(A335582(n)))/2^14.
%e 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) =
%e 824576219.
%p N:= 100: # for a(1)..a(N)
%p count:= 0:
%p R:= NULL:
%p L:= [seq(2*i-33, i=1..16)]:
%p for k from 1 while count < N do
%p for i from 1 to 16 do
%p for x from L[i]+32 by 32 do until isprime(x);
%p L[i]:= x;
%p od;
%p v:= abs(LinearAlgebra:-Determinant(Matrix(4, 4, L)))/2^14;
%p if isprime(v) then count:= count+1; R:= R, v; fi;
%p od:
%p R;
%Y Cf. A335581, A335582.
%K nonn,look
%O 1,1
%A _J. M. Bergot_ and _Robert Israel_, Jan 26 2021