OFFSET
1,1
EXAMPLE
41634068//41634073 = 64524463 * 64524471, where // denotes concatenation.
PROG
(Python)
from itertools import count, islice
from sympy import sqrt_mod
def A116196_gen(): # generator of terms
for j in count(0):
b = 10**j
a = b*10+1
for k in sorted(sqrt_mod(21, a, all_roots=True)):
if a*(b-5) <= k**2-21 < a*(a-6) and k>4:
yield (k**2-21)//a
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Giovanni Resta, Feb 06 2006
STATUS
approved