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

P4177 [CEOI2008]order

时间:2018-07-09 21:29:07      阅读:135      评论:0      收藏:0      [点我收藏+]

标签:algorithm   ++   iostream   index   没有   turn   输出   mem   return   

P4177 [CEOI2008]order

题目描述
有N个工作,M种机器,每种机器你可以租或者买过来. 每个工作包括若干道工序,每道工序需要某种机器来完成,你可以通过购买或租用机器来完成。 现在给出这些参数,求最大利润

输入输出格式
输入格式:
第一行给出 N,M(1<=N<=1200,1<=M<=1200) 下面将有N组数据。

每组数据第一行给出完成这个任务能赚到的钱(其在[1,5000])及有多少道工序

接下来若干行每行两个数,分别描述完成工序所需要的机器编号及租用它的费用(其在[1,20000]) 最后M行,每行给出购买机器的费用(其在[1,20000])

输出格式:
最大利润


若是没有租用一说,本题即为最大权闭合子图求最小割。

当有租用的时候,我们仍按照最大权闭合子图构建模型,与之前不同的是,有关系的两点连边容边不再是INF而是租用费用

我们仍可以用之前的解释:要么舍弃报酬(即正权点),要么买器材(负权点);租用怎么办呢?对于一个实验,我们有了除了之前两种方案的(即不做工作或购买机器)的第三种选择,租用,所以在求最小割的时候,我们要么割断源点到正权之间的边(不做这项工作),要么割断负权到汇点之间的边(买机器),要么割断中间的边(在一个工作中租用仪器),求得的最小割即为舍弃的最小利益

所以最优答案即为 总利益 - 舍弃的最小利益

Code

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
typedef long long ll;
using namespace std;
int RD(){
    int out = 0,flag = 1;char c = getchar();
    while(c < ‘0‘ || c >‘9‘){if(c == ‘-‘)flag = -1;c = getchar();}
    while(c >= ‘0‘ && c <= ‘9‘){out = out * 10 + c - ‘0‘;c = getchar();}
    return flag * out;
    }
const int maxn = 1000019,INF = 1e9;
int numw,numm;
int head[maxn],nume = 1;
struct Node{
    int v,dis,nxt;
    }E[maxn << 3];
void add(int u,int v,int dis){
    E[++nume].nxt = head[u];
    E[nume].v = v;
    E[nume].dis = dis;
    head[u] = nume;
    }
int maxflow,s,t;
int d[maxn];
bool bfs(){
    queue<int>Q;
    memset(d,0,sizeof(d));
    d[s] = 1;
    Q.push(s);
    while(!Q.empty()){
        int u = Q.front();Q.pop();
        for(int i = head[u];i;i = E[i].nxt){
            int v = E[i].v;
            if(!d[v] && E[i].dis){
                d[v] = d[u] + 1;
                Q.push(v);
                if(v == t)return 1;
                }
            }
        }
    return 0;
    }
int Dinic(int u,int flow){
    if(u == t)return flow;
    int rest = flow,k;
    for(int i = head[u];i;i = E[i].nxt){
        int v = E[i].v;
        if(d[v] == d[u] + 1 && E[i].dis){
            k = Dinic(v,min(rest,E[i].dis));
            if(!k)d[v] = 0;
            E[i].dis -= k;
            E[i ^ 1].dis += k;
            rest -= k;
            }
        if(!rest)break;
        }
    return flow - rest;
    }
int tot;
int main(){
    s = 0, t = 10026;
    numw = RD();numm = RD();
    for(int i = 1;i <= numw;i++){
        int money = RD(), num = RD();
        tot += money;
        add(s, i << 1, money);add(i << 1, s, 0);
        for(int j = 1;j <= num;j++){
            int index = RD(), rent = RD();
            add(i << 1, index << 1 | 1, rent);
            add(index << 1 | 1, i << 1, 0);
            }
        }
    for(int i = 1;i <= numm;i++){
        int money = RD();
        add(i << 1 | 1, t, money);
        add(t, i << 1 | 1, 0);
        }
    int flow = 0;
    while(bfs())while(flow = Dinic(s,INF))maxflow += flow;
    printf("%d\n",tot - maxflow);
    return 0;
    }

P4177 [CEOI2008]order

标签:algorithm   ++   iostream   index   没有   turn   输出   mem   return   

原文地址:https://www.cnblogs.com/Tony-Double-Sky/p/9285939.html

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