OFFSET
1,1
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..3000
PROG
(Python)
from itertools import count, islice
from sympy import sqrt_mod
def A116254_gen(): # generator of terms
for j in count(0):
b = 10**j
a = b*10+1
for k in sorted(sqrt_mod(-1, a, all_roots=True)):
m = (k**2+1)//a
if a*(b+4) <= k**2+1 < a*(a+3):
yield k-2
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Giovanni Resta, Feb 06 2006
EXTENSIONS
a(19)-a(22) from Chai Wah Wu, Feb 19 2024
STATUS
approved