洛谷 p3381 最小费最大流 模板题 spfa算法

大体思路是 spfa找能从起点流到终点的最短路，找到不能找了，就找到了最大流，而此时的费用一定是最小费用；

以下给出ac代码：

h[],e[],ne[] 存图， w[],f[],存每条边的流量和单位费用

pre[],last[] 存修改增广路径的点和边

dis[],flow [] 存到点的花费和流量 

 

```

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>

using namespace std;
const int maxn=5010,maxm=1e5+7;
int h[maxn],e[maxm],ne[maxm],w[maxm],f[maxm],idx;
int pre[maxn],last[maxn],dis[maxn],flow[maxn];
bool vis[maxn];
int n,m,s,t,a,b,c,d,maxflow,mincost;

void add(int a,int b,int c,int d){
    e[idx]=b,w[idx]=c,f[idx]=d,ne[idx]=h[a],h[a]=idx++;
}
bool spfa(int s,int t){
    memset(dis,0x3f,sizeof dis);
    memset(flow,0x3f,sizeof flow);
    memset(vis,0,sizeof vis);
    pre[t]=-1;
    queue<int>q;q.push(s);vis[s]=1;dis[s]=0;
    while(q.size()){
        int u=q.front();q.pop();vis[u]=0;
        for(int i=h[u];~i;i=ne[i]){
            int j=e[i];
            if(dis[j]>dis[u]+f[i]&&w[i]>0){
                dis[j]=dis[u]+f[i];
                flow[j]=min(w[i],flow[u]);
                pre[j]=u;
                last[j]=i;
                if(!vis[j])q.push(j),vis[j]=1;
            }
        }
    }
    return pre[t]!=-1;
}
void mcmf(){
    while(spfa(s,t)){
        int p=t;
        maxflow+=flow[t];
        mincost+=flow[t]*dis[t];
        while(p!=s){
            w[last[p]]-=flow[t];
            w[last[p]^1]+=flow[t];
            p=pre[p];
        }
    }
}
int main()
{
    memset(h,-1,sizeof h);
    scanf("%d%d%d%d",&n,&m,&s,&t);
    while(m--){
        scanf("%d%d%d%d",&a,&b,&c,&d);
        add(a,b,c,d);add(b,a,0,-d);
    }
    mcmf();
    printf("%d %d\n",maxflow,mincost);
    return 0;
}

```
貌似有更快的做法，以后懂了再来复更