OFFSET
1,1
COMMENTS
Terms are of the form 10*k^2 + t, with gcd(t, 10) = 1, i.e., {1, 3, 7, 9}.
Sum_{n>=1} 1/a(n) = 0.403068... converges.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
For n=9, a(n)=163 is a term because 16 left 3 is square 4^2=16.
For n=14, a(n)=491 is a term because 49 left 1 is square 7^2=49.
MATHEMATICA
Select[Prime@ Range[5, PrimePi[10^4]], IntegerQ@ Sqrt@ Floor[#/10] &] (* or *)
Select[Union@ Flatten@ Map[10 Range[31]^2 + # &, {1, 3, 7, 9}], PrimeQ] (* Michael De Vlieger, Jul 14 2016 *)
PROG
(PARI) for(n=1, 1e3, forstep(p=10*n^2+1, 10*n^2+9, [2, 4, 2], if(isprime(p), print1(p", ")))) \\ Charles R Greathouse IV, Jul 15 2016
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Dimitris Valianatos, Jul 14 2016
STATUS
approved