HDU 3549 Flow Problem（最大流裸题，EK解法）

EK现在用起来已经比较熟练了，比较简单的题目，直接就可以解决，但发现EK，其实效率有些低，是O（V*E^2）的，搞不好，还会超时。同时看别人的博客，发现用EK的比较少。今天，听学长讲，dinic用的比较多，sap学起来有一定难度，接下来抓紧学一下dinic！

好了，废话不多说了，进入主题！

原题链接：http://acm.hdu.edu.cn/showproblem.php?pid=3549

题意：在给定的有向图中，求点1到n的最大流。

注意点：EK两点之间只能有一条正向边，而顶点最多只有15个，边却最多有1000条，15*14/2<<1000，显然两点间会有重边。所以处理时，既不是取两点间边的最大值，也不是让新加入的边覆盖，而是不断累加。

下面是比较详细的注解的代码:

#include <iostream>
#include <cstdio>
#include <queue>
#include <string.h> 
using namespace std;
#define arraySize 20
#define inf 1000000000
int capacity[arraySize][arraySize],flow[arraySize][arraySize],max_flow[arraySize],pre[arraySize];
//capacity存储点之间最大流量，flow存储点之间当前已经流过的流量 
//max_flow存储每次遍历过程中的值,pre记录查找过程中每个节点的前一节点，用于后续更新

int Edmonds_Karp(int source,int target)//源点,汇点 
{
	//初始化 
	queue <int> store;
	int ans=0,cur;
	//cur当前节点 
	memset(flow,0,sizeof(flow));
	while(true)//一直寻找增广路 
	{
	   memset(max_flow,0,sizeof(max_flow));	
	   memset(pre,0,sizeof(pre));
	   store.push(source);
	   max_flow[source]=inf;
	   while(!store.empty())
	   {
   	     cur=store.front();
   	     store.pop();
   	     for(int next=source;next<=target;next++)
   	     {
   	     	//max_flow[next]恰可以用于标记是否访问过，同时要保证两点之间还有剩余流量 
                //这个过程中，可能会出现多条可行路径，但因为汇点只有一个会被先到达的路径抢占，故每个过程只能找到一条
     	    if(!max_flow[next]&&capacity[cur][next]>flow[cur][next])
            {
                 store.push(next);  
                 //如果这两个点之间的值，比之前的最小值还小，则更新 
                 max_flow[next]=min(max_flow[cur],capacity[cur][next]-flow[cur][next]);
                 //记录前一个节点，用于后续更新 
                 pre[next]=cur;
 			} 	
	     }
   	   }
   	   //说明已经找不到增广路了 
   	   if(max_flow[target]==0)break;
	   //更新操作 
   	   for(int u=target;u!=source;u=pre[u])  
        {  
            flow[pre[u]][u]+=max_flow[target];
			//反向边  
            flow[u][pre[u]]-=max_flow[target];  
        }  
        ans+=max_flow[target]; 
	}
	return ans;
}
int main()
{
	int t,n,m,fm,to,val,cnt=0;
	scanf("%d",&t);
	while(t--)
	{
	  memset(capacity,0,sizeof(capacity));
	  scanf("%d%d",&n,&m);
	  for(int i=1;i<=m;i++)
	  {
  		 scanf("%d%d%d",&fm,&to,&val);
  		 capacity[fm][to]+=val;
  	  }	
  	  printf("Case %d: %d\n",++cnt,Edmonds_Karp(1,n));
	}
	return 0;
}