login
A110588
Squares of the form 2*p+3 that are squares of primes.
1
9, 25, 49, 121, 169, 361, 529, 841, 961, 1369, 1681, 2209, 4489, 5329, 6241, 9409, 10201, 10609, 11881, 22801, 24649, 38809, 44521, 51529, 54289, 57121, 66049, 69169, 85849, 94249, 109561, 113569, 121801, 124609
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
Sequence in context: A176970 A110284 A109367 * A282631 A245694 A045972
KEYWORD
nonn
AUTHOR
Walter Kehowski, Sep 13 2005
EXTENSIONS
More terms from Michel Marcus, Aug 06 2013
STATUS
approved