login
Primes abs(A335581(k))/2^14 for k in A335582.
3

%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