login
A381020
a(n) = A381018(100*n).
2
30, 45, 60, 75, 92, 106, 119, 133, 141, 157, 170, 177, 185, 204, 224, 236, 245, 260, 275, 292, 305, 318, 330, 342, 359, 371, 382, 390, 405, 419, 430, 444, 457, 472, 490, 507, 524, 535, 550, 561, 570, 583, 593, 604, 611, 621, 627, 638, 647, 659, 670, 679, 683, 697
OFFSET
1,1
LINKS
Thomas Ehrenborg, Gauss Circle Primes, arXiv:2502.06804 [math.GM], 2025. See Table 2 p. 4.
PROG
(Python)
from math import isqrt
from sympy import isprime
def A381020(n): return sum(1 for m in range(1, 100*n+1) if isprime(1+(sum(isqrt(k*((m<<1)-k)) for k in range(1, m+1))<<2))) # Chai Wah Wu, Feb 13 2025
CROSSREFS
Sequence in context: A062385 A164803 A152569 * A114944 A075290 A281919
KEYWORD
nonn
AUTHOR
Michel Marcus, Feb 12 2025
EXTENSIONS
More terms from Jinyuan Wang, Feb 13 2025
STATUS
approved