bzoj1877 SDOI2009晨跑 费用流

25 篇文章 0 订阅
4 篇文章 0 订阅

复习的第三道费用流,拆点以后直接跑。
话说现在省选季也快结束了,弥漫着一股悲伤的气氛(雾)
认识的dalao好像都考挂了(大雾)没考挂的发挥好像也不是很好

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e5+5;
const int inf=0x7fffffff/3;
int S,T;
int n,m,k,ans1,ans2;
int dis[N],from[N],tot=1,head[N];
int q[N];
struct node
{
    int from,to,next,v,c;
}e[N<<1];
bool vis[N];
inline void add(int x,int y,int w,int c)
{
    e[++tot].from=x;
    e[tot].to=y;
    e[tot].v=w;
    e[tot].c=c;
    e[tot].next=head[x];
    head[x]=tot;
}
inline void ins(int x,int y,int w,int c)
{
    add(x,y,w,c);
    add(y,x,0,-c);
}
bool spfa()
{
    int t=0,w=1;
    fo(i,S,T)dis[i]=inf;
    q[1]=S;
    dis[S]=0;
    vis[S]=1;
    while (t<w)
    {
        int x=q[++t];
        for(int i=head[x];i;i=e[i].next)
        {
            int v=e[i].to;
            if (e[i].v&&e[i].c+dis[x]<dis[v])
            {
                dis[v]=dis[x]+e[i].c;
                from[v]=i;
                if (!vis[v])
                {
                    vis[v]=1;
                    q[++w]=v;
                }
            }
        }
        vis[x]=0;
    }
    if (dis[T]==inf)return 0;
    else return 1;
}
inline void mcf()
{
    int x=inf;
    int i=from[T];
    while (i)
    {
        x=min(x,e[i].v);
        i=from[e[i].from];
    }
    ans1++;
    i=from[T];
    while (i)
    {
        ans2+=x*e[i].c;
        e[i].v-=x;
        e[i^1].v+=x;
        i=from[e[i].from];
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    S=1,T=n+n;
    fo(i,1,m)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        ins(x+n,y,1,z);
    }
    fo(i,2,n-1)ins(i,i+n,1,0);
    ins(1,S+n,inf,0);
    ins(n,T,inf,0);
    while (spfa())mcf();
    printf("%d %d\n",ans1,ans2);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值