OFFSET
1,1
COMMENTS
Subsequence of A081053. All terms are even except for a(2) = 11.
EXAMPLE
a(1) = 6 as 6 = 2 + 2^2.
a(2) = 11 as 11 = 7 + 2^2 = 2 + 3^2.
a(3) = 56 as 56 = 47 + 3^2 = 31 + 5^2 = 7 + 7^2.
a(4) = 176 as 176 = 167 + 3^2 = 151 + 5^2 = 127 + 7^2 = 7 + 13^2.
a(5) = 188 as 188 = 179 + 3^2 = 163 + 5^2 = 139 + 7^2 = 67 + 11^2 = 19 + 13^2.
a(6) = 362 as 362 = 353 + 3^2 = 337 + 5^2 = 313 + 7^2 = 241 + 11^2 = 193 + 13^2 = 73 + 17^2.
PROG
(Python)
from itertools import count
from math import isqrt
from sympy import isprime, primerange
def A381333(n):
for m in count(1):
c = 0
for p in primerange(isqrt(m)+1):
if isprime(m-p**2):
c += 1
if c>=n:
return m
(PARI) f(k) = my(nb=0); forprime(p=2, sqrtint(k), if (isprime(k-p^2), nb++); ); nb;
a(n) = my(k=1); while (f(k) < n, k++); k; \\ Michel Marcus, Feb 21 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Chai Wah Wu, Feb 20 2025
STATUS
approved