网络流学习(最小费用最大流)

知道了最大流之后这个没什么难度
但是好像全世界都是用的EK算法(SPFA)，要么就来个Dijkstra
好像是说可以证明Dinic跑费用流并没有优化，而且EK算法在费用流里更容易码也更容易理解......
于是乎
luogu模板最小费用最大流

code

换了一种存反边的方式......
话不多说

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iomanip>
#include<algorithm>
#include<ctime>
#include<queue>
#include<stack>
#include<vector>
#define rg register
#define il inline
#define lst long long
#define ldb long double
#define N 5050
#define M 50050
using namespace std;
const int Inf=1e9;
il int read()
{
    rg int s=0,m=0;rg char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')m=1;ch=getchar();}
    while(ch>='0'&&ch<='9')s=(s<<3)+(s<<1)+(ch^48),ch=getchar();
    return m?-s:s;
}

int n,m,S,T;
lst ans,tot;
int hd[N],cnt=1;
int pre[N],vis[N];
int dis[N];
queue<int> Q;
struct EDGE{int to,nxt,c,cost;}ljl[M<<1];
il void add(rg int p,rg int q,rg int o,rg int w)
{
    ljl[++cnt]=(EDGE){q,hd[p],o, w},hd[p]=cnt;
    ljl[++cnt]=(EDGE){p,hd[q],0,-w},hd[q]=cnt;
}

il int SPFA()
{
    memset(dis,63,sizeof(dis));
    while(!Q.empty())Q.pop();
    dis[S]=0,Q.push(S),vis[S]=1;
    while(!Q.empty())
    {
        rg int now=Q.front();Q.pop();vis[now]=0;
        for(rg int i=hd[now];i;i=ljl[i].nxt)
        {
            rg int qw=ljl[i].to;
            if(ljl[i].c&&dis[qw]>dis[now]+ljl[i].cost)
            {
                dis[qw]=dis[now]+ljl[i].cost;
                pre[qw]=i;
                if(!vis[qw])vis[qw]=1,Q.push(qw);
            }
        }
    }
    if(dis[T]==dis[0])return 0;
    rg int flow=Inf;
    for(rg int i=T;i!=S;i=ljl[pre[i]^1].to)
        flow=min(flow,ljl[pre[i]].c);
    tot+=flow,ans+=flow*dis[T];
    for(rg int i=T;i!=S;i=ljl[pre[i]^1].to)
        ljl[pre[i]].c-=flow,ljl[pre[i]^1].c+=flow;
    return 1;
}

int main()
{
    n=read(),m=read(),S=read(),T=read();
    for(rg int i=1;i<=m;++i)
    {
        rg int p=read(),q=read(),o=read(),w=read();
        add(p,q,o,w);
    }
    while(SPFA());
    printf("%lld %lld\n",tot,ans);
    return 0;
}