OFFSET
1,2
COMMENTS
Squares resulting in leading zeros excluded.
(2*10^k-11)/9 are terms, i.e. A165402 is a subsequence. - Chai Wah Wu, Apr 23 2022
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..1838
Eric Weisstein's World of Mathematics, Square Number
EXAMPLE
13516^2 = 18268225{6} -> {6}18268225 = 24865^2.
PROG
(Python)
from itertools import count, islice
from sympy.solvers.diophantine.diophantine import diop_DN
def A045877_gen(): # generator of terms
for l in count(0):
l1, l2 = 10**(l+1), 10**l
yield from sorted(set(abs(x) for z in (diop_DN(10, m*(1-l1)) for m in range(10)) for x, y in z if l1 >= x**2 >= l2))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
EXTENSIONS
More terms from Patrick De Geest, Nov 15 1998
STATUS
approved