【bzoj3672】[Noi2014]购票 斜率优化dp+CDQ分治+树的点分治
题目描述
输入
第 1 行包含2个非负整数 n,t,分别表示城市的个数和数据类型(其意义将在后面提到)。输入文件的第 2 到 n 行,每行描述一个除SZ之外的城市。其中第 v 行包含 5 个非负整数 $f_v,s_v,p_v,q_v,l_v$,分别表示城市 v 的父亲城市,它到父亲城市道路的长度,票价的两个参数和距离限制。请注意:输入不包含编号为 1 的SZ市,第 2 行到第 n 行分别描述的是城市 2 到城市 n。
输出
输出包含 n-1 行,每行包含一个整数。其中第 v 行表示从城市 v+1 出发,到达SZ市最少的购票费用。同样请注意:输出不包含编号为 1 的SZ市。
样例输入
7 3
1 2 20 0 3
1 5 10 100 5
2 4 10 10 10
2 9 1 100 10
3 5 20 100 10
4 4 20 0 10
样例输出
40
150
70
149
300
150
题解
斜率优化dp+CDQ分治+树的点分治
设$f[x]$表示$x$到根的最小代价,那么有dp方程:$f[i]=min\{f[j]+(deep[i]-deep[j])·p[i]+q[i]\}$,其中$j$是$i$的祖先且$deep[i]-deep[j]\le l[i]$。
考虑将这个dp方程变形,得:$f[j]=p[i]·deep[j]+f[i]-p[i]·deep[i]-q[i]$。
这很明显是y=kx+b的形式,因此可以使用斜率优化来解决。其中$f$是y,$p$是k,$deep$是x,$f-p·deep-q$是b。
然而dp在树上进行。。。
由于问题在树上,因此不能直接维护凸包(貌似网上有种“可持久化凸包”的写法,然而不会写。。。 以及一种树剖套凸包的写法,然而是3个log的。。。)
考虑CDQ分治,一条深度递减的链,影响的范围是链底的子树中的所有点(不包括链底端点)。因此可以先选定某个点作为分治中心,先处理其上边的子树,即包含树根部分的子树(如果存在)(CDQ分治处理左区间),处理出其到当前树根的链的凸包。然后对于所有底端子树的点,在凸包中查找(CDQ分治处理左边对右边的影响)。由于斜率不是单调的,因此需要二分查找。最后递归处理其它子树(CDQ分治处理右区间)
这个点的选择依据:每个子树的大小不超过树的1/2,显然选择重心(其实这个过程是树的点分治)。
具体的一些细节:
由于更新的点是不包括重心的,所以如果重心不为根节点,则需要先使用重心以上到根节点的链取更新重心,然后再更新其它节点。
对于题目中的$l_v$的限制,可以先将所有底端的点(右区间)按照$deep-l$,即最小能够接受的深度从大到小排序,然后从下往上加入链的节点时,看有多少点对应的凸包恰好为当前凸包(即下一个点不满足条件)。这些点在当前的凸包中二分更新,然后再维护凸包。
由于链上点的$deep$是自带有序的,因此可以直接使用单调栈维护凸包。
最后再递归处理其它子树即可。
时间复杂度为$O(n\log^2n)$。
代码的solve中$x$代表当前树的根,而$rt$代表当前树的重心,不要弄混。
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 200010
using namespace std;
typedef long long ll;
int head[N] , to[N << 1] , next[N << 1] , cnt , fa[N] , p[N];
int si[N] , ms[N] , sum , root , vis[N] , A[N] , ta , B[N] , tb , sta[N] , tot;
ll deep[N] , q[N] , l[N] , f[N];
inline void add(int x , int y)
{
to[++cnt] = y , next[cnt] = head[x] , head[x] = cnt;
}
bool cmp(const int &a , const int &b)
{
return deep[a] - l[a] > deep[b] - l[b];
}
void getroot(int x , int pre)
{
int i;
si[x] = 1 , ms[x] = 0;
for(i = head[x] ; i ; i = next[i])
if(!vis[to[i]] && to[i] != pre)
getroot(to[i] , x) , si[x] += si[to[i]] , ms[x] = max(ms[x] , si[to[i]]);
ms[x] = max(ms[x] , sum - si[x]);
if(ms[x] < ms[root]) root = x;
}
void fill(int x , int pre)
{
int i;
B[++tb] = x;
for(i = head[x] ; i ; i = next[i])
if(!vis[to[i]] && to[i] != pre)
fill(to[i] , x);
}
inline double slop(int a , int b)
{
return (double)(f[a] - f[b]) / (deep[a] - deep[b]);
}
inline void update(int x)
{
if(!tot) return;
int l = 1 , r = tot - 1 , mid , ret = tot;
while(l <= r)
{
mid = (l + r) >> 1;
if(slop(sta[mid] , sta[mid + 1]) < p[x]) ret = mid , r = mid - 1;
else l = mid + 1;
}
f[x] = min(f[x] , f[sta[ret]] - deep[sta[ret]] * p[x] + q[x]);
}
void solve(int x)
{
int i , j , rt;
sum = si[x] , root = 0 , getroot(x , 0) , vis[rt = root] = 1;
if(x != rt) si[x] -= si[rt] , solve(x);
tot = ta = tb = 0;
A[++ta] = rt;
for(i = rt ; i != x ; i = fa[i])
{
if(deep[rt] - l[rt] <= deep[fa[i]]) f[rt] = min(f[rt] , f[fa[i]] - deep[fa[i]] * p[rt] + q[rt]);
A[++ta] = fa[i];
}
for(i = head[rt] ; i ; i = next[i])
if(!vis[to[i]])
fill(to[i] , 0);
sort(B + 1 , B + tb + 1 , cmp);
for(i = j = 1 ; i <= ta ; i ++ )
{
while(j <= tb && deep[A[i]] < deep[B[j]] - l[B[j]]) update(B[j ++ ]);
while(tot > 1 && slop(sta[tot - 1] , sta[tot]) <= slop(sta[tot] , A[i])) tot -- ;
sta[++tot] = A[i];
}
while(j <= tb) update(B[j ++ ]);
for(i = head[rt] ; i ; i = next[i])
if(!vis[to[i]])
solve(to[i]);
}
int main()
{
int n , i;
scanf("%d%*d" , &n);
for(i = 2 ; i <= n ; i ++ )
{
scanf("%d%lld%d%lld%lld" , &fa[i] , &deep[i] , &p[i] , &q[i] , &l[i]);
deep[i] += deep[fa[i]] , q[i] += deep[i] * p[i];
add(fa[i] , i) , add(i , fa[i]);
}
memset(f , 0x3f , sizeof(f)) , f[1] = 0;
ms[0] = 1 << 30 , si[1] = n , solve(1);
for(i = 2 ; i <= n ; i ++ ) printf("%lld\n" , f[i]);
return 0;
}
【bzoj3672】[Noi2014]购票 斜率优化dp+CDQ分治+树的点分治的更多相关文章
- [Noi2014]购票 斜率优化DP+可持久化凸包
貌似网上大部分题解都是CDQ分治+点分治然后再斜率优化DP,我貌似并没有用这个方法. 这一题跟这题有点像,只不过多了一个l的限制 如果说直接跑斜率优化DP,存储整个序列的话,显然是不行的,如图所示(图 ...
- [BZOJ3672][Noi2014]购票 斜率优化+点分治+cdq分治
3672: [Noi2014]购票 Time Limit: 30 Sec Memory Limit: 512 MBSubmit: 1749 Solved: 885[Submit][Status][ ...
- [NOI2014]购票(斜率优化+线段树)
题目描述 今年夏天,NOI在SZ市迎来了她30周岁的生日.来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. 全国的城市构成了一棵以SZ市为根的有根树,每个城市与它的父亲用道路连接 ...
- [NOI2014]购票 --- 斜率优化 + 树形DP + 数据结构
[NOI2014]购票 题目描述 今年夏天,NOI在SZ市迎来了她30周岁的生日. 来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. 全国的城市构成了一棵以SZ市为根的有根树,每 ...
- P4027 [NOI2007]货币兑换(斜率优化dp+cdq分治)
P4027 [NOI2007]货币兑换 显然,如果某一天要买券,一定是把钱全部花掉.否则不是最优(攒着干啥) 我们设$f[j]$为第$j$天时用户手上最多有多少钱 设$w$为花完钱买到的$B$券数 $ ...
- 【BZOJ2149】拆迁队(斜率优化DP+CDQ分治)
题目: 一个斜率优化+CDQ好题 BZOJ2149 分析: 先吐槽一下题意:保留房子反而要给赔偿金是什么鬼哦-- 第一问是一个经典问题.直接求原序列的最长上升子序列是错误的.比如\(\{1,2,2,3 ...
- HDU 3824/ BZOJ 3963 [WF2011]MachineWorks (斜率优化DP+CDQ分治维护凸包)
题面 BZOJ传送门(中文题面但是权限题) HDU传送门(英文题面) 分析 定义f[i]f[i]f[i]表示在iii时间(离散化之后)卖出手上的机器的最大收益.转移方程式比较好写f[i]=max{f[ ...
- bzoj 1492: [NOI2007]货币兑换Cash【贪心+斜率优化dp+cdq】
参考:http://www.cnblogs.com/lidaxin/p/5240220.html 虽然splay会方便很多,但是懒得写,于是写了cdq 首先要想到贪心的思路,因为如果在某天买入是能得到 ...
- [NOI2014]购票——斜率优化+树链剖分+线段树
建议到UOJ上去交 题解 一眼\(DP\),先把转移方程写出来 设\(dp[i]\)为从点\(i\)出发到点\(1\)的最小费用,那么存在转移 \[f[i]=min\{f[j]+(d[i]-d[j]) ...
随机推荐
- 生产环境 tidb部署实践
TiDB 简介 TiDB 是 PingCAP 公司受 Google Spanner / F1 论文启发而设计的开源分布式 HTAP (Hybrid Transactional and Analytic ...
- spring源码-Aware-3.4
一.Aware接口,这个也是spring的拓展之一,为啥要单独拿出来讲呢,因为他相比于BeanFactoryPostProcessor,BeanPostProcessor的实用性更加高,并且在具体的业 ...
- Drupal 出错的解决办法
今天安装了superfish菜单模块,安装了一个新菜单后.网站突然打不开了.空白! 第一反应看日志,Apache服务器日志没有发现异常. 可以肯定是添加菜单时,在ATTACH BLOCK部分的区块区域 ...
- NB-IOT使用LWM2M移动onenet基础通信套件对接之APN设置
1. 先搞懂APN是做什么的?APN指一种网络接入技术,是通过手机上网时必须配置的一个参数,它决定了手机通过哪种接入方式来访问网络.对于手机用户来说,可以访问的外部网络类型有很多,例如:Interne ...
- libevent学习七(bufferevent)
1. 每个bufferevent 都拥有类型为struct evbuffer的input buffer和out buffer,分别供数据读取和数据写入使用. 2.读取和写入数据是通过编写和设置对应的回 ...
- 刚安装的Linux Centos7使用yum安装firefox时提示:cannot find a valid baseurl for repo
出现这个问题是因为yum在安装包的过程中,虽然已经联网,但是没法解析远程包管理库对应的域名,所以我们只需要在网络配置中添加上DNS对应的ip地址即可. 解决参考链接:https://blog.csdn ...
- Selenium(Python)PageObject页面对象
使用PageObject页面对象的好处是, 当页面元素的位置发生改变时, 只需要去修改Xpath或者ID, 而不用去修改测试用例本身: 本次的思路是: 1.常用方法类 2.页面对象类 3.测试用例类 ...
- selenium自动化测试资源整理
1. 所有版本chrome下载 是不是很难找到老版本的chrome?博主收集了几个下载chrome老版本的网站,其中哪个下载的是原版的就不得而知了. http://www.slimjet.com/ch ...
- Python学习笔记(一)一一一环境安装错误总结
第三方库安装 1 windows存在多个版本的python,pip安装Python库失败 解决方案:进入对应官网下载安装包,步骤:1 下载安装包到C:\Python36\Lib\site-pack ...
- 【SpringCloud 】第八篇: 消息总线(Spring Cloud Bus)
前言: 必需学会SpringBoot基础知识 简介: spring cloud 为开发人员提供了快速构建分布式系统的一些工具,包括配置管理.服务发现.断路器.路由.微代理.事件总线.全局锁.决策竞选. ...