最小费用最大流模板

题目描述
如题，给出一个网络图，以及其源点和汇点，每条边已知其最大流量和单位流量费用，求出其网络最大流和在最大流情况下的最小费用。
输入输出格式
输入格式：
第一行包含四个正整数N、M、S、T，分别表示点的个数、有向边的个数、源点序号、汇点序号。
接下来M行每行包含四个正整数ui、vi、wi、fi，表示第i条有向边从ui出发，到达vi，边权为wi（即该边最大流量为wi），单位流量的费用为fi。
输出格式：
一行，包含两个整数，依次为最大流量和在最大流量情况下的最小费用。
输入输出样例
输入样例#1： 
4 5 4 3
4 2 30 2
4 3 20 3
2 3 20 1
2 1 30 9
1 3 40 5
输出样例#1： 
50 280
说明
时空限制：1000ms,128M
（BYX：最后两个点改成了1200ms）
数据规模：
对于30%的数据：N<=10，M<=10
对于70%的数据：N<=1000，M<=1000
对于100%的数据：N<=5000，M<=50000

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1e5+10;
struct node{
    int u,v,fl,c,ne;
}e[N*2];
int h[N],tot,n,m;
void add1(int u,int v,int fl,int c)
{
    tot++;e[tot]=(node){u,v,fl,c,h[u]};h[u]=tot;
}
void add(int u,int v,int fl,int c)
{
    add1(u,v,fl,c);add1(v,u,0,-c);
}
queue<int>q;
int d[N],S,T,pre[N],v[N];
bool spfa()
{
    for(int i=1;i<=n;++i) 
     d[i]=1e9,pre[i]=0,v[i]=0;
    d[S]=0;q.push(S);v[S]=1;
    while(!q.empty())
    {
        int ff=q.front();q.pop();v[ff]=0;
        for(int i=h[ff];i;i=e[i].ne)
        {
            int rr=e[i].v;
            if(e[i].fl && d[rr]>d[ff]+e[i].c)
            {
                 d[rr]=d[ff]+e[i].c;pre[rr]=i;
                 if(!v[rr]) q.push(rr),v[rr]=1;
            }
        }
    }
    return d[T]<1e9;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&S,&T);
    tot=1;
    for(int i=1,x,y,z,w;i<=m;++i)
    {
        scanf("%d%d%d%d",&x,&y,&w,&z);
        add(x,y,w,z);
    }
    ll ans1=0,ans2=0;
    while(spfa())
    {
        int Min=1e9;
        for(int i=T;i!=S;i=e[pre[i]].u)
         Min=min(Min,e[pre[i]].fl);
        for(int i=T;i!=S;i=e[pre[i]].u)
         e[pre[i]].fl-=Min,e[pre[i]^1].fl+=Min;
        ans1+=Min;ans2+=1ll*Min*d[T];
    }
    cout<<ans1<<" "<<ans2;
    return 0;
}

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1e5+10;
struct node{
    int u,v,fl,c,ne;
}e[N*2];
int h[N],tot,n,m;
void add1(int u,int v,int fl,int c)
{
    tot++;e[tot]=(node){u,v,fl,c,h[u]};h[u]=tot;
}
void add(int u,int v,int fl,int c)
{
    add1(u,v,fl,c);add1(v,u,0,-c);
}
queue<int>q;
int d[N],S,T,H[N],v[N];
bool spfa()
{
    for(int i=1;i<=n;++i) d[i]=1e9,v[i]=0;
    d[S]=0;q.push(S);v[S]=1;
    while(!q.empty())
    {
        int ff=q.front();q.pop();v[ff]=0;
        for(int i=h[ff];i;i=e[i].ne)
        {
            int rr=e[i].v;
            if(e[i].fl && d[rr]>d[ff]+e[i].c)
            {
                 d[rr]=d[ff]+e[i].c;
                 if(!v[rr]) q.push(rr),v[rr]=1;
            }
        }
    }
    return d[T]<1e9;
}
int dfs(int u,int fl)
{
    v[u]=1;
    if(u==T || fl==0) return fl;
    int get=0,f;
    for(int i=H[u];i;i=e[i].ne)
    {
        int rr=e[i].v;
        if(e[i].fl && d[rr]==d[u]+e[i].c && !v[rr])
        {
            f=dfs(rr,min(fl,e[i].fl));
            if(!f) continue;
            get+=f;fl-=f;
            e[i].fl-=f;e[i^1].fl+=f;
            H[u]=i;
            if(fl==0) break;
        }
    }
    if(get==0) d[u]=1e9;
    return get;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&S,&T);
    tot=1;
    for(int i=1,x,y,z,w;i<=m;++i)
    {
        scanf("%d%d%d%d",&x,&y,&w,&z);
        add(x,y,w,z);
    }
    ll ans1=0,ans2=0;
    while(spfa())
    {
        v[T]=1;
        while(v[T])
        {
            for(int i=1;i<=n;++i) H[i]=h[i],v[i]=0;
            int nw=dfs(S,1e9);
            ans1+=nw;ans2+=d[T]*nw;
        }
    }
    cout<<ans1<<" "<<ans2;
    return 0;
}