[JSOI2016]反质数序列 二分图匹配

[JSOI2016]

•  奇数+奇数一定不是质数(1+1除外），偶数+偶数一定不是质数，质数只可能出现在偶数+奇数中
• 把所有的点排成两列，权值为奇数的点在左边，权值为偶数的在右边
• 如果左边的点x+右边的点y是质数，我们就连一条x->y的边
• 最后答案显然是最大独立集=n-最小点覆盖=n-最大匹配数
• 最大匹配数=最大流，Dinic即可
• 由于1比较特殊，考虑到最终答案1的出现次数<=1，所以如果有多个1只保留一个即可

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define per(i,a,b) for(int i=(a);i>=(b);i--)
using namespace std;
const int inf=1e9;
const int N=1e5+100;
struct node{int y,n,v;}e[N*20];
int prime[N],lin[N*20],d[N],v[N<<1],a[N],maxflow,len=1,S,T,n;
inline int read(){
    int num=0;char ch=getchar();
    while(!isdigit(ch))ch=getchar();
    while(isdigit(ch))num=num*10+ch-'0',ch=getchar();
    return num;
}
inline void Print(int x){if(x>9)Print(x/10);putchar(x%10+'0');}
void euler(int n){
    int m=0;
    rep(i,2,n){
        if(!v[i]){
            v[i]=i;
            prime[++m]=i;
        }
        rep(j,1,m){
            if(prime[j]>v[i]||i*prime[j]>n)break;
            v[i*prime[j]]=prime[j];
        }
    }
    memset(v,0,sizeof(v));
    rep(i,1,m)v[prime[i]]=1;
}
void LK(int x,int y,int v)
{e[++len].y=y,e[len].n=lin[x],e[len].v=v,lin[x]=len;}
void add(int x,int y,int v)
{LK(x,y,v),LK(y,x,0);}
bool bfs(){
	queue<int> q;
	memset(d,0,sizeof(d));
	q.push(S),d[S]=1;
	while(q.size()){
		int x=q.front();q.pop();
		for(int i=lin[x];i;i=e[i].n){
			int y=e[i].y;
			if(d[y]||!e[i].v)continue;
			d[y]=d[x]+1;
			q.push(y);
		}
	}return d[T]>0;
}
int dfs(int x,int minf){
	int flow=0,sum=0;
	if(x==T)return minf;
	for(int i=lin[x];i;i=e[i].n){
		int y=e[i].y;
		if(!e[i].v||d[y]!=d[x]+1)continue;
		flow=dfs(y,min(minf,e[i].v));
		if(!flow)d[y]=0;
		minf-=flow,sum+=flow,e[i].v-=flow,e[i^1].v+=flow;
		if(!minf)return sum;
	}return sum;
}
void Dinic(){
	while(bfs())maxflow+=dfs(S,inf);
}
int main(){
    //freopen("a.in","r",stdin);
    euler((N<<1)-1);
    n=read(); S=n+1,T=S+1;
    int flag=0;
    rep(i,1,n){
        a[i]=read();
        if(a[i]==1){
            if(flag) i--,n--;
            else flag=1;
        }
    }
    rep(i,1,n){
        if(a[i]&1)add(S,i,1);
        else add(i,T,1);
    }
 
    rep(i,1,n){
        if(a[i]&1)rep(j,1,n)if((a[j]&1)==0){
            if(v[a[i]+a[j]])add(i,j,1);
        }
    }
    Dinic();
    //cout<<maxflow<<"\n";
    Print(n-maxflow);
    return 0;
}