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 }
  • 相关阅读:
    IDEA在debug时修改变量值
    CSS覆盖公共样式中的某个属性
    POI获取单元格的宽和高
    MySQL将一张表的某些列数据,复制到另外一张表,并且修改某些内容
    哈希是什么?为什么哈希存取比较快?
    工厂模式
    观察者模式
    instanceof判断的对象可以是接口
    JeeSite框架中httpSession.invalidate();无效
    HSQL可视化工具
  • 原文地址:https://www.cnblogs.com/WABoss/p/5281056.html
Copyright © 2011-2022 走看看