【模板】最小费用流

这个板子的细节不少，大概思路就是每次用SPFA找增广路来增广。用SPFA之后可以找的就不仅仅是最短路，所以也会有其他奇奇怪怪的费用问题。

#include<stdio.h>
#include<string.h>
#include<algorithm>
#define inf 1000000000
using namespace std;
int vert,edg,S,T,ansf,ansc,que[100005],dist[100005],pre[5005];
int tot,fr[5005],frr[100005],to[100005],nxt[100005],f[100005],w[100005];
bool inq[100005];
void adde(int p,int q,int F,int W)
{to[tot]=q;frr[tot]=p;f[tot]=F; w[tot]=W;nxt[tot]=fr[p];fr[p]=tot++;}
bool spfa();
int main()
{
    memset(fr,-1,sizeof(fr));
    scanf("%d%d%d%d",&vert,&edg,&S,&T);
    int i,p,q,F,W;
    for(i=1;i<=edg;i++)
    {
        scanf("%d%d%d%d",&p,&q,&F,&W);
        adde(p,q,F,W); adde(q,p,0,-W);
    }
    int minn,now;
    while(spfa())
    {
        minn = inf;
        for(now=T;now!=S;now=frr[pre[now]]) minn = min(minn,f[pre[now]]);
        ansf += minn;
        for(now=T;now!=S;now=frr[pre[now]])
        {
            f[pre[now]]-=minn; f[pre[now]^1] += minn;
            ansc += minn*w[pre[now]];
        }
    }
    printf("%d %d",ansf,ansc);
    return 0;
}
bool spfa()
{
    int i,l=0,r=0,now,t;
    for(i=1;i<=vert;i++) dist[i] = inf;
    memset(inq,false,sizeof(inq));
    dist[S] = 0; que[++r] = S; inq[S] = true;
    while(l<r)
    {
        now = que[++l]; inq[now] = false;
        for(i=fr[now];i!=-1;i=nxt[i])
        {
            t = to[i];
            if(f[i]&&dist[t]>dist[now]+w[i])
            {
                dist[t] = dist[now]+w[i];
                pre[t] = i;
                if(!inq[t])
                {
                    que[++r] = t;
                    inq[t] = true;
                }
            }
        }
    }
    if(dist[T]==inf) return false;    
    return true;
}