OFFSET
1,1
COMMENTS
Conjecture: a(n) exists for all n.
MAPLE
f:= proc(n) local j, p, q, r, s;
p:= ithprime(n);
for j from 1 do
q:= j*2*p^3-1; if not isprime(q) then next fi;
r:= j*p*2*q^2-1; if not isprime(r) then next fi;
s:= j*p*2*r^2-1; if not isprime(s) then next fi;
if isprime(j*p*2*s^2-1) then return j fi;
od
end proc;
map(f, [$1..25]); # Robert Israel, May 15 2025
MATHEMATICA
a[n_] := For[j = 1, j < 10^7, j++, p = Prime[n]; If[PrimeQ[q = j*2*p^3 - 1] && PrimeQ[r = j*2*p*q^2 - 1] && PrimeQ[s = j*2*p*r^2 - 1] && PrimeQ[j*2*p*s^2 - 1], Return[j]]]; Table[Print[an = a[n]]; an, {n, 1, 24}] (* Jean-François Alcover, Apr 12 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Apr 12 2013
EXTENSIONS
More terms from Jean-François Alcover, Apr 12 2013
Name clarified and more terms from Robert Israel, May 15 2025
STATUS
approved
