NOIP2016提高A组集训第3场10.31高维宇宙

Posted 无尽的蓝黄

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了NOIP2016提高A组集训第3场10.31高维宇宙相关的知识,希望对你有一定的参考价值。

题解

技术分享图片

分析

因为只有奇数和偶数配对才有可能得出质数,
暴力求出每一对\(a_i+a_j\)为质数,将其中的奇数想偶数连一条边。
二分图匹配,匈牙利算法。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
const int mo=1000000007;
const int N=45;
using namespace std;
int n,a[N],pt[N*N],used[N*N],b[N][N],ans;
int find(int x)
{
    if(used[x]) return 0;
    used[x]=1;
    for(int i=1;i<=b[x][0];i++)
    {
        if(pt[b[x][i]]==0 || find(pt[b[x][i]]))
        {
            pt[b[x][i]]=x;
            return 1;
        }
    }
    return 0;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if(i!=j)
            {
                bool q=true;
                for(int k=2;k<=sqrt(a[i]+a[j])+1 && q;k++)
                {
                    if((a[i]+a[j])%k==0) q=false;
                }
                if(q) b[i][++b[i][0]]=j;
            }
    for(int i=1;i<=n;i++)
        if(a[i]%2)
        {
            memset(used,0,sizeof(used));
            ans+=find(i);
        }
    cout<<ans<<endl;
}

以上是关于NOIP2016提高A组集训第3场10.31高维宇宙的主要内容,如果未能解决你的问题,请参考以下文章

NOIP2016提高A组集训第1场10.29完美标号

NOIP2016提高A组集训第4场11.1平衡的子集

NOIP2016提高A组集训第14场11.12随机游走

NOIP2016提高A组集训第13场11.11最大匹配

2018.10.2浪在ACM 集训队第二次测试赛

NOIP2016提高A组模拟8.15Password