BZOJ3672/UOJ7 [Noi2014]购票
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作。
本文作者:ljh2000
作者博客:http://www.cnblogs.com/ljh2000-jump/
转载请注明出处,侵权必究,保留最终解释权!
Description
Input
第 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。
Output
输出包含 n-1 行,每行包含一个整数。其中第 v 行表示从城市 v+1 出发,到达SZ市最少的购票费用。同样请注意:输出不包含编号为 1 的SZ市。
Sample Input
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
Sample Output
150
70
149
300
150
HINT
对于所有测试数据,保证 0≤pv≤106,0≤qv≤1012,1≤fv<v;保证 0<sv≤lv≤2×1011,且任意城市到SZ市的总路程长度不超过 2×1011。
输入的 t 表示数据类型,0≤t<4,其中:
当 t=0 或 2 时,对输入的所有城市 v,都有 fv=v-1,即所有城市构成一个以SZ市为终点的链;
当 t=0 或 1 时,对输入的所有城市 v,都有 lv=2×1011,即没有移动的距离限制,每个城市都能到达它的所有祖先;
当 t=3 时,数据没有特殊性质。
n=2×10^5
正解:CDQ分治+斜率优化DP+树分治
解题报告:
这道题是三算法合一经典好题...
我们考虑如果是在一个序列上,就是simple的斜率优化裸题了。但是P[i]的变化不随i变化而规律性变化,所以我每次要求的最优斜率是不一样的,维护好下凸包,每次在凸包上二分找到最优值即可。
对于树上的距离限制,我们肯定对于每个点要重建凸包,如果每次重建显然复杂度无法承受,所以我们需要改变节点的处理顺序,使得祖先上的点顺次加入并且不用重构,这很简单,只要考虑每个点能被最上面的哪些点更新到,排个序就可以扫一遍做完了。
另外树分治是为了保证我整体上的处理复杂度是log次,CDQ分治的最大精髓就是先递归处理一部分,再用处理完的部分来更新未处理的部分,显然当祖先节点递归处理完之后,就可以用来处理子树,而每棵子树内部的点都共用了一条链上的祖先,更新答案,最后再递归处理就可以了。
说的有点抽象,看看代码理解起来还是挺快的。其实是我懒,不想写公式...
//It is made by ljh2000
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
typedef long long LL;
const int MAXN = 200011;
const LL inf = (1LL<<62);
int n,f[MAXN],ecnt,first[MAXN],to[MAXN],next[MAXN],size[MAXN],maxS[MAXN],cnt,dui[MAXN];
LL p[MAXN],q[MAXN],lim[MAXN],dis[MAXN],w[MAXN],dp[MAXN];
bool use[MAXN];
struct node{ LL val; int id; }a[MAXN];
inline bool cmp(node q,node qq){ return q.val>qq.val; }
inline void link(int x,int y,LL z){ next[++ecnt]=first[x]; first[x]=ecnt; to[ecnt]=y; w[ecnt]=z; }
inline LL upd(int i,int j){ return dp[j]+(dis[i]-dis[j])*p[i]+q[i]; }
inline LL K(int x,int y){ return (dp[y]-dp[x])/(dis[y]-dis[x]); }
inline int getint(){
int w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
} inline LL getLL(){
LL w=0; LL q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
} inline void dfs(int x){//预处理每个节点的dis和size
size[x]=1;
for(int i=first[x];i;i=next[i]) {
dis[to[i]]=dis[x]+w[i];
dfs(to[i]); size[x]+=size[to[i]];
}
} inline void find_root(int x,int S,int &rt){
maxS[x]=0; size[x]=1;
for(int i=first[x];i;i=next[i]) {
int v=to[i]; if(use[v]) continue;
find_root(v,S,rt);
size[x]+=size[v]; maxS[x]=max(maxS[x],size[v]);
}
maxS[x]=max(maxS[x],S-size[x]);
if(maxS[x]<maxS[rt] && size[x]>1/*不然会GG!!!*/) rt=x;
} inline void dfs2(int x){
a[++cnt].id=x; a[cnt].val=dis[x]-lim[x];
for(int i=first[x];i;i=next[i]) if(!use[to[i]]) dfs2(to[i]);
} inline void solve(int x,int S){
if(S==1) return ; int rt=0,now;
find_root(x,S,rt);
for(int i=first[rt];i;i=next[i]) use[to[i]]=1;//把重心的儿子节点堵上
//CDQ分治
//先处理除了重心的子树之外的部分
solve(x,S-size[rt]+1);//重心也放入这个部分,方便以后讨论 cnt=0; for(int i=first[rt];i;i=next[i]) dfs2(to[i]);
sort(a+1,a+cnt+1,cmp);//把子树内所有点按能被更新到的最高高度自大往小排序
now=rt;//需要向外更新的是,重心到当前分治根节点的这条链上的所有点
int tail,l,r,mid,pos; tail=0;
for(int i=1;i<=cnt;i++) {
while(now!=f[x] && dis[a[i].id]-lim[a[i].id]<=dis[now]) {
while(tail>1 && K(dui[tail],now)>=K(dui[tail-1],dui[tail])) tail--;
//因为dis[now]递减,相当于是横坐标递减,可以看做是倒着加入,注意斜率判断方向
dui[++tail]=now; now=f[now];
}
if(tail>0) {
//在凸包上二分求得最优解
l=1; r=tail; pos=1;
while(l<=r) {
mid=(l+r)>>1; if(mid==tail) { pos=tail; break; }
if(K(dui[mid],dui[mid+1])>=p[a[i].id]) l=mid+1,pos=mid+1;//当前最优的应该是mid+1!!!
else r=mid-1;
}
dp[a[i].id]=min(dp[a[i].id] , upd(a[i].id,dui[pos]));
}
}
for(int i=first[rt];i;i=next[i]) solve(to[i],size[to[i]]);
} inline void work(){
n=getint(); LL x; x=getint();
for(int i=2;i<=n;i++) {
f[i]=getint(); x=getLL(); link(f[i],i,x);
p[i]=getint(); q[i]=getLL(); lim[i]=getLL();
}
dfs(1); maxS[0]=n+1; for(int i=2;i<=n;i++) dp[i]=inf;
solve(1,size[1]);
for(int i=2;i<=n;i++) printf("%lld\n",dp[i]);
} int main()
{
work();
return 0;
}
BZOJ3672/UOJ7 [Noi2014]购票的更多相关文章
- 【BZOJ3672】[Noi2014]购票 树分治+斜率优化
[BZOJ3672][Noi2014]购票 Description 今年夏天,NOI在SZ市迎来了她30周岁的生日.来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. ...
- 【bzoj3672】[Noi2014]购票 斜率优化dp+CDQ分治+树的点分治
题目描述 给出一棵以1为根的带边权有根树,对于每个根节点以外的点$v$,如果它与其某个祖先$a$的距离$d$不超过$l_v$,则可以花费$p_vd+q_v$的代价从$v$到$a$.问从每个点到1花费 ...
- [BZOJ3672][UOJ#7][NOI2014]购票
[BZOJ3672][UOJ#7][NOI2014]购票 试题描述 今年夏天,NOI在SZ市迎来了她30周岁的生日.来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. ...
- bzoj千题计划251:bzoj3672: [Noi2014]购票
http://www.lydsy.com/JudgeOnline/problem.php?id=3672 法一:线段树维护可持久化单调队列维护凸包 斜率优化DP 设dp[i] 表示i号点到根节点的最少 ...
- [BZOJ3672][Noi2014]购票 斜率优化+点分治+cdq分治
3672: [Noi2014]购票 Time Limit: 30 Sec Memory Limit: 512 MBSubmit: 1749 Solved: 885[Submit][Status][ ...
- bzoj 3672: [Noi2014]购票 树链剖分+维护凸包
3672: [Noi2014]购票 Time Limit: 30 Sec Memory Limit: 512 MBSubmit: 480 Solved: 212[Submit][Status][D ...
- BZOJ 3672: [Noi2014]购票( 树链剖分 + 线段树 + 凸包 )
s弄成前缀和(到根), dp(i) = min(dp(j) + (s(i)-s(j))*p(i)+q(i)). 链的情况大家都会做...就是用栈维护个下凸包, 插入时暴力弹栈, 查询时就在凸包上二分/ ...
- [NOI2014]购票 --- 斜率优化 + 树形DP + 数据结构
[NOI2014]购票 题目描述 今年夏天,NOI在SZ市迎来了她30周岁的生日. 来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. 全国的城市构成了一棵以SZ市为根的有根树,每 ...
- 【BZOJ 3672】 3672: [Noi2014]购票 (CDQ分治+点分治+斜率优化)**
3672: [Noi2014]购票 Description 今年夏天,NOI在SZ市迎来了她30周岁的生日.来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. 全国 ...
随机推荐
- GridView实现编辑删除
前台界面: <asp:GridView ID=" ForeColor="#333333" AutoGenerateColumns=" OnRowCance ...
- jQuery 遍历 (each、map)
jQuery 遍历,意为“移动”,用于根据其相对于其他元素的关系来“查找”(或选取)HTML 元素.以某项选择开始,并沿着这个选择移动,直到抵达您期望的元素为止. jQuery 遍历函数 jQuery ...
- 关于Springboot中dubbo整合注意的误区(不对之处请指出)
这是我的客户端配置! 这是生产的配置, 首先注意一下 scan 我之前尝试这样的客户端配置 然后 果断客户端不能注册接口 @Reference(version="1.0") ...
- Nuxt使用Vuex
Vuex 是一个专为 Vue.js 应用程序开发的状态管理模式.它采用集中式存储管理应用的所有组件的状态,并以相应的规则保证状态以一种可预测的方式发生变化. 基础知识这里不再重述,学习的话请自行到官网 ...
- MySQL事件的先后
今天闲聊之时 提及MySQL事件的执行,发现一些自己之前没有注意的细节 如果在执行事件过程中,如果insert的存储过程发生意外 会如何 USE iot2; CREATE TABLE aaaa (ti ...
- Python3.6全栈开发实例[026]
27.文件a.txt内容:每一行内容分别为商品名字,价钱,个数.apple 10 3tesla 100000 1mac 3000 2lenovo 30000 3chicken 10 3通过代码,将其构 ...
- Python3.6全栈开发实例[020]
20.判断一个数是否是水仙花数, 水仙花数是一个三位数, 三位数的每一位的三次方的和还等于这个数. 那这个数就是一个水仙花数, 例如: 153 = 1**3 + 5**3 + 3**3 num = i ...
- npm报错npm ERR! A complete log of this run can be found in: npm ERR! C:\Users\Administrator\AppData\Roaming\npm-cache\_logs\2018-03-15T01_48_14_769Z-debug.log
全局更新 npm npm i npm -g 就ok了
- Hexo+yilia添加helper-live2d插件宠物动画,很好玩的哦~~
个人主页:https://www.yuehan.online 现在博客:www.wangyurui.top 安装模块: 博客根目录选择cmd命令窗口或者git bash 输入以下代码,安装插件 操作: ...
- 024-Spring Boot 应用的打包和部署
一.概述 二.手工打包[不推荐] 打包命令:maven clean package 打包并导出依赖:maven clean package dependency:copy-dependencies 1 ...