OFFSET
1,1
LINKS
Charles R Greathouse IV, 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: PW||a:=[]; for z from 1 to 1 do for n from 1 to 1000 do p:=ithprime(n); m:=a*p+a+1; if ispower(m, 2) and isprime(sqrt(m)) then PW||a:=[op(PW||a), m] fi od; od; PW||a;
PROG
(PARI) isok(n) = issquare(n) && isprime(sqrtint(n)) && (type(p=(n-3)/2) == "t_INT") && isprime(p) \\ Michel Marcus, Aug 06 2013
(PARI) v=List(); forprime(p=2, 1e4, if(isprime(p^2\2-1), listput(v, p^2))); Vec(v) \\ Charles R Greathouse IV, Aug 06 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Walter Kehowski, Sep 13 2005
EXTENSIONS
More terms from Michel Marcus, Aug 06 2013
STATUS
approved