[LibreOJ β Round #4] 子集

 

    显然是个二分图,直接求最大独立就行了。

 

#include<bits/stdc++.h>
#define ll long long
#define pb push_back
using namespace std;
const int maxn=505;
vector<int> g[maxn];
struct lines{
	int to,flow,cap;
}l[maxn*maxn];
int t=-1,S,T,d[maxn],cur[maxn];
bool v[maxn];

inline void add(int from,int to,int cap){
	l[++t]=(lines){to,0,cap},g[from].pb(t);
	l[++t]=(lines){from,0,0},g[to].pb(t);
}

inline bool BFS(){
	queue<int> q;
	memset(v,0,sizeof(v));
	q.push(S),v[S]=1,d[S]=0;
	int x; lines e;
	
	while(!q.empty()){
		x=q.front(),q.pop();
		for(int i=g[x].size()-1;i>=0;i--){
			e=l[g[x][i]];
			if(e.flow<e.cap&&!v[e.to]){
				v[e.to]=1,d[e.to]=d[x]+1;
				q.push(e.to);
			}
		}
	}
	return v[T];
}

int dfs(int x,int A){
	if(x==T||!A) return A;
	int flow=0,f,sz=g[x].size();
	for(int &i=cur[x];i<sz;i++){
		lines &e=l[g[x][i]];
		if(d[x]==d[e.to]-1&&(f=dfs(e.to,min(e.cap-e.flow,A)))){
			A-=f,flow+=f;
			e.flow+=f,l[g[x][i]^1].flow-=f;
			if(!A) break;
		}
	}
	return flow;
}

inline int max_flow(){
	int an=0;
	while(BFS()){
		memset(cur,0,sizeof(cur));
		an+=dfs(S,1<<30);
	}
	return an;
}

ll gcd(ll x,ll y){ return y?gcd(y,x%y):x;}
ll a[505];
int n;

inline void build(){
	for(int i=1;i<=n;i++) if(a[i]&1) add(S,i,1); else add(i,T,1);
	for(int i=1;i<=n;i++) if(a[i]&1)
	    for(int j=1;j<=n;j++) if(!(a[j]&1))
	        if(gcd(a[i],a[j])==1&&gcd(a[i]+1,a[j]+1)==1) add(i,j,1);
}

int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%lld",a+i);
	S=0,T=n+1,build();
	printf("%d\n",n-max_flow());
	return 0;
}

  

posted @ 2018-04-18 10:15  蒟蒻JHY  阅读(142)  评论(0编辑  收藏  举报