网络最大流-模板（EK，dinic（+当前弧优化））
1. EK（Edmond—Karp）算法
poj1273
#include
#define swap(a,b) (a=a+b,b=a-b,a=a-b)
using namespace std;
typedef long long int LL;
const int MAXL(200);
const int INF(0x3f3f3f3f);
const int MAX(1e6);
const int mod(1e9+7);
int c[MAXL+50][MAXL+50];//残留网络的容量
int f[MAXL+50];//从源点到当前节点实际还剩多少流量可用
int pre[MAXL+50];//标记在这条路径上当前节点的前驱,同时标记该节点是否在队列中
int m,n;//n个定点，m条边
int BFS(int sce,int des)
{
queue<int>q;
memset(pre,-1,sizeof(pre));//初始化pre数组
pre[sce]=0;
f[sce]=INF;
q.push(sce);
while(!q.empty())
{
int index=q.front();
q.pop();
if(index==des)
break;
for(int i=1;i<=n;i++)
{
if(i!=sce&&c[index][i]>0&&pre[i]==-1)
{
pre[i]=index;//记录前驱
f[i]=min(c[index][i],f[index]);//找增量
q.push(i);
}
}
}
return pre[des]==-1? -1:f[des];
}
int MaxFlow(int sce,int des)
{
int increasement=0;
int sumflow=0;
while((increasement=BFS(sce,des))!=-1)
{
int k=des;//利用前驱寻找路径
while(k!=sce)
{
int last=pre[k];
c[last][k]-=increasement;//改变正向边的容量
c[k][last]+=increasement;//改变反向边的容量
k=last;
}
sumflow+=increasement;
}
return sumflow;
}
int main()
{
ios_base::sync_with_stdio(false);
while(cin>>m>>n)
{
memset(c,0,sizeof(c));
memset(f,0,sizeof(f));
for(int i=1;i<=m;i++)
{
int st,ed,power;
cin>>st>>ed>>power;
c[st][ed]+=power;//输入的时候可能会有相同边不同的流量，所以这里最好用+=
}
int ans=MaxFlow(1,n);
cout< 2.dinic算法
HDU3549
邻接表写法
#include": ";
cout< 当前弧优化版本
#include": ";
cout< 这几天做了一道HDU 3416 上面的代码卡时间了。
一直没找出原因，后来不停的debug之后发现是上面的dfs写的其实不好，递归次数可能不是最优的。
那么换一个：
int dfs(int u,int des,int increaseRoad)
{
if(u==des)
return increaseRoad;
int ret=0;
for(int k=head[u]; k!=-1; k=edge[k].next)
{
int v = edge[k].to,w=edge[k].cap;
if(level[v]==level[u]+1&&w!=0)
{
int MIN = min(increaseRoad-ret,w);
w = dfs(v,des,MIN);
if(w > 0)
{
edge[k].cap-=w;
edge[k^1].cap+=w;
ret+=w;
if(ret==increaseRoad)
return ret;
}
}
}
return ret;
}