题目:
https://ac.nowcoder.com/acm/problem/20223
对于一个长度为 L ≥ 2 L\ge2 L≥2的序列 X : x 1 , x 2 . . . x n X:x_1,x_2...x_n X:x1,x2...xn,如果满足对于任意 1 ≤ i < j ≤ L 1\le i<j\le L 1≤i<j≤L,均有 x i + x j x_i+x_j xi+xj不为质数,则认为序列 X X X是一个“反质数序列”。现在有一个长度为 N N N的序列 A : a 1 , a 2 . . . a n A: a_1,a_2...a_n A:a1,a2...an,他希望从中选出一个包含元素最多的子序列,使得这个子序列是一个反质数序列。
思路:
质数
p
=
p=
p=奇数
+
+
+偶数
(
p
>
2
)
(p>2)
(p>2),所以把
A
A
A数组按奇偶分成两类,变成一个二分图,边代表两个数加起来为质数,则最长序列就是最大独立集。因为最长序列最多只有一个1,所以如果原序列1的个数大于等于二,就只算一个。那么这样的话也就不需要考虑
1
+
1
=
2
1+1=2
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;
}