|
MATHEMATICA
|
f[n_] := Block[{l = Floor[ Log[10, n] + 1], rd = RealDigits[ Sqrt[n], 10, 24], id = IntegerDigits[n]}, rdd = Drop[ rd[[1]], rd[[2]]]; While[ rdd[[1]] == 0, rdd = Drop[rdd, 1]]; Take[rdd, l] == id]; Do[ If[ StringPosition[ ToString[ N[ Sqrt[n], 24]], ToString[n]] != {}, If[ f[n], Print[n]]], {n, 2, 6 10^8}] (from Robert G. Wilson v Aug 02 2004)
|