【JZOJ4832】【NOIP2016提高A组集训第3场10.31】高维宇宙
题目描述
数据范围
解法
由于大于4的素数只有可能由奇数和偶数的和得出。
所以根据数的奇偶性可以分出两类数奇数和偶数。
奇数之间不会相互匹配,偶数之间也不会相互匹配。
那么原问题转化为二分图最大匹配。
网络流即可。
代码
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<math.h>
#include<algorithm>
using namespace std;
const char* fin="prime.in";
const char* fout="prime.out";
const int inf=0x7fffffff;
const int maxn=207,maxa=2007;
int n,i,j,k,ans=0,en,ls,rs;
int a[maxn],c[maxa];
int b[maxn];
int tot=1,fi[maxn],la[maxa],va[maxa],ne[maxa],cnt[maxn],bz[maxn];
int gg(int x){
return 1+x;
}
void add_line(int a,int b,int c){
tot++;
ne[tot]=fi[a];
la[tot]=b;
va[tot]=c;
fi[a]=tot;
}
void add(int a,int b){
add_line(a,b,1);
add_line(b,a,0);
}
int gap(int v,int flow){
int i,use=0,k;
if (v==en) return flow;
for (k=fi[v];k;k=ne[k])
if (bz[v]==bz[la[k]]+1 && va[k]){
i=gap(la[k],min(flow-use,va[k]));
use+=i;
va[k]-=i;
va[k^1]+=i;
if (use==flow || bz[1]==en) return use;
}
if (!--cnt[bz[v]]) bz[1]=en;
cnt[++bz[v]]++;
return use;
}
bool p[maxa];
int main(){
freopen(fin,"r",stdin);
freopen(fout,"w",stdout);
for (i=2;i<maxa;i++){
if (!p[i]) c[++c[0]]=i;
for (j=1;j<=c[0];j++){
if (i*c[j]>=maxa) break;
p[i*c[j]]=true;
if (i%c[j]==0) break;
}
}
scanf("%d",&n);
en=2+n;
for (i=1;i<=n;i++){
scanf("%d",&a[i]);
if (a[i]%2) add(1,gg(i));
else add(gg(i),en);
}
for (i=1;i<=n;i++) for (j=1;j<=n;j++)
if (a[i]%2==1 && a[j]%2==0 && !p[a[i]+a[j]]) add(gg(i),gg(j));
cnt[0]=en;
while (bz[1]<en) ans+=gap(1,inf);
printf("%d",ans);
return 0;
}
启发
原本的想法就是网络流,但发现一次匹配会用去两个数的使用次数。
并没有注意到大于四的素数只能由偶数和奇数的和得到这个性质。
两两匹配问题考虑把原数集分为两个独立集。
这道题似乎还可以贪心,脑洞较小无法想出。
想点办法打开脑洞。