login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A335582 Numbers k such that abs(A335581(k))/2^14 is prime. 3
65, 72, 151, 182, 207, 357, 416, 534, 537, 570, 574, 668, 829, 935, 945, 955, 1002, 1008, 1014, 1187, 1196, 1226, 1239, 1551, 1553, 1733, 1864, 1944, 1972, 1978, 2018, 2125, 2263, 2266, 2310, 2340, 2693, 2900, 3167, 3451, 3466, 3475, 3486, 3518, 3691, 3734, 3739, 3838, 3913, 3923, 3969, 4004 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
LINKS
EXAMPLE
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.
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, k fi;
od:
R;
CROSSREFS
Sequence in context: A173379 A095535 A095523 * A282113 A060877 A113688
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Jan 26 2021
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 9 18:10 EDT 2024. Contains 372354 sequences. (Running on oeis4.)