zoukankan      html  css  js  c++  java
  • SPOJ287 Smart Network Administrator(最大流)

    题目大概是说,一个村庄有n间房子,房子间有m条双向路相连。1号房子有网络,有k间房子要通过与1号房子相连联网,且一条路上不能有同样颜色的线缆,问最少要用几种颜色的线缆。

    二分枚举颜色个数,建立容量网络跑最大流判断解是否成立:源点是1,所有要联网的房子向汇点连容量1的边,所有双向边改为容量为枚举的解的边,如果最大流等于k那这个解就成立。

    注意重边。。

      1 #include<cstdio>
      2 #include<cstring>
      3 #include<queue>
      4 #include<algorithm>
      5 using namespace std;
      6 #define INF (1<<30)
      7 #define MAXN 555
      8 #define MAXM 555*555*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 n,m,k,h[555];
     92 bool G[555][555];
     93 bool isok(int mx){
     94     vs=1; vt=0; NV=n; NE=0;
     95     memset(head,-1,sizeof(head));
     96     for(int i=0; i<k; ++i) addEdge(h[i],vt,1);
     97     for(int i=1; i<=n; ++i){
     98         for(int j=1; j<=n; ++j){
     99             if(G[i][j]) addEdge(i,j,mx);
    100         } 
    101     }
    102     return ISAP()==k;
    103 }
    104 int main(){
    105     int t,a,b;
    106     scanf("%d",&t);
    107     while(t--){
    108         memset(G,0,sizeof(G));
    109         scanf("%d%d%d",&n,&m,&k);
    110         for(int i=0; i<k; ++i) scanf("%d",h+i);
    111         for(int i=0; i<m; ++i){
    112             scanf("%d%d",&a,&b);
    113             G[a][b]=G[b][a]=1;
    114         }
    115         int l=0,r=k;
    116         while(l<r){
    117             int mid=l+r>>1;
    118             if(isok(mid)) r=mid;
    119             else l=mid+1;
    120         }
    121         printf("%d
    ",l);
    122     }
    123     return 0;
    124 } 
  • 相关阅读:
    Python3.x基础学习-静态方法和类方法
    Python3.x基础学习-@property
    Python3.x基础学习-装饰器
    SpringMVC学习(七)
    SpringMVC学习(六)
    SpringMVC学习(五)
    整合SSM(二)
    整合SSM(一)
    SpringMVC学习(四)
    SpringMVC学习(三)
  • 原文地址:https://www.cnblogs.com/WABoss/p/5285038.html
Copyright © 2011-2022 走看看