hdu 4280 最大流 sap模板

给你岛的坐标求最西边到最东边的最大流

 /*
 最大流模板
 sap
 */
 #include<stdio.h>
 #include<string.h>
 #include<algorithm>
 #include<iostream>
 using namespace std;
 
 const int MAXN=;//点数的最大值
 const int MAXM=;//边数的最大值
 const int INF=0x3f3f3f3f;
 
 struct Node
 {
     int from,to,next;
     int cap;
 }G[MAXM];
 int tol;
 int head[MAXN];
 int dep[MAXN];
 int gap[MAXN];//gap[x]=y :说明残留网络中dep[i]==x的个数为y
 
 int n;//n是总的点的个数，包括源点和汇点
 
 void init()
 {
     tol=;
     memset(head,-,sizeof(head));
 }
 
 void build(int u,int v,int w)
 {
     G[tol].from=u;
     G[tol].to=v;
     G[tol].cap=w;
     G[tol].next=head[u];
     head[u]=tol++;
     G[tol].from=v;
     G[tol].to=u;
     G[tol].cap=;
     G[tol].next=head[v];
     head[v]=tol++;
 }
 void BFS(int st,int ed)
 {
     memset(dep,-,sizeof(dep));
     memset(gap,,sizeof(gap));
     gap[]=;
     int que[MAXN];
     int left,right;
     left=right=;
     dep[ed]=;
     que[right++]=ed;
     while(left!=right)
     {
         int u=que[left++];
         if(left==MAXN)left=;
         for(int i=head[u];i!=-;i=G[i].next)
         {
             int v=G[i].to;
             if(dep[v]!=-)continue;
             que[right++]=v;
             if(right==MAXN)right=;
             dep[v]=dep[u]+;
             ++gap[dep[v]];
         }
     }
 }
 int SAP(int st,int ed)
 {
     int res=;
     BFS(st,ed);
     int cur[MAXN];
     int S[MAXN];
     int top=;
     memcpy(cur,head,sizeof(head));
     int u=st;
     int i;
     while(dep[st]<n){
         if(u==ed){
             int temp=INF;
             int inser;
             for(i=;i<top;i++)
                if(temp>G[S[i]].cap){
                    temp=G[S[i]].cap;
                    inser=i;
                }
             for(i=;i<top;i++){
                 G[S[i]].cap-=temp;
                 G[S[i]^].cap+=temp;
             }
             res+=temp;
             top=inser;
             u=G[S[top]].from;
         }
         if(u!=ed&&gap[dep[u]-]==)//出现断层，无增广路
           break;
         for(i=cur[u];i!=-;i=G[i].next)
             if(G[i].cap!=&&dep[u]==dep[G[i].to]+)
              break;
         if(i!=-){
             cur[u]=i;
             S[top++]=i;
             u=G[i].to;
         }
         else{
             int min=n;
             for(i=head[u];i!=-;i=G[i].next){
                 if(G[i].cap==)continue;
                 if(min>dep[G[i].to]){
                     min=dep[G[i].to];
                     cur[u]=i;
                 }
             }
             --gap[dep[u]];
             dep[u]=min+;
             ++gap[dep[u]];
             if(u!=st)u=G[S[--top]].from;
         }
     }
     return res;
 }
 
 int main()
 {
     int st,ed;
     int m;
     int u,v,z;
     int T;
     scanf("%d",&T);
 
     while(T--){
         init();
         scanf("%d%d",&n,&m);
         int minx=;
         int maxx=-;
         int x,y;
         for(int i=;i<=n;i++){
             scanf("%d%d",&x,&y);
             if(minx>x){
                 minx=x;
                 st=i;
             }
             if(maxx<x){
                 maxx=x;
                 ed=i;
             }
         }
         while(m--){
             scanf("%d%d%d",&u,&v,&z);
             build(u,v,z);
             build(v,u,z);
         }
         //n一定是点的总数，这是使用SAP模板需要注意的
         int ans=SAP(st,ed);
         printf("%d\n",ans);
     }
     return ;
 }