POJ 2455 二分+最大流

题意：

问从起点到终点有T条不相交路径的最大边的最小值。

题解：

典型的二分答案！最大流还真爱和二分一起出啊~

 

View Code 
  1 #include <iostream>
  2 #include <cstring>
  3 #include <cstdlib>
  4 #include <cstdio>
  5 #include <algorithm>
  6 
  7 #define N 1000
  8 #define M 1000000
  9 #define INF 1e9
 10 
 11 using namespace std;
 12 
 13 int head[N],next[M],to[M],len[M];
 14 int layer[N],q[M*10];
 15 int n,m,cnt,tot,S,T,l,r,mid,sm;
 16 
 17 struct PX
 18 {
 19     int a,c,b;
 20 }px[M];
 21 
 22 inline bool cmp(const PX &a,const PX &b)
 23 {
 24     return a.c<b.c;
 25 }
 26 
 27 inline void add(int u,int v,int w)
 28 {
 29     to[cnt]=v; len[cnt]=w; next[cnt]=head[u]; head[u]=cnt++;
 30     to[cnt]=u; len[cnt]=0; next[cnt]=head[v]; head[v]=cnt++;
 31 }
 32 
 33 inline void read()
 34 {
 35     l=0;r=0; S=1; T=n; tot=0;
 36     for(int i=1;i<=m;i++)
 37     {
 38         tot++;
 39         scanf("%d%d%d",&px[tot].a,&px[tot].b,&px[tot].c);
 40         ++tot; px[tot].a=px[tot-1].b; px[tot].b=px[tot-1].a; px[tot].c=px[tot-1].c;
 41         r=max(r,px[tot].c);
 42     }
 43     sort(px+1,px+1+tot,cmp);
 44 }
 45 
 46 inline void build()
 47 {
 48     memset(head,-1,sizeof head); cnt=0;
 49     for(int i=1;i<=tot;i++)
 50     {
 51         if(px[i].c>mid) break;
 52         add(px[i].a,px[i].b,1);
 53     }
 54 }
 55 
 56 inline bool bfs()
 57 {
 58     memset(layer,-1,sizeof layer);
 59     int h=1,t=2,sta;
 60     q[1]=S; layer[S]=0;
 61     while(h<t)
 62     {
 63         sta=q[h++];
 64         for(int i=head[sta];~i;i=next[i])
 65             if(len[i]&&layer[to[i]]<0)
 66             {
 67                 layer[to[i]]=layer[sta]+1;
 68                 q[t++]=to[i];
 69             }
 70     }
 71     return layer[T]!=-1;
 72 }
 73 
 74 inline int find(int u,int cur_flow)
 75 {
 76     if(u==T) return cur_flow;
 77     int res=0,tmp;
 78     for(int i=head[u];~i&&res<cur_flow;i=next[i])
 79         if(len[i]&&layer[to[i]]==layer[u]+1)
 80         {
 81             tmp=find(to[i],min(cur_flow-res,len[i]));
 82             len[i]-=tmp; len[i^1]+=tmp; res+=tmp;
 83         }
 84     if(!res) layer[u]=-1;
 85     return res;
 86 }
 87 
 88 inline int dinic()
 89 {
 90     int ans=0;
 91     while(bfs()) ans+=find(S,INF);
 92     return ans;
 93 }
 94 
 95 inline void go()
 96 {
 97     int ans;
 98     while(l<=r)
 99     {
100         mid=(l+r)>>1;
101         build();
102         if(dinic()>=sm) ans=mid,r=mid-1;
103         else l=mid+1;
104     }
105     printf("%d\n",ans);
106 }
107 
108 int main()
109 {
110     while(scanf("%d%d%d",&n,&m,&sm)!=EOF) read(),go();
111     return 0;
112 }

 

 

posted @ 2013-01-07 23:45  proverbs  阅读(141)  评论(0)  编辑  收藏  举报