亲和数

一：什么是亲和数？

亲和数问题最早是由毕达哥拉斯学派发现和研究的。他们在研究数字的规律的时候发现有以下性质特点的两个数：

220的真因子是：1、2、4、5、10、11、20、22、44、55、110；
284的真因子是：1、2、4、71、142。

而这两个数恰恰等于对方的真因子各自加起来的和（sum[i]表示数i 的各个真因子的和），即
220=1+2+4+71+142=,
284=1+2+4+5+10+11+20+22+44+55+110=sum[220]。
得284的真因子之和sum[284]=220，且220的真因子之和sum[220]=284，即有sum[220]=sum[sum[284]]=284。

二：如何求500万以内的亲和数？

假如我们构建一个数组sum，数组的下标 i 和sum[i] 之间满足sum[i]=i的所有真因子之和，这个数组是不是很容易构建？

vector<int> sum(N+1,1);
    for(int i=2;i<=N/2;i++){
    for(int j=2*i;j<=N;j+=i){
<span style="white-space:pre">	</span>sum[j]+=i;
    }
}

简单的说，就是

1. 因为1是所有数的真因数，初始化sum[N]为1

2. 从2到N/2的所有数i，都可能是N内其他数的真因子，遍历一遍，给所有以sum[i]加上i。此时注意一个数不可能是它自己的真因数。

三：代码

<pre name="code" class="cpp">void f(int N){
    vector<int> sum(N+1,1);
    for(int i=2;2*i<=N;i++){
        for(int j=2*i;j<=N;j+=i){
            sum[j]+=i;
        }
    }
    for(int i=0;i<=N;i++){
        if(sum[i]<=N && sum[i]>i && sum[sum[i]]==i)
            cout<<i<<" "<<sum[i]<<endl;
    }
}