cp's OEIS Frontend

A337146 Numbers k such that abs(A337145(k))/8 is prime.

%I A337146 #11 Jan 27 2021 22:37:02
%S A337146 1,13,41,50,53,62,67,76,89,98,108,113,137,180,211,225,236,240,250,281,
%T A337146 293,300,303,308,355,362,384,392,393,400,414,425,434,458,468,477,489,
%U A337146 525,588,589,593,625,653,662,664,671,673,674,696,698,732,758,765,795,800,819,831,851,880,916,933,938
%N A337146 Numbers k such that abs(A337145(k))/8 is prime.
%H A337146 Robert Israel, <a href="/A337146/b337146.txt">Table of n, a(n) for n = 1..10000</a>
%e A337146 a(3) = 41 is a term because A337145(41) = det(1097, 883; 877, 967) = 286408 = 8*35801 and 35801 is prime.
%p A337146 R:= NULL:
%p A337146 count:= 0:
%p A337146 L:= [-7, -5, -3, -1]:
%p A337146 for k from 1 while count < 100 do
%p A337146   for i from 1 to 4 do
%p A337146    for x from L[i]+8 by 8 do until isprime(x);
%p A337146    L[i]:= x;
%p A337146   od;
%p A337146   v:= abs(L[1]*L[4]-L[2]*L[3])/8;
%p A337146   if isprime(v) then count:= count+1; R:= R, k; fi
%p A337146 od:
%p A337146 R;
%Y A337146 Cf. A335593, A337145, A337147.
%K A337146 nonn
%O A337146 1,2
%A A337146 _J. M. Bergot_ and _Robert Israel_, Jan 27 2021