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

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

#include
#define Ll long long
using 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