最小费用最大流

最小费最大流

#include<bits/stdc++.h>
#define Ll long long
using namespace std;
const int N=5e3+5;
queue<int>Q;//v是流量，c是花费
struct cs
{
    int to,nxt,v,c;
} a[N*20];
int head[N],ll=1;
int dist[N],flow[N],g[N],xb[N];//flow是到i点的流量，dist是花费，g存spfa路径，xb是下标
bool in[N];
int n,m,S,E,x,y,z,v,ans1,ans2,ans3;
void init(int x,int y,int v,int c)
{
    a[++ll].to=y;
    a[ll].v=v;
    a[ll].c=c;
    a[ll].nxt=head[x];
    head[x]=ll;
}
bool spfa()
{
    memset(dist,0x3f,sizeof(dist));

    int INF=dist[0];
    flow[S]=1e9;
    dist[S]=0;
    Q.push(S);
    while(!Q.empty())
    {
        int x=Q.front();
        Q.pop();
        in[x]=0;
        for(int k=head[x]; k; k=a[k].nxt)
            if(a[k].v>0&&dist[a[k].to]>dist[x]+a[k].c)
            {
                dist[a[k].to]=dist[x]+a[k].c;
                g[a[k].to]=x;
                xb[a[k].to]=k;
                if(!in[a[k].to])in[a[k].to]=1,Q.push(a[k].to);
            }
    }
    return dist[E]!=INF;
}
void work()
{
    while(spfa())
    {
        int sd = 1e9;
         for(int k=E; k!=S; k=a[xb[k]^1].to)//反边就是点也可以直接用g
         {
             sd = min(sd,a[xb[k]].v);
         }
        for(int k=E; k!=S; k=a[xb[k]^1].to)
        {
            a[xb[k]].v-=sd,a[xb[k]^1].v+=sd;
        }
        ans1+=sd;
        ans3+=sd*dist[E];
    }
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&S,&E);
    for(int i=1; i<=m; i++)
    {
        scanf("%d%d%d%d",&x,&y,&z,&v);
        init(x,y,z,v);
        init(y,x,0,-v);
    }
    work();
    printf("%d %d\n",ans1,ans3);
}

![输入
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