自守数问题
如果一个自然数的平方数的尾部仍然为该自然数本身,则称其为自守数。 例如:
5 x 5 = 25
76 x 76 = 5776
625 x 625 = 390625
下面代码的目的是寻找出 2 千万以内的所有自守数。
注意,2 千万的平方已经超出了整数表达的最大范围,所以该程序使用了一个巧妙的方案。
如果我们仔细观察乘法的计算过程,就会发现实际上对乘积的尾数有贡献的环节,从而不用真正计算出整个乘积。
#include <stdio.h>
#include <stdlib.h>
void zishou()
{
int n;
for(n=1; n<20 * 1000 * 1000; n++)
{
int n2 = n;
int m = 0;
for(;;)
{
if(n2==0)
{
printf("%d\n", n);
break;
}
int k = n2 % 10; //取当前尾声
m += k * n;
if(m % 10 != k) break; //尾数必须相同
m = m / 10;
n2 = n2 / 10;
}
}
}
int main(int argc, char* argv[])
{
zishou();
return 0;
}