OFFSET
1,1
COMMENTS
Squares of the form (10^(d+1)+1)*x + 10^d where 10^(d-1) <= x < 10^d.
EXAMPLE
a(3) = 12722025112722025 is a term because it is 112791955^2 and is the concatenation of 12722025, 1 and 12722025.
MAPLE
Res:= NULL: count:= 0:
for d from 2 to 14 do
S:= [isolve(10^d + x*(10^(d+1)+1) = y^2)];
S:= map(t -> subs(subs(_Z1=0, t), x), S);
S:= sort(select(t -> t >= 10^(d-1) and t < 10^d, S));
for s in S do
t:= (10^(d+1)+1)*s+10^d; Res:= Res, t; count:= count+1;
od
od:
Res;
PROG
(Python)
from math import isqrt
from itertools import count, islice
def b(n): return 10**len(str(n))*(10*n+1) + n # A392225
def agen(): # generator of terms
yield from (filter(lambda x: isqrt(x)**2 == x, (b(n) for n in count(1) if n%10 in {0, 1, 4, 5, 6, 9})))
g = agen()
[print(next(g), end=", ") for n in count(1)] # Michael S. Branicky, Jan 04 2026
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Robert Israel, Jan 03 2026
STATUS
approved
