【模板】最小费用最大流

【模板】最小费用最大流

题目链接：https://www.luogu.org/problem/P3381

推荐阅读：https://www.luogu.org/problemnew/solution/P3381

题目描述

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

输入格式

第一行包含四个正整数N、M、S、T，分别表示点的个数、有向边的个数、源点序号、汇点序号。

接下来M行每行包含四个正整数ui、vi、wi、fi，表示第i条有向边从ui出发，到达vi，边权为wi（即该边最大流量为wi），单位流量的费用为fi。

输出格式

一行，包含两个整数，依次为最大流量和在最大流量情况下的最小费用。

输入输出样例

输入

4 5 4 3
4 2 30 2
4 3 20 3
2 3 20 1
2 1 30 9
1 3 40 5

输出

50 280

代码：

#include <bits/stdc++.h>
#define inf 0x3f3f3f3f

using namespace std;

const int maxn = 1e5+10;
struct node {
    int to,w,f,nxt;
} e[maxn];
int n,m,s,t,maxflow,mincost;
int dis[maxn],flow[maxn],via[maxn],pre[maxn],last[maxn];
int head[maxn],cnt=0;

void addedge( int u, int v, int f, int w )
{
    e[cnt].to = v;
    e[cnt].f = f;
    e[cnt].w = w;
    e[cnt].nxt = head[u];
    head[u] = cnt++;
}

int spfa()
{
    memset(dis,inf,sizeof(dis));
    memset(flow,inf,sizeof(flow));
    memset(via,0,sizeof(via));
    queue <int> Q;
    Q.push(s); via[s]=1; dis[s]=0; pre[t]=-1;
    while ( !Q.empty() ) {
        int x = Q.front(); Q.pop(); via[x]=0;
        for ( int i=head[x]; i!=-1; i=e[i].nxt ) {
            int y = e[i].to, f=e[i].f, w=e[i].w;
            if ( f && dis[y]>dis[x]+w ) {
                dis[y] = dis[x] + w;
                pre[y] = x;
                last[y] = i;
                flow[y] = min(flow[x],f);
                if ( via[y]==0 ) {
                    Q.push(y); via[y]=1;
                }
            }
        }
    }
    return pre[t]!=-1;
}

void MCMF()
{
    while ( spfa() ) {
        int x = t;
        maxflow += flow[t];
        mincost += flow[t]*dis[t];
        while ( x!=s ) {
            e[last[x]].f -= flow[t];
            e[last[x]^1].f += flow[t];
            x = pre[x];
        }
    }
}

int main()
{
    int x,y,f,w;
    memset(head,-1,sizeof(head));
    cin >> n >> m >> s >> t;
    for ( int i=0; i<m; i++ ) {
        cin>>x>>y>>f>>w;
        addedge( x,y,f,w );
        addedge( y,x,0,-w );
    }
    MCMF();
    cout << maxflow << " " << mincost << endl;

    return 0;
}