dinic求最大流//
const int INF = 1e8;
int n,m,S,T;
int e[M],ne[M],h[N],f[M],idx=0;
int q[N],d[N],cur[N]; //层数,当前弧优化
void add(int a,int b,int c){
e[idx] = b, f[idx] = c, ne[idx] = h[a] , h[a] = idx ++ ;
e[idx] = a, f[idx] = 0, ne[idx] = h[b] , h[b] = idx ++ ;
}
bool bfs(){
int hh = 0 , tt = 0;
memset(d,-1,sizeof d);
q[0] = S,d[S] = 0,cur[S] = h[S];
while(hh <= tt){
int u = q[ hh ++ ];
for(int i=h[u];~i;i=ne[i]){
int ver = e[i];
if(d[ver] == -1 && f[i]){
d[ver] = d[u] + 1; //初始化分层图
cur[ver] = h[ver]; //当前弧初始化为当前链式前向星的第一个边
if(ver == T) return true;
q[++ tt] = ver;
}
}
}
return false; //没有找到增广路径
}
int find(int u,int limit){
if(u == T) return limit;
int flow = 0; //准备从u点向后找最大的增广路上的残留容量
for(int i = cur[u]; ~i && flow < limit;i = ne[i]){ //flow < limit 必加优化
cur[u] = i;
int ver = e[i];
if(d[ver] == d[u] + 1 && f[i]){
int t = find(ver,min(f[i],limit-flow));
if(!t) d[ver] = -1; //如果此边到达不了终点,就直接删除
f[i] -= t, f[i ^ 1] += t , flow += t;
}
}
return flow;
}
int dinic(){
int ans = 0, flow = 0;
while(bfs()) while(flow = find(S,INF)) ans += flow;
return ans;
}
最小费用最大流
最小费用其实就是一个最短路 我们使用spfa求解
#include <bits/stdc++.h>
using namespace std;
const int N=5010 ,M=100010,INF=1e8
int n,m,S,T;
int h[N],ne[M],e[M],cnt,w[M],f[M];
int q[N],d[N],pre[N],incf[N];
bool st[N];
void add(int a,int b,int c,int d)
{
e[cnt]=b,f[cnt]=c,w[cnt]=d,ne[cnt]=h[a],h[a]=cnt++;
e[cnt]=a,f[cnt]=0,w[cnt]=-d,ne[cnt]=h[b],h[b]=cnt++;
}
bool spfa ()
{
int hh=0,tt=1;
memset(d,0x3f3f,sizeof d);
memset(incf,0,sizeof incf);
q[0]=S;
d[S]=0;
incf[S]=INF;
while(hh!=tt)
{
int t=q[hh++];
if(hh==N)hh=0;
st[t]=0;
for(int i=h[t];~i;i=ne[i])
{
int ver=e[i];
if(f[i]&&d[ver]>d[t]+w[i])
{
d[ver]=d[t]+w[i];
pre[ver]=i;
incf[ver]=min(f[i],incf[t]);
if(!st[ver])
{
q[tt++]=ver;
if(tt==N)tt=0;
st[ver]=true;
}
}
}
}
return incf[T]>0;
}
}
void EK(int &flow,int &cost)
{
flow=cost=0;
while(spfa())
{
int t=incf[T];
flow+=t;
cost+=d[T]*t;
for(int i=T;i!=S;i=e[pre[i]^1])
{
f[pre[i]]-=t;
f[pre[i]^1]+=t;
}
}
}
int main()
{
scanf("%d%d%d%d", &n, &m, &S, &T);
memset(h, -1, sizeof h);
while (m -- )
{
int a, b, c, d;
scanf("%d%d%d%d", &a, &b, &c, &d);
add(a, b, c, d);
}
int flow, cost;
EK(flow, cost);
printf("%d %d\n", flow, cost);
return 0;
}
变成最大流 就是把流量返回去 把权值边负
for (int i = 0; i < idx; i += 2)
{
f[i] += f[i ^ 1], f[i ^ 1] = 0;
w[i] = -w[i], w[i ^ 1] = -w[i ^ 1];
}