|
MAPLE
|
f:=proc(n) local i, j, k; i:=cat(n, n); j:=convert(i, decimal, 10); issqr(j); end;
with(numtheory): Digits:=50:for d from 1 to 22 do tendp1:=10^d+1: tendp1fact:=ifactors(tendp1)[2]: n:=mul(piecewise(tendp1fact[i][2] mod 2=1, tendp1fact[i][1], 1), i=1..nops(tendp1fact)):for i from ceil(sqrt((10^(d-1))/n)) to floor(sqrt((10^d-1)/n)) do printf("%d, ", tendp1*n*i^2) od: od: (C. Ronaldo)
|
