OFFSET
1,1
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..500
FORMULA
a(n) = prime(A064370(n+1)). - Michel Marcus, Jul 11 2014
EXAMPLE
37 is in the sequence because primepi(37) = 12, and 37 - 12 = 5^2.
541 is in the sequence because primepi(541) = 100, and 541 - 100 = 21^2.
547 is not in the sequence because primepi(547) = 101, and 547 - 101 = 446, which is not a perfect square.
MAPLE
with(numtheory): A245061:=n->`if`(type(sqrt(n-pi(n)), integer) and type(n, prime), n, NULL): seq(A245061(n), n=2..10^5); # Wesley Ivan Hurt, Jul 10 2014
MATHEMATICA
Select[Prime[Range[200]], IntegerQ[Sqrt[# - PrimePi[#]]] &] (* Alonso del Arte, Jul 11 2014 *)
PROG
(PARI) select(p->issquare(p-primepi(p)), primes(15000)) \\ Michel Marcus, Jul 11 2014
(Python)
import sympy, gmpy2
[sympy.prime(n) for n in range(1, 10**6) if gmpy2.is_square(sympy.prime(n)-n)] # Chai Wah Wu, Jul 11 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Chai Wah Wu, Jul 10 2014
EXTENSIONS
More terms from Michel Marcus, Jul 11 2014
STATUS
approved