# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a335591
Showing 1-1 of 1

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

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE