hdu 4280 最大流 sap模板
给你岛的坐标求最西边到最东边的最大流
1 /* 2 最大流模板 3 sap 4 */ 5 #include<stdio.h> 6 #include<string.h> 7 #include<algorithm> 8 #include<iostream> 9 using namespace std; 10 11 const int MAXN=100010;//点数的最大值 12 const int MAXM=400010;//边数的最大值 13 const int INF=0x3f3f3f3f; 14 15 struct Node 16 { 17 int from,to,next; 18 int cap; 19 }G[MAXM]; 20 int tol; 21 int head[MAXN]; 22 int dep[MAXN]; 23 int gap[MAXN];//gap[x]=y :说明残留网络中dep[i]==x的个数为y 24 25 int n;//n是总的点的个数，包括源点和汇点 26 27 void init() 28 { 29 tol=0; 30 memset(head,-1,sizeof(head)); 31 } 32 33 void build(int u,int v,int w) 34 { 35 G[tol].from=u; 36 G[tol].to=v; 37 G[tol].cap=w; 38 G[tol].next=head[u]; 39 head[u]=tol++; 40 G[tol].from=v; 41 G[tol].to=u; 42 G[tol].cap=0; 43 G[tol].next=head[v]; 44 head[v]=tol++; 45 } 46 void BFS(int st,int ed) 47 { 48 memset(dep,-1,sizeof(dep)); 49 memset(gap,0,sizeof(gap)); 50 gap[0]=1; 51 int que[MAXN]; 52 int left,right; 53 left=right=0; 54 dep[ed]=0; 55 que[right++]=ed; 56 while(left!=right) 57 { 58 int u=que[left++]; 59 if(left==MAXN)left=0; 60 for(int i=head[u];i!=-1;i=G[i].next) 61 { 62 int v=G[i].to; 63 if(dep[v]!=-1)continue; 64 que[right++]=v; 65 if(right==MAXN)right=0; 66 dep[v]=dep[u]+1; 67 ++gap[dep[v]]; 68 } 69 } 70 } 71 int SAP(int st,int ed) 72 { 73 int res=0; 74 BFS(st,ed); 75 int cur[MAXN]; 76 int S[MAXN]; 77 int top=0; 78 memcpy(cur,head,sizeof(head)); 79 int u=st; 80 int i; 81 while(dep[st]<n){ 82 if(u==ed){ 83 int temp=INF; 84 int inser; 85 for(i=0;i<top;i++) 86 if(temp>G[S[i]].cap){ 87 temp=G[S[i]].cap; 88 inser=i; 89 } 90 for(i=0;i<top;i++){ 91 G[S[i]].cap-=temp; 92 G[S[i]^1].cap+=temp; 93 } 94 res+=temp; 95 top=inser; 96 u=G[S[top]].from; 97 } 98 if(u!=ed&&gap[dep[u]-1]==0)//出现断层，无增广路 99 break; 100 for(i=cur[u];i!=-1;i=G[i].next) 101 if(G[i].cap!=0&&dep[u]==dep[G[i].to]+1) 102 break; 103 if(i!=-1){ 104 cur[u]=i; 105 S[top++]=i; 106 u=G[i].to; 107 } 108 else{ 109 int min=n; 110 for(i=head[u];i!=-1;i=G[i].next){ 111 if(G[i].cap==0)continue; 112 if(min>dep[G[i].to]){ 113 min=dep[G[i].to]; 114 cur[u]=i; 115 } 116 } 117 --gap[dep[u]]; 118 dep[u]=min+1; 119 ++gap[dep[u]]; 120 if(u!=st)u=G[S[--top]].from; 121 } 122 } 123 return res; 124 } 125 126 int main() 127 { 128 int st,ed; 129 int m; 130 int u,v,z; 131 int T; 132 scanf("%d",&T); 133 134 while(T--){ 135 init(); 136 scanf("%d%d",&n,&m); 137 int minx=10000000; 138 int maxx=-10000000; 139 int x,y; 140 for(int i=1;i<=n;i++){ 141 scanf("%d%d",&x,&y); 142 if(minx>x){ 143 minx=x; 144 st=i; 145 } 146 if(maxx<x){ 147 maxx=x; 148 ed=i; 149 } 150 } 151 while(m--){ 152 scanf("%d%d%d",&u,&v,&z); 153 build(u,v,z); 154 build(v,u,z); 155 } 156 //n一定是点的总数，这是使用SAP模板需要注意的 157 int ans=SAP(st,ed); 158 printf("%d\n",ans); 159 } 160 return 0; 161 }
