[NOI2014]购票——斜率优化+树链剖分+线段树
建议到UOJ上去交
题解
一眼\(DP\),先把转移方程写出来
设\(dp[i]\)为从点\(i\)出发到点\(1\)的最小费用,那么存在转移
\]
这个式子看起来可以斜率优化啊,往下化几步,可以得到类似下面的东西:
若\(d[j] < d[k]\),则当\(\frac{dp[j]-dp[k]}{d[j]-d[k]}\geqslant p[i]\)时从\(j\)转移更优,否则从\(k\)转移更优
这表明我们只需要维护一个下凸壳,转移时二分一下就行了
假设这个问题是在序列上且没有距离限制,你就已经可以\(O(nlogn)\)地\(A\)掉它了
加上距离限制时,我们可以拿个线段树维护一下每一小段的凸壳,查询时取个最大值
挪到树上时,只需要上个树剖
托上面两个东西的福,复杂度也变成了\(O(nlog^3n)\)[手动滑稽]
然后就是码码码
#include <bits/stdc++.h>
using namespace std;
#define MAXN 200000
#define vi vector<int>
#define pb push_back
#define ll long long
#define INF 0x7f7f7f7f7f7f7f7f
#define LIM 17
int n, t;
vi G[MAXN + 5];
int summit[MAXN + 5], f[20][MAXN + 5], fa[MAXN + 5], top[MAXN + 5], sz[MAXN + 5], hson[MAXN + 5], id[MAXN + 5], dfn[MAXN + 5], dfn_clk;
ll d[MAXN + 5], s[MAXN + 5], p[MAXN + 5], q[MAXN + 5], l[MAXN + 5], dp[MAXN + 5];
namespace HLD {
void dfs1(int u) {
sz[u] = 1;
f[0][u] = fa[u];
for (int i = 1; i <= LIM; ++i) f[i][u] = f[i - 1][f[i - 1][u]];
for (auto v : G[u]) {
if (v == fa[u]) continue;
d[v] = d[u] + s[v];
dfs1(v);
sz[u] += sz[v];
if (sz[v] > sz[hson[u]]) hson[u] = v;
}
}
void dfs2(int u, int tp) {
top[u] = tp;
dfn[u] = ++dfn_clk;
id[dfn_clk] = u;
if (hson[u]) dfs2(hson[u], tp);
for (auto v : G[u]) {
if (v == fa[u] || v == hson[u]) continue;
dfs2(v, v);
}
}
}
double slope(int x, int y) {
return static_cast<double>(dp[y] - dp[x])/(d[y] - d[x]);
}
void pre() {
for (int i = 2; i <= n; ++i) {
int u = i;
for (int j = LIM; ~j; --j)
if (d[i] - d[f[j][u]] <= l[i])
u = f[j][u];
summit[i] = !u ? 1 : u;
}
}
namespace SegementTree {
#define mid ((l + r) >> 1)
#define lson (o << 1)
#define rson (o << 1 | 1)
vi ch[4 * MAXN + 5];
void addPoint(vi &v, int x) { // 把点x扔进下凸壳
while (v.size() >= 2 && slope(x, v[v.size() - 2]) < slope(v[v.size() - 1], v[v.size() - 2])) v.pop_back();
v.push_back(x);
}
ll get(vi &v, ll p) { // 在凸壳上二分斜率
if (v.size() == 1) return dp[v[0]] - d[v[0]] * p;
ll ret = INF;
int l = 0, r = v.size() - 1, m;
double s1, s2;
while (l <= r) {
m = (l + r) >> 1;
if (m == v.size() - 1) {
s1 = slope(v[m - 1], v[m]);
ret = min(ret, dp[v[m]] - d[v[m]] * p);
if (s1 <= p) l = m + 1;
else r = m - 1;
}
else if (!m) {
s2 = slope(v[m], v[m + 1]);
ret = min(ret, dp[v[m]] - d[v[m]] * p);
if (s2 <= p) l = m + 1;
else r = m - 1;
}
else {
s1 = slope(v[m - 1], v[m]);
s2 = slope(v[m], v[m + 1]);
ret = min(ret, dp[v[m]] - d[v[m]] * p);
if (s1 <= p && s2 <= p) l = m + 1;
else if(s1 <= p && s2 > p) return min(ret, dp[v[m]] - d[v[m]] * p);
else r = m - 1;
}
}
return ret;
}
void insert(int o, int l, int r, int x, int u) { // 插入一个点
addPoint(ch[o], u);
if (l == r) return ;
if (x <= mid) insert(lson, l, mid, x, u);
else insert(rson, mid + 1, r, x, u);
}
ll query(int o, int l, int r, int L, int R, ll p) { // 在那一堆凸壳中找最小值
if (L <= l && r <= R) return get(ch[o], p);
ll ret = INF;
if (L <= mid) ret = min(ret, query(lson, l, mid, L, R, p));
if (R > mid) ret = min(ret, query(rson, mid + 1, r, L, R, p));
return ret;
}
#undef mid
#undef lson
#undef rson
}
using namespace SegementTree;
void update(int u) { // 求点u的DP值
dp[u] = INF;
int x = fa[u];
while (d[top[x]] >= d[summit[u]]) {
dp[u] = min(dp[u], query(1, 1, n, dfn[top[x]], dfn[x], p[u]) + d[u] * p[u] + q[u]);
x = fa[top[x]];
}
if (d[x] >= d[summit[u]]) dp[u] = min(dp[u], query(1, 1, n, dfn[summit[u]], dfn[x], p[u]) + d[u] * p[u] + q[u]);
}
void work(int u) {
if (u != 1) update(u);
insert(1, 1, n, dfn[u], u);
for (auto v : G[u]) {
if (v == fa[u]) continue;
work(v);
}
}
int main() {
scanf("%d%d", &n, &t);
for (int i = 2; i <= n; ++i) {
scanf("%d%lld%lld%lld%lld", &fa[i], &s[i], &p[i], &q[i], &l[i]);
G[fa[i]].pb(i);
}
HLD::dfs1(1);
HLD::dfs2(1, 1);
d[0] = -1;
pre();
work(1);
for (int i = 2; i <= n; ++i) printf("%lld\n", dp[i]);
return 0;
}
[NOI2014]购票——斜率优化+树链剖分+线段树的更多相关文章
- BZOJ 3672[NOI2014]购票(树链剖分+线段树维护凸包+斜率优化) + BZOJ 2402 陶陶的难题II (树链剖分+线段树维护凸包+分数规划+斜率优化)
前言 刚开始看着两道题感觉头皮发麻,后来看看题解,发现挺好理解,只是代码有点长. BZOJ 3672[NOI2014]购票 中文题面,题意略: BZOJ 3672[NOI2014]购票 设f(i)f( ...
- 【bzoj2402】陶陶的难题II 分数规划+树链剖分+线段树+STL-vector+凸包+二分
题目描述 输入 第一行包含一个正整数N,表示树中结点的个数.第二行包含N个正实数,第i个数表示xi (1<=xi<=10^5).第三行包含N个正实数,第i个数表示yi (1<=yi& ...
- 【BZOJ-2325】道馆之战 树链剖分 + 线段树
2325: [ZJOI2011]道馆之战 Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 1153 Solved: 421[Submit][Statu ...
- 【BZOJ2243】[SDOI2011]染色 树链剖分+线段树
[BZOJ2243][SDOI2011]染色 Description 给定一棵有n个节点的无根树和m个操作,操作有2类: 1.将节点a到节点b路径上所有点都染成颜色c: 2.询问节点a到节点b路径上的 ...
- BZOJ2243 (树链剖分+线段树)
Problem 染色(BZOJ2243) 题目大意 给定一颗树,每个节点上有一种颜色. 要求支持两种操作: 操作1:将a->b上所有点染成一种颜色. 操作2:询问a->b上的颜色段数量. ...
- POJ3237 (树链剖分+线段树)
Problem Tree (POJ3237) 题目大意 给定一颗树,有边权. 要求支持三种操作: 操作一:更改某条边的权值. 操作二:将某条路径上的边权取反. 操作三:询问某条路径上的最大权值. 解题 ...
- bzoj4034 (树链剖分+线段树)
Problem T2 (bzoj4034 HAOI2015) 题目大意 给定一颗树,1为根节点,要求支持三种操作. 操作 1 :把某个节点 x 的点权增加 a . 操作 2 :把某个节点 x 为根的子 ...
- HDU4897 (树链剖分+线段树)
Problem Little Devil I (HDU4897) 题目大意 给定一棵树,每条边的颜色为黑或白,起始时均为白. 支持3种操作: 操作1:将a->b的路径中的所有边的颜色翻转. 操作 ...
- Aizu 2450 Do use segment tree 树链剖分+线段树
Do use segment tree Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://www.bnuoj.com/v3/problem_show ...
- 【POJ3237】Tree(树链剖分+线段树)
Description You are given a tree with N nodes. The tree’s nodes are numbered 1 through N and its edg ...
随机推荐
- oracle 的安装 及环境的配置...
前言 最近这段时间和香港同事一起做项目 负责给日本客户做一个产品 使用的是Oracle 数据库 在一开始项目启动时一直报和oracle数据库有关 本机在这之前就已经安装好了 由于使 ...
- 《MIT 6.828 Lab 1 Exercise 3》实验报告
本实验的网站链接:mit 6.828 lab1 Exercise 3. 题目 Exercise 3. Take a look at the lab tools guide, especially th ...
- TIME_WAIT和CLOSE_WAIT的区别
系统上线之后,通过如下语句查看服务器时,发现有不少TIME_WAIT和CLOSE_WAIT. netstat -n | awk '/^tcp/ {++S[$NF]} END {for(a in S) ...
- SQL中的关键词
AS的用法及妙用 https://www.cnblogs.com/zhaotiancheng/p/6692553.html
- ajax提交 的编码小结
今天用ajax的post方法提交数据给webservice,由于有email的缘故发现,注册方法调用总是报错,打开报文一看提交的 发现@符号被变成了%40 , abc@126.com 即abc%401 ...
- dfs/bfs专项训练
A.棋盘问题——poj1321 在一个给定形状的棋盘(形状可能是不规则的)上面摆放棋子,棋子没有区别.要求摆放时任意的两个棋子不能放在棋盘中的同一行或者同一列,请编程求解对于给定形状和大小的棋盘,摆放 ...
- kmp跑两串的最大相同前后缀 HDU2594
题意:http://acm.hdu.edu.cn/showproblem.php?pid=2594 如题. 思路: Next数组记录的是pos位置失配时要跑到哪里,所以最后得再添加一个字符‘#’. 连 ...
- Python 异常处理与反射机制
Python 的创始人为吉多·范罗苏姆(Guido van Rossum).1989年的圣诞节期间,吉多·范罗苏姆为了在阿姆斯特丹打发时间,决心开发一个新的脚本解释程序,作为ABC语言的一种继承.Py ...
- oracle学习笔记03
一:表空间 /* 创建表空间:逻辑单位,通常我们新建一个项目,就会去创建表空间,在表空间中创建用户,用户去创建表. 语法:create tablespace 表空间名字 datafile '文件的路径 ...
- java方法可变参数研究
1 问题引出 (1)缘由 最近在研究如何在项目中引入Redis缓存,于是遇到可变参数这个疑惑点,之前没有好好研究过,为了避免项目后期出现问题. (2)项目相关技术 SpringBoot Redis K ...