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

题目：
https://ac.nowcoder.com/acm/problem/20223

对于一个长度为$L\ge 2$的序列$X:x_1,x_2...x_n$，如果满足对于任意$1\le i<j\le L$，均有$x_i+x_j$不为质数，则认为序列$X$是一个“反质数序列”。现在有一个长度为$N$的序列$A:a_1,a_2...a_n$，他希望从中选出一个包含元素最多的子序列，使得这个子序列是一个反质数序列。

思路：
质数$p=$奇数$+$偶数$(p>2)$，所以把$A$数组按奇偶分成两类，变成一个二分图，边代表两个数加起来为质数，则最长序列就是最大独立集。因为最长序列最多只有一个1，所以如果原序列1的个数大于等于二，就只算一个。那么这样的话也就不需要考虑$1+1=2$这种情况。

#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;
}