OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MAPLE
ispower := proc(n, b) andmap(proc(w) evalb(w[2] mod b = 0) end, ifactors(n)[2]) end: a:=2: SQRTP||a:=[]; for z from 1 to 1 do for n from 1 to 10000 do p:=ithprime(n); m:=a*p+a+1; if ispower(m, 2) and isprime(sqrt(m)) then SQRTP||a:=[op(SQRTP||a), sqrt(m)] fi od; od; SQRTP||a;
MATHEMATICA
fQ[n_] := PrimeQ[(n^2 - 3)/2]; Select[ Prime@ Range@129, fQ@# &] (* Robert G. Wilson v, Jun 19 2006 *)
Select[Table[Sqrt[2 Prime[n] + 3], {n, 1, 30000}], PrimeQ] (* Vincenzo Librandi, Mar 19 2015 *)
PROG
(Magma) [p: p in PrimesUpTo(1000) | IsPrime((p^2-3) div 2)]; // Vincenzo Librandi, Mar 19 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Walter Kehowski, Sep 13 2005
EXTENSIONS
More terms from Robert G. Wilson v, Jun 19 2006
STATUS
approved