%I #15 Aug 17 2019 17:04:14
%S 3,13,31,53,59,97,127,137,139,179,241,293,313,389,457,547,563,569,647,
%T 733,769,853,857,877,881,977,1019,1117,1129,1153,1283,1291,1409,1433,
%U 1531,1543,1567,1579,1597,1699,1741,1747,1759,1889,1907,2053,2083,2221
%N Primes p such that p + ceiling(sqrt(p)) is also prime.
%H Charles R Greathouse IV, <a href="/A145907/b145907.txt">Table of n, a(n) for n = 1..10000</a>
%e For p = 13; 13 + ceiling(sqrt(13)) = 17, which is prime.
%t Select[Prime[Range[400]],PrimeQ[#+Ceiling[Sqrt[#]]]&] (* _Harvey P. Dale_, Aug 17 2019 *)
%o (PARI) isok(n) = isprime(n) && isprime(n + sqrtint(n) + 1); \\ _Michel Marcus_, Aug 15 2013
%Y Cf. A086085, A145908.
%K nonn
%O 1,1
%A _Kyle D. Balliet_, Oct 24 2008, Nov 06 2008
%E More terms from _R. J. Mathar_, Oct 31 2008