OFFSET
1,1
COMMENTS
Subsequence of primes of A028881. - Michel Marcus, Apr 11 2015
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..3000
Patrick De Geest, Palindromic Quasipronics of the form n(n+x)
MAPLE
A028883:=n->`if`(isprime(n^2-7), n^2-7, NULL): seq(A028883(n), n=1..500); # Wesley Ivan Hurt, Apr 11 2015
MATHEMATICA
Select[Range[3, 410]^2 - 7, PrimeQ] (* Harvey P. Dale, Sep 20 2011 *)
PROG
(Magma) [a: n in [3..500] | IsPrime(a) where a is n^2-7]; // Vincenzo Librandi, Dec 01 2011
(PARI) lista(nn) = forprime (n=1, nn, if (issquare(n+7), print1(n, ", "))) \\ Michel Marcus, Apr 11 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
More terms from Michel Marcus, Apr 11 2015
STATUS
approved