poj 3228(二分+最大流)

题目链接：http://poj.org/problem?id=3228

思路：增设一个超级源点和一个超级汇点，源点与每一个gold相连，容量为gold数量，汇点与仓库相连，容量为仓库的容量，然后就是二分最小的最大相邻距离，跑最大流验证即可。最大流用的是别人的Dinic模版。

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<queue>
 using namespace std;
 #define MAXN 222
 #define inf 1<<30

 int vs,vt;
 int dep[MAXN];
 int map[MAXN][MAXN];

 int gold[MAXN],store[MAXN];
 int dist[MAXN][MAXN];
 void Build_Map(int limit,int n)
 {
     memset(map,,sizeof(map));
     for(int i=;i<=n;i++)
         for(int j=;j<=n;j++)
             if(dist[i][j]<=limit)map[i][j]=inf;
     for(int i=;i<=n;i++){
         map[vs][i]=gold[i];
         map[i][vt]=store[i];
     }
 }

 int BFS(){
     queue<int> q;
     while(!q.empty())
         q.pop();
     memset(dep,-,sizeof(dep));
     dep[vs]=;
     q.push(vs);
     while(!q.empty()){
         int u=q.front();
         q.pop();
         for(int v=vs;v<=vt;v++)
             if(map[u][v]> && dep[v]==-){
                 dep[v]=dep[u]+;
                 q.push(v);
             }
     }
     return dep[vt]!=-;
 }

 int DFS(int u,int minx){
     if(u==vt)
         return minx;
     int tmp;
     for(int v=vs;v<=vt;v++)
         if(map[u][v]> && dep[v]==dep[u]+ && (tmp=DFS(v,min(minx,map[u][v])))){
             map[u][v]-=tmp;
             map[v][u]+=tmp;
             return tmp;
         }
     dep[u]=-;
     return ;
 }

 int Dinic(){
     int ans=,tmp;
     while(BFS()){
         while(){
             tmp=DFS(vs,inf);
             if(tmp==)
                 break;
             ans+=tmp;
         }
     }
     return ans;
 }

 int main()
 {
     int total,n,m,a,b,c;
     while(~scanf("%d",&n)){
         if(n==)break;
         total=,vs=,vt=n+;
         for(int i=;i<=n;i++)
             for(int j=;j<=n;j++)
                 dist[i][j]=inf;
         for(int i=;i<=n;i++){
             scanf("%d",&gold[i]);
             total+=gold[i];
         }
         for(int i=;i<=n;i++)scanf("%d",&store[i]);
         scanf("%d",&m);
         while(m--){
             scanf("%d%d%d",&a,&b,&c);
             dist[a][b]=dist[b][a]=c;
         }

         int low=,high=,mid,ans=-;
         while(low<=high){
             mid=(low+high)>>;
             Build_Map(mid,n);
             if(Dinic()==total){
                 ans=mid;
                 high=mid-;
             }else
                 low=mid+;
         }
         if(ans==-){
             printf("No Solution\n");
         }else
             printf("%d\n",ans);
     }
     return ;
 }