OFFSET
1,1
MATHEMATICA
a[n_]:=(Module[{k=1}, While[!PrimeQ[m=2k^2+2k+1]||IntegerLength[m]<n, k++]; m]-10^(n-1)); Array[a, 15]
PROG
(Python)
from math import isqrt
from itertools import count
from sympy import isprime
def A376993(n):
for k in count(isqrt(((a:=10**(n-1))<<1)-1>>2)):
m = 2*k*(k+1)+1
if m >= a and isprime(m):
return m-a # Chai Wah Wu, Oct 13 2024
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Stefano Spezia, Oct 11 2024
EXTENSIONS
a(21)-a(26) from Chai Wah Wu, Oct 13 2024
STATUS
approved