OFFSET
0,2
COMMENTS
MATHEMATICA
c = 0; k = 1; Do[ While[k <= 10^n, If[ PrimeQ@ Round@ Sqrt[k(k + 1)(k + 2)(k + 3) + 1], c++ ]; k++ ]; Print@c, {n, 0, 9}] (* Robert G. Wilson v, May 01 2006 *)
PROG
(PARI) g(n)= { for(j=0, n, c=0; for(x=0, 10^j, y=round(sqrt(x*(x+1)*(x+2)*(x+3)+1)); if(ispseudoprime(y), c++)); print1(c", ") ) }
(Python)
from sympy import isprime
def A115366(n): return sum(1 for k in range(1, 10**n+1) if isprime(k*(k+3)+1)) # Chai Wah Wu, Jun 19 2024
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Cino Hilliard, Mar 07 2006
EXTENSIONS
Edited by Don Reble, Apr 24 2006
a(10) from Michael S. Branicky, Jun 20 2024
a(11) from Michael S. Branicky, Jun 23 2024
STATUS
approved