HDOJ 4280 - Island Transport 裸最大流(Dinic) ..申请更多栈空间...

               题意:

                        有很多小岛在XY坐标系上...有些岛之间是有路径的..并且每条路劲有容量上限..问做最西(左)侧的点到最东(点)侧的点..最多能运载多少旅客过去..

               题解:

                       这道题....我去年做网络赛的时候..写dinic..爆栈...然后就没有然后了...然后队友写sap超时....给跪..还好其他题做的不错.进了现场赛...

                       爆栈是HDOJ喜闻乐见的情况..现在知道在前面加语句了....我估计我以前那个dinic就算不爆栈也会超时..现在这个效率不错..不过要是数据抢了也得跪..比较和谐的方法是化成对偶图来做...还要学习...

Program:

#pragma comment(linker,"/STACK:102400000,102400000") 
#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<queue>
#define MAXN 200005
#define MAXM 200005
#define oo 1000000007
#define ll long long
using namespace std;
int X[MAXN];
struct Dinic
{
      struct node 
      {
             int x,y,c,next;
      }line[MAXM];   
      int Lnum,_next[MAXN],dis[MAXN];
      void initial(int n) 
      {
             for (int i=0;i<=n;i++) _next[i]=-1;
             Lnum=-1;
      } 
      void addline(int x,int y,int c)
      {
             line[++Lnum].next=_next[x],_next[x]=Lnum;
             line[Lnum].x=x,line[Lnum].y=y,line[Lnum].c=c;
             line[++Lnum].next=_next[y],_next[y]=Lnum;
             line[Lnum].x=y,line[Lnum].y=x,line[Lnum].c=c;             
      }
      bool BFS(int s,int e)
      { 
             queue<int> Q;
             while (!Q.empty()) Q.pop();
             memset(dis,0,sizeof(dis));
             dis[s]=1;
             Q.push(s);
             while (!Q.empty())
             {
                   int h,k;
                   h=Q.front(),Q.pop();
                   if (h==e) return dis[e];
                   for (k=_next[h];k!=-1;k=line[k].next)
                      if (line[k].c && !dis[line[k].y])
                         dis[line[k].y]=dis[h]+1,Q.push(line[k].y);                 
             } 
             return false;
      }
      int dfs(int x,int flow,int e)  
      {     
             if (x==e) return flow;   
             int temp,cost=0;  
             for (int k=_next[x];k!=-1;k=line[k].next)  
             if (line[k].c && dis[line[k].y]==dis[x]+1)  
             {  
                    temp=dfs(line[k].y,min(flow-cost,line[k].c),e);   
                    if (temp)  
                    {  
                           line[k].c-=temp,line[k^1].c+=temp;  
                           cost+=temp;  
                           if (flow==cost) return cost;  
                    }else dis[line[k].y]=-1;  
             }  
             return cost;  
      }  
      int MaxFlow(int s,int e)
      {
             int MaxFlow=0;
             while (BFS(s,e)) 
                MaxFlow+=dfs(s,oo,e); 
             return MaxFlow;
      }
}T;
int main()  
{          
       int n,m,s,e,i,x,y,c,cases; 
       scanf("%d",&cases);
       while (cases--)  
       {
               scanf("%d%d",&n,&m);
               s=e=1;
               for (i=1;i<=n;i++)
               {
                      scanf("%d%d",&x,&y);
                      X[i]=x;
                      if (x<X[s]) s=i;
                      if (x>X[e]) e=i;
               }
               T.initial(n);
               while (m--)
               {
                       scanf("%d%d%d",&x,&y,&c);
                       T.addline(x,y,c); 
               }
               printf("%d\n",T.MaxFlow(s,e));
       }
       return 0;  
}