华为机试HJ28:素数伴侣

简介: 华为机试HJ28:素数伴侣

题目描述:

若两个正整数的和为素数,则这两个正整数称之为“素数伴侣”,如2和5、6和13,它们能应用于通信加密。现在密码学会请你设计一个程序,从已有的N(N为偶数)个正整数中挑选出若干对组成“素数伴侣”,挑选方案多种多样,例如有4个正整数:2,5,6,13,如果将5和6分为一组中只能得到一组“素数伴侣”,而将2和5、6和13编组将得到两组“素数伴侣”,能组成“素数伴侣”最多的方案称为“最佳方案”,当然密码学会希望你寻找出“最佳方案”。

输入描述:

输入说明

1 输入一个正偶数n

2 输入n个整数

注意:数据可能有多组

输出描述:

求得的“最佳方案”组成“素数伴侣”的对数。

示例:

输入:

4

2 5 6 13

2

3 6

输出:

2

0

解题思路:

本题是一个匈牙利算法的标准应用题。将输入的奇数偶数分两组,建立二分图,两两相组合能得到素数的就赋1;再建立全局的连接表modd和meven,记录该奇数(偶数)连接的偶数(奇数);然后遍历每个奇数,每个奇数在寻找合适的偶数时,建立一个visit记录在该轮查找中已经搜过的偶数,匈牙利的原则就是“先到先得,能让就让,不能让就拜拜”,所以在确定合适偶数对象时,如果该偶数没有被配对过,直接绑定,若它被配对了,那就让配对它的那个奇数再去找其他的偶数进行配对,以此类推进行递归,若到了某个位置发现找到了,则这条新的配对方案建立,比之前的配对个数多了1,若没找到,则保持原先的方案。另外因为有多组输入,别忘了清空容器和图。


我认为讲解匈牙利算法比较好的一篇文章:


匈牙利算法-看这篇绝对就够了!_夜阑听风-CSDN博客_匈牙利算法


供大家学习,一同进步。

测试代码:

#include <iostream>
#include <vector>
using namespace std;
int t[100][100]={0};
vector<int> modd(100,-1);
vector<int> meven(100,-1);
vector<int> visit(100,0);
bool IsPrimer(int n)
{
    for(int i=2;i*i<=n;i++)
    {
        if(n%i==0)
            return false;
    }
    return true;
}
int dfs(int i)
{
    for(int j=0;j<100;++j)
    {
        if(t[i][j]==1&&!visit[j])
        {
            visit[j]=1;
            if(meven[j]==-1||dfs(meven[j]))
            {
                modd[i]=j;
                meven[j]=i;
                return 1;
            }
        }
    }
    return 0;
}
int main()
{
    int number;
    vector<int> even,odd;
    while(cin>>number)
    {
        for(int i=0;i<number;++i)
        {
            int temp;cin>>temp;
            if(temp%2!=0)
                odd.push_back(temp);
            else
                even.push_back(temp);
        }
        for(int i=0;i<odd.size();++i)
        {
            for(int j=0;j<even.size();++j)
            {
                if(IsPrimer(odd[i]+even[j]))
                    t[i][j]=1;
            }
        }
        // 匈牙利算法
        int res=0;
        for(int i=0;i<odd.size();++i)
        {
            if(modd[i]==-1)
            {
                for(int j=0;j<visit.size();++j)
                {
                    visit[j]=0;
                }
                res+=dfs(i);
            }
        }
        cout<<res<<endl;
        even.clear();
        odd.clear();
        for(int i=0;i<100;i++)
        {
            for(int j=0;j<100;j++)
            {
                t[i][j]=0;
            }
        }
        for(int i=0;i<100;i++)
        {
            meven[i]=-1;
            modd[i]=-1;
        }
    }
    return 0;
}


相关文章
|
2月前
|
存储
每日一题啦(● ̄(エ) ̄●)(尼克切斯定理,等差数列)
每日一题啦(● ̄(エ) ̄●)(尼克切斯定理,等差数列)
10 0
|
7月前
|
算法
华为机试HJ76:尼科彻斯定理
华为机试HJ76:尼科彻斯定理
|
7月前
华为机试HJ88:扑克牌大小
华为机试HJ88:扑克牌大小
|
7月前
|
机器学习/深度学习 测试技术
华为机试HJ55:挑7
华为机试HJ55:挑7
|
7月前
华为机试HJ105:记负均正II
华为机试HJ105:记负均正II
|
7月前
|
Serverless 测试技术
华为机试HJ97:记负均正
华为机试HJ97:记负均正
100 1
|
7月前
|
测试技术
华为机试HJ77:火车进站
华为机试HJ77:火车进站
|
7月前
|
机器学习/深度学习
华为机试HJ35:蛇形矩阵
华为机试HJ35:蛇形矩阵
|
7月前
华为机试HJ100:等差数列
华为机试HJ100:等差数列
|
7月前
|
算法
华为机试HJ108:求最小公倍数
华为机试HJ108:求最小公倍数