#include <stdio.h>
#include <stdlib.h>
#include <math.h>
main()
{
    int i, j, k, m, n;
    for (i=1; i<=20000000; ++i)
    {
        k = i%10;
        if (k!=1 && k!=5 && k!=6)
            continue;
        k = 0;
        for (n=1,m=i,j=log10(i)+1; j>0; --j,m/=10,n*=10)
            k += (i%(int)pow(10,j)) * n*(m%10);
        if (k%n == i)
            printf("%d\n", i);
    }
}

看懂了~就是按位计算平方的部分尾数来避免数据溢出~~~

主要脑筋都是动在避免数据溢出上。为免溢出所做的动作比算答案本身还多，结果也就拖慢了。

#include <stdio.h>
main()
{
    int nn = 100000001;
    int n, c, k, m;
    for (n=1; n<nn; ++n)
    {
        c = n;
        m = 0;
        for (;;)
        {
            if (c == 0)
            {
                printf("%d\n", n);
                break;
            }
            k = c % 10;
            m += k * n;
            if (m%10 != k)
                break;
            m = m / 10;
            c = c/10;
        }
    }
}

能简单说下思路吗？谢谢

好像看出来了

能简单说下思路吗？谢谢

就是逐位乘比对个位
如：

1、m=0
   n=625
   c=625, k=c%10=5
 ---------
   k*n=625*5=3125  
   m+=k*n=0+3125=3125, m%10=5 与k比对，不等n不是自守

2、m=m/10=312
   n=625
   c=c/10=62, k=c%10=2
   -----------
   k*n=2*625=1250
   m+=k*n=312+1250=1562, m%10=2 与k比对，不等n不是自守

3、m=m/10=156
   n=625
   c=c/10=6, k=c%10=6
   -----------
   k*n=6*625=3750
   m+=k*n=156+3750=3906, m%10=6 与k比对，不等n不是自守

4、c=c/10=0 结束，每一位都对上，n是自守。

谢谢。很受益