A110588 Squares of the form 2*p+3 that are squares of primes.

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

Author

Walter Kehowski, Sep 13 2005

Keywords

• nonn

Links

• Charles R Greathouse IV, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A098828, A109358.

Programs

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;

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

Extensions

More terms from Michel Marcus, Aug 06 2013