【模板】费用流

传送门

//by sdfzchy
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
typedef long long LL;
const int inf=(1<<30),N=100010,mod=1e9+7;
int n,m,ecnt=1,s,t,head[N];
inline int in()
{
    char tmp=getchar();
    int res=0,f=1;
    while((tmp<'0'||tmp>'9')&&tmp!='-')tmp=getchar();
    if(tmp=='-') f=-1,tmp=getchar();
    while(tmp>='0'&&tmp<='9')   res=(res<<1)+(res<<3)+(tmp^48),tmp=getchar();
    return res*f;
}
struct E
{
    int to,nxt,c,w;
}e[N];
void add(int u,int v,int c,int w)
{
    ++ecnt; e[ecnt]=(E){v,head[u],c,w}; head[u]=ecnt;
    ++ecnt; e[ecnt]=(E){u,head[v],0,-w};head[v]=ecnt;
}
int dis[N];
bool vis[N];
bool bfs()
{
    for(int i=1;i<=n;i++) dis[i]=inf,vis[i]=0;
    deque<int> q; q.push_front(s);
    vis[s]=1,dis[s]=0;
    while(!q.empty())
    {
        int u=q.front(); q.pop_front(); vis[u]=0;
        for(int i=head[u];i;i=e[i].nxt)
        {
            int v=e[i].to;
            if(e[i].c>0&&dis[v]>dis[u]+e[i].w)
            {
                dis[v]=dis[u]+e[i].w;
                if(!vis[v])
                {
                    if(q.empty()) q.push_front(v);
                    else if(dis[v]<dis[q.front()]) q.push_front(v);
                    else q.push_back(v);
                    vis[v]=1;
                }
            }           
        } 
    }
    return dis[t]^inf;
}
int allc;
int dfs(int u,int F)
{
    if(vis[u]) return 0;
    if(u==t||!F) return allc+=dis[t]*F,F;
    vis[u]=1;
    int ret=0,tmp;
    for(int i=head[u];i;i=e[i].nxt)
    {
        int v=e[i].to;
        if(dis[v]==dis[u]+e[i].w&&e[i].c)
            if(tmp=dfs(v,min(F,e[i].c)))
            {
                e[i].c-=tmp,e[i^1].c+=tmp;
                ret+=tmp,F-=tmp;
            }
    }
    vis[u]=0;
    return ret;
}

int mcmf()
{
    int ret=0;
    while(bfs())
    ret+=dfs(s,inf);
    return ret;
}

int main()
{
    n=in(),m=in(),s=1,t=n;
    for(int i=1,u,v,c,w;i<=m;i++)
    {
        u=in(),v=in(),c=in(),w=in();
        add(u,v,c,w);
    }
    printf("%d ",mcmf());   
    printf("%d",allc);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值