本蒟蒻又来写解题报告了。这次的题目是codevs 1035 火车停留。

题目大意就是给m个火车的到达时间、停留时间和车载货物的价值,车站有n个车道,而火车停留一次车站就会从车载货物价值中获得1%的利润,让你来求一种安排方法,使得车站收益最大,并输出收益值。

蒟蒻的思路是这样的:

一眼看出:最大费用最大流(MCMF)
显然cost:表示车站收益
然后……
以车站为点建立图求流?同一个车站可能经过好几辆火车……,貌似很麻烦……;
那么以什么建图、连边,还有怎么连?
貌似有点类似于方格取数2之中的拆点……;
那么这个就可以……以火车为点,把一个点拆成两个,然后建立流的关系。
Reach[i]和stay[i]作为建立不同火车是否可以建边的判断条件
假设火车i,将其分为点2*i-1和2*i,连一条流量为1,费用为-cost[i]的边
如果reach[i] + stay[i] < reach[j],在2*i和2*j-1之间连一条流量为1,费用为0的边。
建立汇点和起点S,T,从S向每个2*i-1连一条流量为1,费用为0的边。从每个2*i向T连一条流量为1,费用为0的边。

那么怎么限制n个车道呢?其实很简单,只需要建立超级汇点ST,然后从T向ST连一条流量为n,费用为0的边。

这样这个题目就大功告成了,代码不长,刚98行。建议大家还是去刷一下代码能力题,因为NOIP2015蒟蒻就被代码能力题给坑惨了。

废话不多说,上代码

 #include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
using namespace std;
const int INF = ;
const int maxe = ;
const int maxn = ;
int n,m,reach[maxn],stay[maxn],cost[maxn],vis[maxn<<],d[maxn<<],h[maxn<<],pre[maxn<<],rid[maxn],cid[maxn],now;
double ans;
struct edge{
int to,cost,cap,next;
}tr[maxe];
inline void init(){
now = ;
memset(h,-,sizeof(h));
memset(tr,,sizeof(tr));
memset(reach,,sizeof(reach));
memset(stay,,sizeof(stay));
memset(cost,,sizeof(cost));
}
inline void add(int u,int v,int cap,int cost){
tr[now].to = v;tr[now].cap = cap;tr[now].cost = cost;tr[now].next = h[u];
h[u] = now++;
tr[now].to = u;tr[now].cap = ;tr[now].cost = -cost;tr[now].next = h[v];
h[v] = now++;
}
bool SPFA(int s,int t,int &flow,int &cost){
for(int i = s;i <= t;++i){
d[i] = INF;
}
int minflow = INF;
memset(vis,,sizeof(vis));
memset(pre,-,sizeof(pre));
deque<int>q;
d[s] = ;
vis[s] = ;
q.push_back(s);
while(!q.empty()){
int x = q.front();q.pop_front();
vis[x] = ;
for(int i = h[x];i != -;i = tr[i].next){
edge e = tr[i];
if(d[e.to] > d[x] + e.cost && e.cap){
d[e.to] = d[x] + e.cost;
pre[e.to] = i;
minflow = min(minflow,e.cap);
if(!vis[e.to]){
vis[e.to] = ;
q.push_back(e.to);
}
}
}
}
if(d[t] == INF)return false;
flow += minflow;
cost += d[t] * minflow;
for(int i = t;i != s;i = tr[pre[i]^].to){
tr[pre[i]].cap -= minflow;
tr[pre[i]^].cap += minflow;
}
return true;
}
int MCMF(int s,int t){
int flow = ,cost = ;
while(SPFA(s,t,flow,cost));
return cost;
}
int main(){
scanf("%d%d",&n,&m);
init();
for(int i = ;i <= m;++i){
scanf("%d%d%d",&reach[i],&cost[i],&stay[i]);
rid[i] = *i;cid[i] = *i-;
add(cid[i],rid[i],,-cost[i]);
}
for(int i = ;i <= m;++i){
for(int j = ;j <= m;++j){
if(j != i)
if(reach[i] + stay[i] < reach[j]){
add(rid[i],cid[j],,);
}
}
}
int S = ,T = *m+,ST = *m+;
for(int i = ;i <= m;++i){
add(S,cid[i],,);
}
for(int i = ;i <= m;++i){
add(rid[i],T,,);
}
add(T,ST,n,);
ans = double(MCMF(S,ST))/;
printf("%0.2lf\n",-ans);
return ;
}

如果大家认为有用的话,欢迎转载。如果大家有意见的话,请在下方评论栏中给我留言,或者给我的E-mail PC-worker@outlook.com发邮件留言。谢谢大家。

习题:codevs 1035 火车停留解题报告的更多相关文章

  1. wikioi 1035 火车停留 裸费用流

    链接:http://wikioi.com/problem/1035/ 怎么说呢,只能说这个建图很有意思.因为只有m条道,然后能互相接在一起的连通,对每个点进行拆点,很有意思的一道裸费用留题. 代码: ...

  2. 习题:codevs 2822 爱在心中 解题报告

    这次的解题报告是有关tarjan算法的一道思维量比较大的题目(真的是原创文章,希望管理员不要再把文章移出首页). 这道题蒟蒻以前做过,但是今天由于要复习tarjan算法,于是就看到codevs分类强联 ...

  3. 习题: codevs 2492 上帝造题的七分钟2 解题报告

    这道题是受到大犇MagHSK的启发我才得以想出来的,蒟蒻觉得自己的代码跟MagHSK大犇的代码完全比不上,所以这里蒟蒻就套用了MagHSK大犇的代码(大家可以关注下我的博客,友情链接就是大犇MagHS ...

  4. 习题:codevs 1519 过路费 解题报告

    今天拿了这道题目练练手,感觉自己代码能力又增强了不少: 我的思路跟别人可能不一样. 首先我们很容易就能看出,我们需要的边就是最小生成树算法kruskal算法求出来的边,其余的边都可以删掉,于是就有了这 ...

  5. [置顶] 刘汝佳《训练指南》动态规划::Beginner (25题)解题报告汇总

    本文出自   http://blog.csdn.net/shuangde800 刘汝佳<算法竞赛入门经典-训练指南>的动态规划部分的习题Beginner  打开 这个专题一共有25题,刷完 ...

  6. LeetCode: Sort Colors 解题报告

    Sort ColorsGiven an array with n objects colored red, white or blue, sort them so that objects of th ...

  7. CH Round #56 - 国庆节欢乐赛解题报告

    最近CH上的比赛很多,在此会全部写出解题报告,与大家交流一下解题方法与技巧. T1 魔幻森林 描述 Cortana来到了一片魔幻森林,这片森林可以被视作一个N*M的矩阵,矩阵中的每个位置上都长着一棵树 ...

  8. 二模13day1解题报告

    二模13day1解题报告 T1.发射站(station) N个发射站,每个发射站有高度hi,发射信号强度vi,每个发射站的信号只会被左和右第一个比他高的收到.现在求收到信号最强的发射站. 我用了时间复 ...

  9. BZOJ 1051 最受欢迎的牛 解题报告

    题目直接摆在这里! 1051: [HAOI2006]受欢迎的牛 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 4438  Solved: 2353[S ...

随机推荐

  1. 1027 HTML的初学

       HTML   是一种超文本标记语言内容(Hyper Text  Markup Language)    CSS       网页美化    Javascript    脚本语言(JS)    H ...

  2. 破解Excel密码保护文件

    首先打开vba编辑器,输入代码: Public Sub AllInternalPasswords() ' Breaks worksheet and workbook structure passwor ...

  3. 事件异步(EAP)使用事件异步处理一些耗时操作

    比如需要下载一些比较大的文件,如果使用会UI卡顿,使用异步可以节省一些时间 下面是一些例子: using System; using System.Collections.Generic; using ...

  4. KnockoutJS 3.X API 第四章 数据绑定(2) 控制流foreach绑定

    foreach绑定 foreach绑定主要用于循环展示监控数组属性中的每一个元素,一般用于table标签中 假设你有一个监控属性数组,每当您添加,删除或重新排序数组项时,绑定将有效地更新UI的DOM- ...

  5. centos-lynx

    1.官网 http://lynx.isc.org 2.稳定版本 http://invisible-mirror.net/archives/lynx/tarballs/lynx2.8.8rel.2.ta ...

  6. WPF自定义控件与样式(8)-ComboBox与自定义多选控件MultComboBox

    一.前言 申明:WPF自定义控件与样式是一个系列文章,前后是有些关联的,但大多是按照由简到繁的顺序逐步发布的等,若有不明白的地方可以参考本系列前面的文章,文末附有部分文章链接. 本文主要内容: 下拉选 ...

  7. javascript运算符——关系运算符

    × 目录 [1]恒等 [2]相等 [3]大于[4]小于 前面的话 关系运算符用于测试两个值之间的关系,根据关系是否存在而返回true或false,关系表达式总是返回一个布尔值,通常在if.while或 ...

  8. tomcat实现ServletContext的addListener方法的源码解说(原创)

    tomcat 8.0.36 知识点: 动态监听器有七类: ServletContextAttributeListener ServletRequestListener ServletRequestAt ...

  9. [Node.js] Node.js中的流

    原文地址:http://www.moye.me/2015/03/29/streaming_in_node/ 什么是流? 说到流,就涉及到一个*nix的概念:管道——在*nix中,流在Shell中被实现 ...

  10. Spring Setter Injection and Constructor Injection

    Setter Injection AppContext.xml <?xml version="1.0" encoding="UTF-8"?> < ...