C语言自守数

#include <stdio.h>
#include <math.h>

int main(){
   int i,a,k,m;
   for(i=1;i<10000;i++)
   {
       a=i;        
       a*=a;      //求出a的平方即i的平方
       k=log10(i)+1;//求出i的位数k
       for(m=1;k;k--)
           m*=10;
       a%=m;//通过取余求出尾数
       if(a==i)
           printf("%d\n",i);
   }
   return 0;
}