• 相亲数问题


    相亲数(Amicable Pair),又称亲和数、友爱数、友好数,指两个正整数中,彼此的全部约数之和(本身除外)与另一方相等。毕达哥拉斯曾说:“朋友是你灵魂的倩影,要像220与284一样亲密。” 例如220与284: 220的全部因数(除掉本身)相加是:1+2+4+5+10+11+20+22+44+55+110=284 284的全部因数(除掉本身)相加的和是:1+2+4+71+142=220 亲和数中可轻易推出,一方的全部约数之和与另一方的全部约数之和相等。(此叙述不可逆,不能用来判断是否为亲和数) 220的全部约数之和是:1+2+4+5+10+11+20+22+44+55+110+220 = 284+220 = 504 284的全部约数之和是:1+2+4+71+142+284 = 220+284 = 504

    问题:求3000以内的相亲数 简单方法:

    #include <iostream>
    #include <stdio.h>
    
    const int maxn=30000+3;
    
    int sum[maxn];
    //计算m除本身外全部约数之和
    int sumnum(int m){
        int n=0;
        for(int i=1;i<=m/2;i++){
            if(!(m%i)){
                n+=i;
            }
        }
        return n;
    }
    
    int main()
    {
        sum[1]=0;
        for(int i=2;i<=30000;i++){
            sum[i]=sumnum(i);
        }
        for(int i=2;i<=30000;i++){
            //防止数组越界
            if(sum[i]>30000) continue;
            //防止出现像6这样全部约数之和等于本身
            if(sum[i]!=i && sum[sum[i]]==i){
                cout<< i<<" "<<sum[i]<<endl;
                //防止重复打印输出
                sum[i]=1;
            }
        }
        return 0;
    }
    复制代码

    该方法的时间复杂度比较高。
    实际计算 m 的所有约数,只需要计算sqrt m次即可 例如:

    100
    1 * 100
    2 * 50
    4 * 25
    5 * 20
    10 * 10
    sum[100] = 1+2+50+4+25+5+20+10
    复制代码

    除了平方以外不会有重复出现
    完整程序

    #include <iostream>
    #include <stdio.h>
    
    const int maxn=30000+3;
    
    int sum[maxn];
    //计算m除本身外全部约数之和
    int sumnum(int m){
        //自2开始每个数都有一个为1的约数
        int n=1;
        int i;
        for(i=2;i*i<m;i++){
            if(!(m%i)){
                n+=i;
                n+=(m/i);
            }
        }
        //如果m为平方数,则其开方只能使用一次
        if(i*i==m){
            n+=i;
        }
        return n;
    }
    
    int main()
    {
        sum[1]=0;
        for(int i=2;i<=30000;i++){
            sum[i]=sumnum(i);
        }
        for(int i=2;i<=30000;i++){
            //防止数组越界
            if(sum[i]>30000) continue;
            //防止出现像6这样全部约数之和等于本身
            if(sum[i]!=i && sum[sum[i]]==i){
                cout<< i<<" "<<sum[i]<<endl;
                //防止重复打印输出
                sum[i]=1;
            }
        }
        return 0;
    }
    复制代码

    转载于:https://juejin.im/post/5aadc9975188255579186fcf

  • 相关阅读:
    POJ 3713 Transferring Sylla【Tarjan求割点】
    Tarjan算法【强连通分量】
    POJ 1273 Drainage Ditches【最大流模版】
    CDOJ 1960 构造哈密顿路径
    HDU 1384 Intervals【差分约束-SPFA】
    POJ 1364 / HDU 3666 【差分约束-SPFA】
    【SPFA与Dijkstra的对比】CDOJ 1961 咸鱼睡觉觉【差分约束-负权最短路径SPFA】
    CDOJ 1965 连通域统计【DFS】
    CDOJ 1964 命运石之门【最短路径Dijkstra/BFS】
    最小生成树模板【kruskal & prim】
  • 原文地址:https://www.cnblogs.com/twodog/p/12137335.html
Copyright © 2020-2023  润新知