|
|
A003226——Automorphic numbers: m^2 ends with m.
1, 5, 6, 25, 76, 376, 625, 9376, 90625, 109376, 890625, 2890625, 7109376, 12890625, 87109376, 212890625, 787109376, 1787109376, 8212890625, 18212890625, 81787109376, 918212890625, 9918212890625,
40081787109376, 59918212890625, 259918212890625, 740081787109376, 3740081787109376, 6259918212890625, 43740081787109376, 56259918212890625, 256259918212890625, 743740081787109376,
2256259918212890625, 7743740081787109376, 92256259918212890625, 392256259918212890625, 607743740081787109376, 2607743740081787109376, 7392256259918212890625, 22607743740081787109376,
77392256259918212890625, 977392256259918212890625, 9977392256259918212890625, 19977392256259918212890625, 80022607743740081787109376, 380022607743740081787109376, ......}
Union[{1}, Table[PowerMod[16, 5^n, 10^n], {n, 27}], Table[PowerMod[5, 2^n, 10^n], {n, 99}]]——代码是长一些。但它是在 "生成" 答案。10^99, 电脑很快就出来了。
{1, 5, 6, 25, 76, 376, 625, 9376, 90625, 109376, 890625, 2890625, 7109376, 12890625, 87109376}
Select[Range[10^8], PowerMod[#, 2, 10^IntegerLength[#]] == # &]——代码是短一些。但它是在 "寻找" 答案。10^9, 电脑就罢工了。
典型的数学优化胜过工程优化的案例。 体现了:理解问题的数学结构往往能带来数量级的性能提升!
再来2个也不行。——OEIS不会给你说这些。
Select[Range[10^8], Mod[#^2, 10^Ceiling[Log10[#]]] == # &]——代码是短一些。但它是在 "寻找" 答案。10^9, 电脑就罢工了。
Select[Range[10^8], 10^IntegerExponent[#^2 - #, 10] > # &]——代码是短一些。但它是在 "寻找" 答案。10^9, 电脑就罢工了。 |
|