OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
FORMULA
(a(n)^2-1)/8-a(n) = A165567(n).
EXAMPLE
For p=17, (17^2-1)/8-17=19 is prime, which adds p=17 to the sequence.
For p=23, (23^2-1)/8-23=43 is prime, which adds p=23 to the sequence.
MATHEMATICA
Select[Prime[Range[5, 1000]], PrimeQ[(#^2 - 1) / 8 - #]&] (* Vincenzo Librandi, Apr 10 2013 *)
PROG
(Magma) [p: p in PrimesInInterval(7, 1500) | IsPrime(((p^2-1) div 8) - p)]; // Vincenzo Librandi, Apr 10 2013
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Sep 22 2009
EXTENSIONS
Edited and extended by R. J. Mathar, Sep 25 2009
STATUS
approved