https://zhuanlan.zhihu.com/p/363996239
https://www.luogu.com.cn/blog/ONE-PIECE/wang-lao-liu-jiang-xie-zhi-dinic
这两篇文章讲的很好。
const int maxn=500;
const int inf1=1e9+5;
struct edge
{
int from,to,cap,flow;
};
struct dinic
{
int n, m,s,t;
vector <edge> gg;
vector <int> ggg[maxn];
bool vis[maxn];//访问标记
int d[maxn];//记录bfs的层次
int cur[maxn];//记录dfs访问的边数
void init(int l)
{
gg.clear();
for(int i=0; i<=l; i++)
{
ggg[i].clear();
}
}
void add(int from,int to,int cap)
{
gg.push_back((edge)
{
from,to,cap,0
});
gg.push_back((edge)
{
to,from,0,0
});
m=gg.size();
ggg[from].push_back(m-2);
ggg[to].push_back(m-1);
}
bool bfs()
{
memset(vis,0,sizeof(vis));
queue<int> q;
q.push(s);
d[s]=0;
vis[s]=1;//和ek不一样 dinic每次bfs一个点只访问一次
while(!q.empty())
{
int x=q.front();
q.pop();
for(int i=0; i<ggg[x].size(); i++)
{
edge &e=gg[ggg[x][i]];
if(!vis[e.to]&&e.cap>e.flow)
{
vis[e.to]=1;
d[e.to]=d[x]+1;//记录层次
q.push(e.to);
}
}
}
return vis[t];
}
int dfs(int x,int a)//当前正在处理的点 已经经过的边所能达到的最大流量
{
if(x==t||a==0)//达到终点或者无法找到流量大于0的增广路
{
return a;
}
int flow=0,f;
for(int &i=cur[x]; i<ggg[x].size(); i++)//一定注意i是引用,i增大的同时cur[x]也会增大
{ //这样当一轮dfs第二次访问x时可以直接从第一次访问
edge &e=gg[ggg[x][i]]; //x的地方接着遍历,而不用从头开始
if(d[x]+1==d[e.to]&&(f=dfs(e.to,min(a,e.cap-e.flow)))>0)//一个可行流
{
e.flow+=f;//修正正反边
gg[ggg[x][i]^1].flow-=f;
flow+=f;//修正流量
a-=f;
if(a==0)
break;
}
}
return flow;
}
int maxflow(int s,int t)
{
this->s=s;
this->t=t;
int flow=0;
while(bfs())//这是一个可行的增广路
{
memset(cur,0,sizeof(cur));//清除标记
flow+=dfs(s,inf1);
}
return flow;
}
};