最大流与费用流

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];
    }