#include <stdio.h> //自守数算法 //ep : 25 ^ 2 = 625 76 ^ 2 = 5776 9376 ^ 2 = 87909376 /*ep : * 376 被乘数 * *376 乘数 * ------ --------- * 2256 第一个部分积=被乘数*乘数的倒数第一位 * 2632 第二个部分积=被乘数*乘数的倒数第三位 * 1125 第三个部分积=被乘数*乘数的倒数第三位 *-------- * 141376 将以上的部分积的后3位求和后截取后3位就是3位数乘积的后3位。 */ /*由number的位数确定截取数字进行乘法时的系数k*/ #define forech_bit_num(mul,number,k) \ for(mul=number,k=1;(mul/=10)>0;k*=10) ; //在0~xxxx这些数中寻找自守数 #define forech_number(number,num) \ for(number=0;number<num;number++) //自守数核心算法:(部分积+截取被乘数的后N位*截取乘数的第M位),%kk再截取部分积 #define automorphic_number(mul,number,k,ll,kk) \ mul=(mul+(number%(k*10))*(number%ll-number%(ll/10)))%kk; long print_automorphic_number(long num) { long mul,number,k,ll,kk; forech_number(number,num) { forech_bit_num(mul,number,k); kk=k*10; /*kk为截取部分积时的系数*/ mul=0; /*积的最后n位*/ ll=10; /*ll为截取乘数相应位时的系数*/ while(k > 0) { automorphic_number(mul,number , k ,ll ,kk); k/=10; /*k为截取被乘数时的系数*/ ll*=10; } if(number==mul){ /*判断若为自守数则输出*/ printf("%ld ", number); } } } int main(void) { print_automorphic_number(1000); return 0 ; }