OFFSET
1,1
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
PROG
(Python)
from itertools import count, islice
from sympy import sqrt_mod
def A116267_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)):
if a*(b+3) <= k**2-1 < a*(a+2):
yield k-2
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Giovanni Resta, Feb 06 2006
EXTENSIONS
a(29)-a(32) from Chai Wah Wu, Feb 19 2024
STATUS
approved