Gems

zoj2332:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=2332

题意：这一道题的题意，我看了很久，也没有看明白，最终还是理解错了。题目的意思是有n*m种石头，有n种形状，m种颜色。一个人对于每一种形状有一定的容忍度，另外一个人对于每一种颜色有一定容忍度，然后问你能把这些石头分给这两个人，使得每个人能够满意。

题解：每一种石头作为一个点，与源点连边，边的容量就是石头的数量，然后每一行作为一个点，对应每一种形状的，容量是INF，然后是每一列一个点，与宝石连接，INF，然后每个行与BF连接，就是每一种形状的容忍度，每一列与GF连接，就是每一种颜色的容忍度，然后BF，GF与汇点连接，INF。

 #include<iostream>
 #include<cstring>
 #include<algorithm>
 #include<cstdio>
 #include<queue>
 #define INF 100000000
 using namespace std;
 const int N=;
 const int M=;
 struct Node{
    int v;
    int f;
    int next;
 }edge[M];
 int n,m,u,v,cnt,sx,ex;
 int head[N],pre[N];
 int mp[N][N],a1[N],a2[N];
 void init(){
     cnt=;
     memset(head,-,sizeof(head));
 }
 void add(int u,int v,int w){
     edge[cnt].v=v;
     edge[cnt].f=w;
     edge[cnt].next=head[u];
     head[u]=cnt++;
     edge[cnt].f=;
     edge[cnt].v=u;
     edge[cnt].next=head[v];
     head[v]=cnt++;
 }
 bool BFS(){
   memset(pre,,sizeof(pre));
   pre[sx]=;
   queue<int>Q;
   Q.push(sx);
  while(!Q.empty()){
      int d=Q.front();
      Q.pop();
      for(int i=head[d];i!=-;i=edge[i].next    ){
         if(edge[i].f&&!pre[edge[i].v]){
             pre[edge[i].v]=pre[d]+;
             Q.push(edge[i].v);
         }
      }
   }
  return pre[ex]>;
 }
 int dinic(int flow,int ps){
     int f=flow;
      if(ps==ex)return f;
      for(int i=head[ps];i!=-;i=edge[i].next){
         if(edge[i].f&&pre[edge[i].v]==pre[ps]+){
             int a=edge[i].f;
             int t=dinic(min(a,flow),edge[i].v);
               edge[i].f-=t;
               edge[i^].f+=t;
               flow-=t;
              if(flow<=)break;
         }

      }
       if(f-flow<=)pre[ps]=-;
       return f-flow;
 }
 int solve(){
     int sum=;
     while(BFS())
         sum+=dinic(INF,sx);
     return sum;
 }
 int main() {
     int T,k,temp,sum,d,t1,t2,t3,t4,tt=;
     scanf("%d",&T);
     while(T--) {
          scanf("%d%d",&n,&m);
          init();sum=;
         for(int i=;i<=n;i++)
         for(int j=;j<=m;j++){
             scanf("%d",&temp);
             sum+=temp;
             add(,(i-)*m+j,temp);
         }

         scanf("%d",&d);
         for(int i=;i<=d;i++){
             scanf("%d%d%d%d",&t1,&t2,&t3,&t4);
             t2++;t4++;
             add(t1*m+t2,t3*m+t4,INF);
             add(t3*m+t4,t1*m+t2,INF);
         }
         for(int i=;i<=n;i++){
             scanf("%d",&temp);
             for(int j=;j<=m;j++){
                 add((i-)*m+j,n*m+i,INF);
             }
             add(n*m+i,n*m+m+n+,temp);
         }
         for(int i=;i<=m;i++){
             scanf("%d",&temp);
             for(int j=;j<=n;j++){
             add((j-)*m+i,n*m+n+i,INF);
             }
             add(n*m+n+i,n*m+m+n+,temp);
         }
         add(n*m+n+m+,n*m+n+m+,INF);
         add(n*m+n+m+,n*m+n+m+,INF);
         sx=,ex=n*m+n+m+;
         int td=solve();
         if(td==sum)printf("Yes\n");
         else
             printf("No\n");
     }
     return ;
 }