|
MAPLE
|
with(linalg): P:=proc(q) local a, b, c, d, i, j, k, n, t;
for i from 1 to q do for n from 1 to q do
d:=ilog10(n)+1; a:=convert(n, base, 10); c:=[];
for k from 1 to nops(a) do c:=[op(c), a[-k]]; od; t:=[op([]), c];
for k from 2 to d do b:=[op([]), c[nops(c)]]; for j from 1 to nops(c)-1 do
b:=[op(b), c[j]]; od; c:=b; t:=[op(t), c]; od; if i*n=det(t) then
print(n); break; fi; od; od; end: P(10^7);
|