BZOJ 1458 士兵占领

http://www.lydsy.com/JudgeOnline/problem.php?id=1458

题意:n x m的棋盘，k个位置不能放，每行和每列都有要求至少的士兵，求能否有最少的满足条件的士兵放法是多少。

思路:先全放满求能否满足，再尽量删掉士兵:

对于每行:能放m[i],至少放c[i],就从S连向i:m[i]-c[i],代表能删的最大

对于每列:能放m[i],至少放c[i],就从i+n连向T:m[i]-c[i],代表能删的最大

这样对于i,j这个位置如果可以放士兵，那就从i行连向j+n，流量为1，共享删掉1个

求最大流，就等于最小割，得到的就是能删的最多，然后答案就是Σm1[i]-最大流

 #include<algorithm>
 #include<cstdio>
 #include<cmath>
 #include<cstring>
 #include<iostream>
 #define inf 0x7fffffff
 int tot,go[],next[],flow[],first[];
 int op[],nodes,T,S,dis[],cnt[],pd[][],n,m,k;
 int m1[],m2[],c1[],c2[];
 int read(){
     int t=,f=;char ch=getchar();
     while (ch<''||ch>''){if (ch=='-') f=-;ch=getchar();}
     while (''<=ch&&ch<=''){t=t*+ch-'';ch=getchar();}
     return t*f;
 }
 void insert(int x,int y,int z){
     tot++;
     go[tot]=y;
     next[tot]=first[x];
     first[x]=tot;
     flow[tot]=z;
 }
 void add(int x,int y,int z){
     insert(x,y,z);op[tot]=tot+;
     insert(y,x,);op[tot]=tot-;
 }
 int dfs(int x,int f){
     if (x==T) return f;
     int mn=nodes,sum=;
     for (int i=first[x];i;i=next[i]){
         int pur=go[i];
         if (flow[i]&&dis[pur]+==dis[x]){
             int save=dfs(pur,std::min(f-sum,flow[i]));
             flow[i]-=save;
             flow[op[i]]+=save;
             sum+=save;
             if (sum==f||dis[S]>=nodes) return sum;
         }
         if (flow[i]) mn=std::min(mn,dis[pur]);
     }
     if (sum==){
         cnt[dis[x]]--;
         if (cnt[dis[x]]==){
             dis[S]=nodes;
         }else{
             dis[x]=mn+;
             cnt[dis[x]]++;
         }
     }
     return sum;
 }
 int main(){
     n=read();m=read();k=read();
     S=,T=n+m+;nodes=T+;
     int sum1=,sum2=;
     for (int i=;i<=n;i++){
         int x=read();
         sum1+=x;
         m1[i]=x;
     }
     for (int i=;i<=m;i++){
         int x=read();
         sum2+=x;
         m2[i]=x;
     }
     for (int i=;i<=k;i++){
         int x=read(),y=read();
         pd[x][y]=;
     }
     int ans=;
     for (int i=;i<=n;i++)
      for (int j=;j<=m;j++)
       if (!pd[i][j])
        c1[i]++,c2[j]++;
     int sx=;
     for (int i=;i<=n;i++)
      sx+=c1[i];
     for (int i=;i<=n;i++)
      if (c1[i]<m1[i]) {
             printf("JIONG!\n");return ;
      }
      for (int i=;i<=m;i++)
       if (c2[i]<m2[i]){
             printf("JIONG!\n");return ;
       }
     for (int i=;i<=n;i++)
      add(S,i,c1[i]-m1[i]);
     for (int i=;i<=m;i++)
      add(i+n,T,c2[i]-m2[i]);
     for (int i=;i<=n;i++)
      for (int j=;j<=m;j++)
       if (!pd[i][j])
        add(i,j+n,);
     while (dis[S]<nodes) ans+=dfs(S,inf);
     ans=sx-ans;
     printf("%d\n",ans);
 }