【模板】最小费用最大流

传送门：洛谷-【模板】最小费用最大流

---

题意

给出一个网络图，以及其源点和汇点，每条边已知其最大流量和单位流量费用，求出其网络最大流和在最大流情况下的最小费用。

---

数据范围

对于30%的数据：N<=10，M<=10

对于70%的数据：N<=1000，M<=1000

对于100%的数据：N<=5000，M<=50000

---

题解

没有讲解，只有代码

#include<bits/stdc++.h>
using namespace std;
const int inf=0x7fffffff;
const int N=5050,M=4e5+100;
int n,m,S,T,flow,cost,vis[N],dis[N];
int head[N],to[M],nxt[M],w[M],ct[M],tot=1;
char c;

deque<int>Q;

template<class TT>
inline void rd(TT &x)
{
    c=getchar();x=0;int f=0;
    for(;!isdigit(c);c=getchar()) if(c=='-') f=1;
    for(;isdigit(c);c=getchar()) x=x*10+(c^48);
    if(f) x=-x;
}

inline void lk(int u,int v,int vv,int f)
{to[++tot]=v;nxt[tot]=head[u];head[u]=tot;w[tot]=vv;ct[tot]=f;}

inline bool spfa()
{
    int i,j,x;
    for(i=1;i<=n;++i) vis[i]=0,dis[i]=inf;
    vis[T]=1;dis[T]=0;Q.push_back(T);
    while(!Q.empty()){
        x=Q.front();Q.pop_front();
        for(i=head[x];i;i=nxt[i]){
            j=to[i];if((!w[i^1])||(dis[j]<=dis[x]-ct[i])) continue;
            dis[j]=dis[x]-ct[i];
            if(!vis[j]){
               if(Q.empty() || dis[Q.front()]<=dis[j]) Q.push_back(j);
               else Q.push_front(j);vis[j]=1;
            }
        }
        vis[x]=0;
    }
    return dis[S]<inf;
}

inline int dfs(int x,int f)
{
    vis[x]=1;
    if(x==T) return f;
    int i,j,res,ss=0;
    for(i=head[x];i;i=nxt[i]){
        j=to[i];if(vis[j] || !w[i] || dis[j]!=dis[x]-ct[i]) continue;
        res=dfs(j,min(w[i],f-ss));if(!res) continue;
        w[i]-=res;w[i^1]+=res;cost+=res*ct[i];ss+=res;
        if(ss==f) return f;
    }
    if(!ss) dis[x]=inf;
    return ss;
}

int main(){
    int i,j,ix,iy,iz,iu;
    rd(n);rd(m);rd(S);rd(T);
    for(i=1;i<=m;++i){
        rd(ix);rd(iy);rd(iz);rd(iu);
        lk(ix,iy,iz,iu);
        lk(iy,ix,0,-iu);
    }
    for(;spfa();){
        vis[T]=1;
        for(;vis[T];){
            for(i=1;i<=n;++i) vis[i]=0;
            flow+=dfs(S,inf);
        }
    }
    printf("%d %d",flow,cost);
}