OFFSET
1,1
COMMENTS
These are the primes that are located between an oblong number and its following square number. - Charles Kusniec, Apr 17 2020
Primes in A063657. - Charles Kusniec, Sep 04 2022
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
3, seq(op(select(isprime, [$(t^2+t+1)..(t^2+2*t-1)])), t=1..31); # Robert Israel, Apr 20 2020
MATHEMATICA
Select[Range[1000], PrimeQ[#] && ContinuedFraction[Sqrt[#]][[2, 1]] == 1 &] (* Amiram Eldar, Apr 17 2020 *)
Select[Prime[Range[150]], ContinuedFraction[Sqrt[#]][[2, 1]]==1&] (* Harvey P. Dale, Mar 04 2022 *)
PROG
(PARI) isok(p) = isprime(p) && contfrac(sqrt(p))[2] == 1;
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Apr 17 2020
STATUS
approved