大体思路是 spfa找能从起点流到终点的最短路,找到不能找了,就找到了最大流,而此时的费用一定是最小费用;
以下给出ac代码:
h[],e[],ne[] 存图, w[],f[],存每条边的流量和单位费用
pre[],last[] 存修改增广路径的点和边
dis[],flow [] 存到点的花费和流量
```
#include
#include
#include
#include
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
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;
}
```
貌似有更快的做法,以后懂了再来复更