自守数

一、问题描述

  自守数是指一个数的平方的尾数等于该数自身的自然数，如25*25=625，76*76=5776.要求求出一定范围内的所有自守数。求100000以内的自守数

二、设计思路

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

三、程序流程图

 

四、代码实现

#include<stdio.h>
void main()
{
long mul,number,k,a,b;
printf("It exists following automorphic nmbers small than 100000:\n");
for (number=0;number<100000;nimber++)
{
for(mul=number,k=1;(mul/=10)>0;k*=10);
a=k*10;
mul=0;
b=10;
while(k>0)
{
mul=(mul+(number%(k*10))*(number%b-number%(b/10)))%a;
k/=10;
b*=10;
}
if(number==mul)
printf("%1d   ",number);
}
printf("\n");
}