zoukankan      html  css  js  c++  java
  • ZOJ2332 Gems(最大流)

    题目大概说,alsomagic有宝石,宝石有颜色和形状两个属性;他有一种法力可以将某些颜色形状的宝石转化成另一种颜色形状的宝石;他的女朋友对各个颜色都有一定的容忍数量,而他自己也对各个形状都有一定的容忍数量,超过这个数量就不行;问能否有成功把宝石分成两份,一份给他,一份给他女朋友。

    题目这句话“But each way of transform can only be used ONCE!”理解不太清晰。。反正容量网络是这样建的:

    • 源点、各种宝石、各种颜色、各种形状、汇点
    • 源点向每种宝石连容量为这种宝石数量的边
    • 可以互相转化的两种宝石互相间连一条容量INF的边
    • 各种宝石各自分别向其对应的颜色、形状连容量INF的边
    • 各个颜色向汇点连容量为各自的容忍数量的边,各个形状同样如此
    • 跑最大流,如果结果等于所有宝石数量,就是Yes,否则No。
      1 #include<cstdio>
      2 #include<cstring>
      3 #include<queue>
      4 #include<algorithm>
      5 using namespace std;
      6 #define INF (1<<30)
      7 #define MAXN 222
      8 #define MAXM 222*222*2
      9 
     10 struct Edge{
     11     int v,cap,flow,next;
     12 }edge[MAXM];
     13 int vs,vt,NE,NV;
     14 int head[MAXN];
     15 
     16 void addEdge(int u,int v,int cap){
     17     edge[NE].v=v; edge[NE].cap=cap; edge[NE].flow=0;
     18     edge[NE].next=head[u]; head[u]=NE++;
     19     edge[NE].v=u; edge[NE].cap=0; edge[NE].flow=0;
     20     edge[NE].next=head[v]; head[v]=NE++;
     21 }
     22 
     23 int level[MAXN];
     24 int gap[MAXN];
     25 void bfs(){
     26     memset(level,-1,sizeof(level));
     27     memset(gap,0,sizeof(gap));
     28     level[vt]=0;
     29     gap[level[vt]]++;
     30     queue<int> que;
     31     que.push(vt);
     32     while(!que.empty()){
     33         int u=que.front(); que.pop();
     34         for(int i=head[u]; i!=-1; i=edge[i].next){
     35             int v=edge[i].v;
     36             if(level[v]!=-1) continue;
     37             level[v]=level[u]+1;
     38             gap[level[v]]++;
     39             que.push(v);
     40         }
     41     }
     42 }
     43 
     44 int pre[MAXN];
     45 int cur[MAXN];
     46 int ISAP(){
     47     bfs();
     48     memset(pre,-1,sizeof(pre));
     49     memcpy(cur,head,sizeof(head));
     50     int u=pre[vs]=vs,flow=0,aug=INF;
     51     gap[0]=NV;
     52     while(level[vs]<NV){
     53         bool flag=false;
     54         for(int &i=cur[u]; i!=-1; i=edge[i].next){
     55             int v=edge[i].v;
     56             if(edge[i].cap!=edge[i].flow && level[u]==level[v]+1){
     57                 flag=true;
     58                 pre[v]=u;
     59                 u=v;
     60                 //aug=(aug==-1?edge[i].cap:min(aug,edge[i].cap));
     61                 aug=min(aug,edge[i].cap-edge[i].flow);
     62                 if(v==vt){
     63                     flow+=aug;
     64                     for(u=pre[v]; v!=vs; v=u,u=pre[u]){
     65                         edge[cur[u]].flow+=aug;
     66                         edge[cur[u]^1].flow-=aug;
     67                     }
     68                     //aug=-1;
     69                     aug=INF;
     70                 }
     71                 break;
     72             }
     73         }
     74         if(flag) continue;
     75         int minlevel=NV;
     76         for(int i=head[u]; i!=-1; i=edge[i].next){
     77             int v=edge[i].v;
     78             if(edge[i].cap!=edge[i].flow && level[v]<minlevel){
     79                 minlevel=level[v];
     80                 cur[u]=i;
     81             }
     82         }
     83         if(--gap[level[u]]==0) break;
     84         level[u]=minlevel+1;
     85         gap[level[u]]++;
     86         u=pre[u];
     87     }
     88     return flow;
     89 }
     90 
     91 int main(){
     92     int t,n,m,a,r1,c1,r2,c2;
     93     scanf("%d",&t);
     94     while(t--){
     95         scanf("%d%d",&n,&m);
     96         vs=n*m+n+m; vt=vs+1; NV=vt+1; NE=0;
     97         memset(head,-1,sizeof(head));
     98         int tot=0;
     99         for(int i=0; i<n; ++i){
    100             for(int j=0; j<m; ++j){
    101                 scanf("%d",&a); tot+=a;
    102                 addEdge(vs,i*m+j,a);
    103                 addEdge(i*m+j,i+n*m,INF);
    104                 addEdge(i*m+j,j+n*m+n,INF);
    105             }
    106         }
    107         scanf("%d",&a);
    108         while(a--){
    109             scanf("%d%d%d%d",&r1,&c1,&r2,&c2);
    110             addEdge(r1*m+c1,r2*m+c2,INF);
    111             addEdge(r2*m+c2,r1*m+c1,INF);
    112         }
    113         for(int i=0; i<n; ++i){
    114             scanf("%d",&a);
    115             addEdge(i+n*m,vt,a);
    116         }
    117         for(int i=0; i<m; ++i){
    118             scanf("%d",&a);
    119             addEdge(i+n*m+n,vt,a);
    120         }
    121         if(ISAP()==tot) puts("Yes");
    122         else puts("No");
    123     } 
    124     return 0;
    125 }
  • 相关阅读:
    异常
    C++中的mutable,volatile,explicit关键字
    Vi配置文件--Vimrc
    结构体和类的区别
    [转]恢复视力的方法(500度以下)
    与struct相关的宏定义 ---今Tencent笔试用到的
    如何在C++中调用C的代码
    C中如何调用C++函数?
    技术博走起
    Shell常见命令实践
  • 原文地址:https://www.cnblogs.com/WABoss/p/5281056.html
Copyright © 2011-2022 走看看