|
MAPLE
|
g:= n -> convert(convert(n, base, 10), set) subset {2, 3, 4, 9}:
Res:= 2, 3:
extend:= proc(n, d) local B; global Res;
B:= {seq(x*10^d+n, x=[2, 3, 4, 9])};
Res:= Res, op(select(t -> g(t^2), B));
op(select(t -> g(t^2 mod 10^(d+1)), B))
end proc:
Agenda:= {2, 3, 4, 9}:
for d from 1 to 25 do Agenda:= map(extend, Agenda, d) od:
|