bzoj4857反质数序列

本题是网络流。。然而我还不想不出。。

看到数据范围，综合时间范围考虑，可以想到网络流，但是具体怎么建图，其实里面挺有学问的。

考虑整个序列内最多只有一个1的情况。这时，我们将所有和为质数的数对连边，发现一定不会出现奇环，这个可以利用奇偶性证明。

这时。。就成了二分图。。

然后，对于一个序列，如果有许多个1，那便只留下一个1即可，其他的1不会给出更好的答案，直接舍弃。

这样，跑最大独立集即可。

总结：

1、我本来想写dinic的，写完wa了2次，原因在于不知道二分图的连边方向（究竟是a连向b还是b连向a）最后直接爆炸。后来改成直接匹配就好了。。

2、建图无论如何都想不出。。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<string>
#include<cstring>
#include<cctype>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<iomanip>
#include<sstream>
#include<cstdlib>
#include<ctime>
#include<list>
#include<deque>
#include<bitset>
#include<fstream>
#define ld double
#define ull unsigned long long
#define ll long long
#define pii pair<int,int >
#define iiii pair<int,pii >
#define mp make_pair
#define INF 1000000000
#define MOD 1000000007
#define rep(i,x) for(int (i)=0;(i)<(x);(i)++)
inline int getint()
{
	int x=0,p=1;char c=getchar();
	while (c<=32)c=getchar();
	if(c==45)p=-p,c=getchar();
	while (c>32)x=x*10+c-48,c=getchar();
	return x*p;
}
using namespace std;
//
const int maxn=3000+5;
const int maxs=2e5+10;
vector<int>G[maxn],v;
queue<int>qu;
int ok[maxs],match[maxn],a[maxn],dis[maxn],iter[maxn],n;
bool vis[maxn];
//
inline bool isok(int x)
{
	if(ok[x]>=0)return ok[x];
	for(int i=2;i*i<=x;i++)if(x%i==0)return ok[x]=0;
	return ok[x]=1;
}
inline void faddedge(int x,int y)
{
	G[x].push_back(y);
	G[y].push_back(x);
}
bool dfs(int v)
{
	vis[v]=true;
	rep(i,G[v].size())
	{
		int u=G[v][i],w=match[u];
		if(w<0||!vis[w]&&dfs(w))
		{
			match[u]=v;
			match[v]=u;
			return true;
		}
	}
	return false;
}
int work()
{
	int ans=0;	
	memset(match,-1,sizeof(match));
	rep(i,n)if(match[i]<0)
	{
		memset(vis,0,sizeof(vis));
		if(dfs(i))ans++;
	}
	return ans;
} 
void finit()
{
	n=getint();
	rep(i,n)a[i]=getint();
	bool flag=false;
	rep(i,n)
	{
		if(flag&&a[i]==1)continue;
		v.push_back(a[i]);
		if(a[i]==1)flag=true;
	}
	n=v.size();
	rep(i,n)a[i]=v[i];
	memset(ok,-1,sizeof(ok));
	rep(i,n)for(int j=i+1;j<n;j++)if(isok(a[i]+a[j]))
		faddedge(i,j);	
}
int main()
{
	finit();
	printf("%d\n",n-work());
	return 0;
}