OFFSET
1,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
FORMULA
a(n) is in this sequence iff A106482(a(n)) = 2.
a(n) is in this sequence iff a(n) is prime and 2*a(n)^2-1 is also prime.
a(n) = prime(A092058(n)). - R. J. Mathar, Aug 20 2019
MAPLE
q:= p-> andmap(isprime, [p, 2*p^2-1]):
select(q, [$2..2000])[]; # Alois P. Heinz, Jun 21 2022
MATHEMATICA
Select[Table[Prime[n], {n, 500}], PrimeQ[2*#^2 - 1] &] (* Ray Chandler, May 03 2005 *)
PROG
(Magma) [p: p in PrimesUpTo(2500)| IsPrime(2*p^2-1)]; // Vincenzo Librandi, Jan 29 2011
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jonathan Vos Post, May 03 2005
EXTENSIONS
Extended by Ray Chandler, May 03 2005
STATUS
approved