POJ 3228 - Gold Transportation 构图最大流..继续用超强dinic模板..

                   题意:

                            有一些城市....有些城市之间有路径..某些城市可能有金子..某些城市可能有仓库.现在想让所有的金子都运到仓库中并且所经过所有边的最大值最小...

                   题解:

                            最大值最小..果断是二分..而二分枚举后的判断.就用最大流...超级源点向所有有金子的城市做边.容量为其金子数量....然后添加所有边长度不大于当前枚举值的边(无向的).容量为无穷大..再所有有仓库的往超级汇点做边..容量为其仓库容量...跑最大流..若最大流等于金子总数量..则说明在该长度限制下是可以的..

Program:

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<queue>
#define MAXN 305
#define MAXM 80005
#define oo 1000000007
#define ll long long
using namespace std;  
struct Dinic    
{    
       struct node  
       {   
             int c,u,v,next;  
       }edge[MAXM];  
       int ne,head[MAXN];  
       int cur[MAXN], ps[MAXN], dep[MAXN];  
       void initial()  
       {  
             ne=2;  
             memset(head,0,sizeof(head));   
       }  
       void addedge(int u, int v,int c)  
       {   
             edge[ne].u=u,edge[ne].v=v,edge[ne].c=c,edge[ne].next=head[u];  
             head[u]=ne++;  
             edge[ne].u=v,edge[ne].v=u,edge[ne].c=c,edge[ne].next=head[v];  
             head[v]=ne++;  
       }  
       int MaxFlow(int s,int t)  
       {                       
             int tr, res = 0;  
             int i,j,k,f,r,top;  
             while(1)  
             {  
                    memset(dep, -1, sizeof(dep));  
                    for(f=dep[ps[0]=s]=0,r=1;f!= r;)  
                       for(i=ps[f++],j=head[i];j;j=edge[j].next)  
                         if(edge[j].c&&dep[k=edge[j].v]==-1)  
                         {  
                               dep[k]=dep[i]+1;  
                               ps[r++]=k;  
                               if(k == t){  f=r; break;  }  
                         }  
                    if(dep[t]==-1) break;  
                    memcpy(cur,head,sizeof(cur));  
                    i=s,top=0;  
                    while(1)  
                    {  
                         if(i==t)  
                         {  
                               for(tr=oo,k=0;k<top;k++)  
                                  if(edge[ps[k]].c<tr)  
                                     tr=edge[ps[f=k]].c;  
                               for(k=0;k<top;k++)  
                               {  
                                     edge[ps[k]].c-=tr;  
                                     edge[ps[k]^1].c+=tr;  
                               }  
                               i=edge[ps[top=f]].u;  
                               res+= tr;  
                         }  
                         for(j=cur[i];cur[i];j=cur[i]=edge[cur[i]].next)   
                             if(edge[j].c && dep[i]+1==dep[edge[j].v]) break;   
                         if(cur[i])  ps[top++]=cur[i],i=edge[cur[i]].v;   
                         else  
                         {  
                                 if(!top) break;  
                                 dep[i]=-1;  
                                 i=edge[ps[--top]].u;  
                         }  
                   }  
             }  
             return res;  
      }  
}T;
struct node
{
      int x,y,d;
}L[MAXM];
int S[MAXN],E[MAXN];
bool judge(int R,int n,int m)
{
      int s=n+1,e=s+1,sum=0,i;
      T.initial();
      for (i=1;i<=n;i++)
      {
           if (S[i]) T.addedge(s,i,S[i]),sum+=S[i];
           if (E[i]) T.addedge(i,e,E[i]);
      }
      for (i=1;i<=m;i++)
         if (L[i].d<=R) T.addedge(L[i].x,L[i].y,oo);      
      return T.MaxFlow(s,e)==sum;
}
int main() 
{       
      int n,m,i;
      while (~scanf("%d",&n) && n)
      {
              for (i=1;i<=n;i++) scanf("%d",&S[i]);
              for (i=1;i<=n;i++) scanf("%d",&E[i]);
              scanf("%d",&m);
              for (i=1;i<=m;i++) scanf("%d%d%d",&L[i].x,&L[i].y,&L[i].d);
              int l=-1,r=10005,mid;
              while (r-l>1)
              {
                      mid=l+r>>1;
                      if (judge(mid,n,m)) r=mid;
                                    else  l=mid;
              }
              if (r==10005) printf("No Solution\n");
                      else  printf("%d\n",r);
      }
      return 0;
}