参考代码：

#include <stdio.h>
#define MAXNUMBER 1000
int main()
{
    int sum, x, i;
    for(x = 1; x <= MAXNUMBER; x++){
        sum = 0;
        for(i = 1; i <= x/2; i++)
            if(x%i == 0)
                sum += i;
        if(sum == x)
            printf("%d\n", x);
    }
    return 0;
}

参考答案中的迭代次数利用了以下结论：假设有一个正整数N，对于 N/2 < a < N，a 不可能是 n 的因子。因为可以推导出 1<N/a<2，所以迭代范围缩小为 i<=N/2。

可以对以上情形进行推广，假设有一个正整数N，N 的所有因子为 a0 a1 ... an，且 a0<a1<...<an，则有以下结论：若 a1 为 N 的因子，则 N/a1 也是 n 的因子，且 N=a0*an=a1*an-1=...。对于 an-1<i<an，i 不可能是 N 的因子；因为可以推导出 a0=N/an<N/i<N/an-1=a1，a0 和 a1 范围内不存在N的因子。所以在确定了 a 是 N 的因子之后，则可以确定 N 的另一个因子 N/a，且迭代范围可以缩小为i<N/a。

#include <stdio.h>

int perfect(int n)
{
    int s=1, t=n, i;
    for(i=2; i<t; i++)
    {
        if(n % i == 0)
        {
            t = n / i;
            s += i;
            if(t != i) s += t;
        }
    }
    if(n==s) return 1;
    else return 0;
}

int main()
{
    int i, j, t;
    for(i=1; i<1001; i++)
    {
        if(perfect(i))
        {
            printf("%d=1", i);
            t = i;
            for(j=2; j<t; j++)
            {
                if(i % j == 0)
                {
                    t = i / j;
                    printf("+%d+%d", j, t);
                }
            }
            printf("\n");
        }
    }
}