码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj3876

时间:2017-08-12 10:27:42      阅读:258      评论:0      收藏:0      [点我收藏+]

标签:clu   std   bzoj   style   ems   amp   上界   pac   src   

有上下界的费用流

建图:每条边连接两个节点,边的费用为花费时间,下界为1,上界正无穷,源点为1,所有点连向1,下界0,上界正无穷,表示随时可以返回,于是我们想求一个最小费用可行流

具体做法是先建立超级源汇,对于每条有上下界的边(u,v),u->t,流量为u的入度,费用为0,s->v,流量为下界1,费用为(u,v)费用,u->v连边,流量inf,花费为(u,v)的花费,这样强迫流量流入流出,流满了下界

其他边和之前一样连

u->t的目的是让流入u的流量至少流过到v的边的下界,这样不会让流到u的流量小于下界,s->v的边的目的是强制流入v点至少下界的流量,花费为原来边的花费,这样保证了对于一条边(u,v),流入u的流量不少于下界,从v流出的流量至少为下界,且花费为(u,v)的花费。保证对于点u,流入u的流量至少为所有u的出边的下界之和,这样流入u,也就是流入u出边的下界可以被满足,然后点v满足了到v边下界的流量肯定能流入v,且花费分别为所有流入v的边

技术分享
#include<bits/stdc++.h>
using namespace std;
const int N = 10010, inf = 1 << 29;
struct edge {
    int nxt, to, f, c;
} e[N << 2];
int n, cnt = 1, source, sink;
int head[N], d[N], pree[N], prev[N], vis[N];
inline void link(int u, int v, int f, int c)
{
    e[++cnt].nxt = head[u];
    head[u] = cnt;
    e[cnt].to = v;
    e[cnt].f = f;
    e[cnt].c = c;
}
inline void insert(int u, int v, int f, int c)
{
    link(u, v, f, c);
    link(v, u, 0, -c);
}
bool spfa()
{
    queue<int> q;
    memset(d, -1, sizeof(d));
    q.push(source);
    d[source] = 0;
    while(!q.empty())
    {
        int u = q.front();
        q.pop();
        vis[u] = 0;
        for(int i = head[u]; i; i = e[i].nxt) if(e[i].f && (d[e[i].to] > d[u] + e[i].c || d[e[i].to] == -1))
        {
            d[e[i].to] = d[u] + e[i].c;
            pree[e[i].to] = i;
            prev[e[i].to] = u;
            if(vis[e[i].to] == 0)
            {
                vis[e[i].to] = 1;
                q.push(e[i].to);
            }
        } 
    }
    return d[sink] != -1;
}
int Edmonds_Karp()
{
    int ret = 0;
    while(spfa())
    {
        int now = sink, delta = inf;
        while(now != source)
        {
            delta = min(delta, e[pree[now]].f);
            now = prev[now];
        }
        now = sink;
        while(now != source)
        {
            e[pree[now]].f -= delta;
            e[pree[now] ^ 1].f += delta;
            now = prev[now];
        }
        ret += d[sink] * delta;
    }
    return ret;
}
int main()
{
    scanf("%d", &n);
    sink = n + 1;
    for(int i = 1; i <= n; ++i)
    {
        int m, b, t;
        scanf("%d", &m);
        insert(i, sink, m, 0);
        while(m--)
        {
            scanf("%d%d", &b, &t);
            insert(source, b, 1, t);
            insert(i, b, inf, t);
        }
        if(i != 1) insert(i, 1, inf, 0);    
    }
    printf("%d\n", Edmonds_Karp());
    return 0;
}
View Code

 

bzoj3876

标签:clu   std   bzoj   style   ems   amp   上界   pac   src   

原文地址:http://www.cnblogs.com/19992147orz/p/7348725.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!