费用流——洛谷P3381 【模板】最小费用最大流

来源：互联网发布：ucsc数据库编辑：程序博客网时间：2024/05/16 15:02

https://daniu.luogu.org/problem/show?pid=3381
我学习了EK算法之后发现别人都学zkw的我靠；
跑的比我快我靠；
思想就是每次sfpa找到一条可以增广并且费用最小的路；
然后暴力修改这条路上的信息；

#include<bits/stdc++.h>#define Ll long longusing 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;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,63,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;                flow[a[k].to]=min(flow[x],a[k].v);                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()){        for(int k=E;k!=S;k=g[k])a[xb[k]].v-=flow[E],a[xb[k]^1].v+=flow[E];        ans1+=flow[E]; ans2+=flow[E]*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",ans1,ans2);}

zkw费用流模版在这里
http://blog.csdn.net/jzq233jzq/article/details/73123089

阅读全文

2 0