BZOJ 3158 千钧一发 最小割

分析：

　　偶数对满足条件2，所有奇数对满足条件1。

　　如果你能一眼看出这个规律，这道题就完成了一半。

　　我们只需要将数分为两类，a值为奇数，就从S向这个点连容量为b值的边，a值为偶数，就从这个点向T连容量为b值的边。

　　暴力枚举，对于奇集合和偶集合中不能共存的两个数，连容量为无穷大的边。

　　求出最小割，代表这个割集要被我们舍弃。

　　然后直接用b值总和减去最小割就好。

代码：

 #include<bits/stdc++.h>
 #define ms(a,x) memset(a,x,sizeof(a))
 #define pf(x) (1LL*x*x)
 #define ll long long
 using namespace std;
 const int N=,M=,inf=0x3f3f3f3f;
 struct node{int y,z,nxt;}e[M*];int S,T,tot=;
 int d[N],h[N],c=,a[N],b[N],q[N],n,m;ll s[N];
 int gcd(int x,int y){
     return y?gcd(y,x%y):x;
 } bool pd(ll x){
     ll y=sqrt(x);return y*y!=x;
 } void add(int x,int y,int z){
     e[++c]=(node){y,z,h[x]};h[x]=c;
     e[++c]=(node){x,,h[y]};h[y]=c;
 } bool bfs(){
     int f=,t=;ms(d,-);
     q[++t]=S;d[S]=;
     while(f<=t){
         int x=q[f++];
         for(int i=h[x],y;i;i=e[i].nxt)
         if(d[y=e[i].y]==-&&e[i].z)
         d[y]=d[x]+,q[++t]=y;
     } return (d[T]!=-);
 } int dfs(int x,int f){
     if(x==T) return f;int w,tmp=;
     for(int i=h[x],y;i;i=e[i].nxt)
     if(d[y=e[i].y]==d[x]+&&e[i].z){
         w=dfs(y,min(e[i].z,f-tmp));
         if(!w) d[y]=-;
         e[i].z-=w;e[i^].z+=w;tmp+=w;
         if(tmp==f) return f;
     } return tmp;
 } void dinic(){
     while(bfs()) tot+=dfs(S,inf);
 } signed main(){
     scanf("%d",&n);int ans=;S=;T=n+;
     for(int i=;i<=n;i++)
     scanf("%d",&a[i]),s[i]=pf(a[i]);
     for(int i=;i<=n;i++)
     scanf("%d",&b[i]),ans+=b[i];
     for(int i=;i<=n;i++)
     for(int j=i+;j<=n;j++)
     if(gcd(a[i],a[j])==&&!pd(s[i]+s[j]))
     a[i]&?add(i,j,inf):add(j,i,inf);
     for(int i=;i<=n;i++)
     a[i]&?add(S,i,b[i]):add(i,T,b[i]);
     dinic();printf("%d\n",ans-tot);
     return ;
 }

最小割