最小费用最大网络流

要打网络赛了，抓紧时间复习。
这里我先介绍一下费用流全名：最小费用最大网络流。所以我先介绍最大网络流。

最大网络流

就是给定一个图，每一条边都有一个上限，表示这条边的最大流量，求从源点流到汇点的最大流量。
对于这种题目，我们一般采用dinic算法(复杂度$n^2\ast m$).

dinic

对于dinic算法，每次增广前，我们先用 BFS 来将图分层。设源点的层数为 0 ，那么一个点的层数便是它离源点的最近距离。
分层：

1. 确定每一个点的层数，那么只要每一次都走到下一层，则在寻找增广路时保证搜索深度最小。
2. 确定是否还存在有增广路。

对于反向边，我们每一次选择一条增广路之后都要建立一条反向边，因为你所选的这条边不一定是最优的，那么你就可能要返回，所以建立反向边来解决最优问题。

dinic优化

1. 多路增广 ：每次找到一条增广路的时候，如果残余流量没有用完怎么办呢？我们可以利用残余部分流量，再找出一条增广路。这样就可以在一次 DFS 中找出多条增广路，大大提高了算法的效率。
2. 当前弧优化 ：如果一条边已经被增广过，那么它就没有可能被增广第二次。那么，我们下一次进行增广的时候，就可以不必再走那些已经被增广过的边。

费用流

在最大网络流的图里面再加上一个单位费用，求花费最小的最大网络流。
这里我只介绍类dinic算法。

类dinic算法

对于这个类dinic算法和dinic算法有啥区别呢？

1. 关于dinic算法，我们是使用BFS来做的分层，而在类dinic算法里面我们采用spfa来判断增广路，同时没有分层。
2. 在寻找增广路时，按照spfa形成的最短路去找增广路，只要能找到增广路，那么就可以放到最小费用流里面去。
   板子题

#include<stdio.h>
#include<algorithm>
#include<queue>
#include<map>
#include<set>
#include<vector>
#include<string.h>
using namespace std;
#define rep(a,b,c) for(int a = (b) ;a<=(c);a++)
typedef long long ll;
const int N = 2e5+7;
const int INF = 0x3f3f3f3f;
struct ED{
    int pre,id,w,f;
}ed[N];
int dis[N],vis[N],lin[N],ret,tot,head[N];
void add(int u,int v,int w,int f){
    ed[++tot].pre = head[u];
    ed[tot].id = v;
    ed[tot].w = w;
    ed[tot].f=f;
    head[u] = tot;
}
void init(){
    tot = 1;
    ret = 0;
}
bool spfa(int s,int t){
    memset(dis,0x3f,sizeof dis);
    memcpy(lin,head,sizeof head);
    queue<int> p;
    p.push(s);
    dis[s] = 0,vis[s] = 1;
    while(p.size()){
        int u = p.front();
        p.pop(),vis[u]=0;
        for(int i = head[u];i;i=ed[i].pre){
            int v = ed[i].id;
            if(ed[i].w>0&&dis[v]>dis[u]+ed[i].f){
                dis[v] = dis[u]+ed[i].f;
                if(!vis[v]) vis[v]=1,p.push(v);
            }
        }
    }
    return dis[t]!=INF;
}
int dfs(int u,int t,int flow){
    if(u==t||flow==0) return flow;
    vis[u] = 1;
    int ans = 0;
    for (int &i = lin[u]; i && ans < flow; i = ed[i].pre) {
        int v = ed[i].id;
        if (!vis[v] && ed[i].w && dis[v] == dis[u] + ed[i].f) {
            int x = dfs(v, t, min(ed[i].w, flow - ans));
            if (x) ret += x * ed[i].f, ed[i].w -= x, ed[i ^ 1].w += x, ans += x;
        }
    }
    vis[u] = 0;
    return ans;
}
void mcmf(int s,int t){
    int ans = 0;
    while(spfa(s,t)){
        ans += dfs(s,t,INF);
    }
    printf("%d %d\n",ans,ret);
}
int main(){
    int n,m,s,t;
    scanf("%d %d %d %d",&n,&m,&s,&t);
    init();
    rep(i,1,m){
        int u,v,w,f;
        scanf("%d %d %d %d",&u,&v,&w,&f);
        add(u,v,w,f);
        add(v,u,0,-f);
    }
    mcmf(s,t);
    //system("pause");
}