OFFSET
0,1
LINKS
Robert Israel, Table of n, a(n) for n = 0..1000
FORMULA
5^n*a(n) = A358426(n)^2 - 6.
EXAMPLE
a(3) = 1031 because 5^3 * 1031 + 6 = 128881 = 359^2 where 1031 and 359 are prime, and no smaller prime works.
MAPLE
f:= proc(n) local v, a, b, k, p, q;
v:= 5^n;
a:= numtheory:-msqrt(6, v);
if a < v/2 then b:= v-a
else b:= a; a:= v-a
fi;
for k from 0 do
for q in [k*v+a, k*v+b] do
if isprime(q) then
p:= (q^2-6)/v;
if isprime(p) then return p fi;
fi
od od
end proc:
map(f, [$0..30]);
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Nov 15 2022
STATUS
approved