【C语言训练】自守数问题

题目链接：http://115.28.203.224/problem.php?cid=1010&pid=24

题目描述

自守数是指一个数的平方的尾数等于该数自身的自然数。
例如：
25^2=625
76^2=5776
9376^2=87909376
请求出200000以内的自守数?

输入

输出

200000以内的自守数(包括0, 数之间用两个空格分开)

样例输出

0 1 5 6 25 76 376 625 9376 90625

提示

若采用“求出一个数的平方后再截取最后相应位数”的方法显然是不可取的，因为计算机无法表示过大的整数。
分析手工方式下整数平方(乘法)的计算过程，以376为例：
376 被乘数
X 376 乘数
----------
2256 第一个部分积=被乘数*乘数的倒数第一位
2632 第二个部分积=被乘数*乘数的倒数第二位
1128 第三个部分积=被乘数*乘数的倒数第三位
----------
141376 积
本问题所关心的是积的最后三位。分析产生积的后三位的过程，可以看出，在每一次的部分积中，并不是它的每一位都会对积的后三位产生影响。总结规律可以得到：在三位数乘法中，对积的后三位产生影响的部分积分别为：
第一个部分积中：被乘数最后三位*乘数的倒数第一位
第二个部分积中：被乘数最后二位*乘数的倒数第二位
第三个部分积中：被乘数最后一位*乘数的倒数第三位
将以上的部分积的后三位求和后截取后三位就是三位数乘积的后三位。这样的规律可以推广到同样问题的不同位数乘积。
按照手工计算的过程可以设计算法编写程序。

代码

1 #include <iostream>
2 using namespace std;
3 int main()
4 {
5     int i, j, k, l, m;
6     for (i = 0; i < 200000; i++)
7     {
8         for (j = i, k = 1; (j /= 10) > 0; k *= 10);
9             l = k * 10;
10         j = 0;
11         m = 10;
12         while (k > 0)
13         {
14             j = (j + (i % (k * 10)) * (i % m - i % (m / 10))) % l;
15             k /= 10;
16             m *= 10;
17         }
18         if (i == j)
19         {
20             if (i == 0)
21                 cout<<i;
22             else
23                 cout<<" "<<i;
24         }
25     }
26     return 0;
27 }