login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A110589 Primes p such that 2*q+3 = p^2, where q is prime. 8
3, 5, 7, 11, 13, 19, 23, 29, 31, 37, 41, 47, 67, 73, 79, 97, 101, 103, 109, 151, 157, 197, 211, 227, 233, 239, 257, 263, 293, 307, 331, 337, 349, 353, 359, 367, 389, 397, 409, 443, 449, 463, 487, 491, 509, 521, 547, 569, 571, 587, 613, 619, 653, 661, 673, 727 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
LINKS
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
Sequence in context: A109358 A216285 A192862 * A138004 A045395 A191377
KEYWORD
nonn
AUTHOR
Walter Kehowski, Sep 13 2005
EXTENSIONS
More terms from Robert G. Wilson v, Jun 19 2006
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 15:38 EDT 2024. Contains 371254 sequences. (Running on oeis4.)