筛法求亲和数对

【输入格式】

一个整数M 

【输出格式】

输出所有小于M 的亲和数对的和

 

300<=M<=2000000

 

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;

const int NR=2000000;
int arr[NR];
//数组中存储下标i的真因数之和
int main(){
    int i,j;
    //1是每个数的真因数，因此数组中均加上1
    for(i = 1;i <= NR;i++) {
        arr[i] = 1;
    }
    
    //i代表因子，从2开始遍历，依次加到符合条件的数组中
    for(i = 2;i + i <= NR;i++){
        j = i + i;//j为数字的下标
        //依次将因子加到相应的数组值中
        //数组中存放因子之和
        while(j <= NR){
            arr[j] += i;
            j += i;
        }
    }
 
    int n;
    scanf("%d",&n);
    int sum=0;
    //遍历亲和数对，
    for(i = 2;i <= n;i++){
        //printf("%d   %d\n",i,arr[i]);
        //a，b两个数不相等
        //不能重复累加a和b，因此筛选出a<b的数对
        //b小于n
        if((arr[i] < i) && (i < n) && (arr[arr[i]] == i)) {
            //printf("----%d   %d\n",i,arr[i]);
            sum = sum + arr[i] + i;
        }
    }
    printf("%d",sum);
    return 0;
}