[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; }
我爱学习,学习使我快乐