%I #13 Jan 27 2021 18:42:41
%S 65,72,151,182,207,357,416,534,537,570,574,668,829,935,945,955,1002,
%T 1008,1014,1187,1196,1226,1239,1551,1553,1733,1864,1944,1972,1978,
%U 2018,2125,2263,2266,2310,2340,2693,2900,3167,3451,3466,3475,3486,3518,3691,3734,3739,3838,3913,3923,3969,4004
%N Numbers k such that abs(A335581(k))/2^14 is prime.
%H Robert Israel, <a href="/A335582/b335582.txt">Table of n, a(n) for n = 1..10000</a>
%e 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.
%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, k fi;
%p od:
%p R;
%Y Cf. A335581, A335591.
%K nonn
%O 1,1
%A _J. M. Bergot_ and _Robert Israel_, Jan 26 2021
|