关于自守数

wayne

其实我前面替Mathematica多想了一步，效率反而低了。下面的代码更简洁，效率更高：

1. m = 100; n = 21; sol = a /. Solve[a^n - a == 10^m k && 10^(m - 1) <= a < 10^m, {a, k},Integers]

复制代码