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