[Noi2014]购票
来自FallDream 的博客,未经允许,请勿转载,谢谢。
今年夏天,NOI在SZ市迎来了她30周岁的生日。来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会。
dp式子比较容易列出来 f[i]=min f[j]+q[i]+p[i]*(d[i]-d[j)
考虑斜率优化,发现j比k优,当且仅当p<(f[j]-f[k])/(d[j]-d[k]) 所以考虑维护凸包
可以直接线段树+无旋treap 大力维护凸包
也可以cdq分治上树,每次像点分治那样做,然后用包含根的那个子树更新其它子树。
更新时候按照深度(用于更新的点直接放深度,否则放深度减去limit)排序从下往上做
复杂度nlogn^2
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
#define MN 200000
#define INF 2000000000
using namespace std;
inline ll read()
{
ll x = , f = ; char ch = getchar();
while(ch < '' || ch > ''){ if(ch == '-') f = -; ch = getchar();}
while(ch >= '' && ch <= ''){x = x * + ch - '';ch = getchar();}
return x * f;
} ll q[MN+],p[MN+],l[MN+],dep[MN+],f[MN+];
int n,head[MN+],cnt=,top,size[MN+],b[MN+],fa[MN+],mx[MN+],Mx[MN+];
struct edge{int to,next;ll w;}e[MN*+];
inline double Getslop(int x,int y){return !y?-INF:(double)(f[y]-f[x])/(dep[y]-dep[x]);}
struct data{int x,kind;ll dep;}a[MN+];
struct CH
{
int s[MN+],tp;
void clear() {tp=;}
void ins(int x)
{
while(tp>&&Getslop(s[tp],s[tp-])<Getslop(x,s[tp])) --tp;
s[++tp]=x;
}
int query(ll p)
{
int l=,r=tp,mid,ans;s[tp+]=;
while(l<=r)
{
mid=l+r>>;
if(Getslop(s[mid],s[mid+])<(double)p) ans=mid,r=mid-;
else l=mid+;
}
return s[ans];
}
}H;
inline void ins(int f,int t,ll w){e[++cnt]=(edge){t,head[f],w};head[f]=cnt;} void Pre(int x)
{
size[x]=;mx[x]=;
for(int i=head[x];i;i=e[i].next)
dep[e[i].to]=dep[x]+e[i].w,Pre(e[i].to),size[x]+=size[e[i].to],mx[x]=max(mx[x],size[e[i].to]);
} void GetRt(int x,int tot,int &rt)
{
size[x]=;mx[x]=;
for(int i=head[x];i;i=e[i].next)
if(!b[e[i].to]) GetRt(e[i].to,tot,rt),size[x]+=size[e[i].to],mx[x]=max(mx[x],size[e[i].to]);
Mx[x]=max(mx[x],tot-size[x]);
if(size[x]>&&Mx[x]<Mx[rt]) rt=x;
} void Dfs(int x)
{
a[++top]=(data){x,,dep[x]-l[x]};
for(int i=head[x];i;i=e[i].next)
if(!b[e[i].to]) Dfs(e[i].to);
}
bool cmp(data x,data y){return x.dep==y.dep?x.kind>y.kind:x.dep>y.dep;}
ll GetAns(int x,int y){return f[y]+p[x]*(dep[x]-dep[y])+q[x];}
void Solve(int x,int Sz)
{
if(Sz<=) return;int rt=;
GetRt(x,Sz,rt);
for(int i=head[rt];i;i=e[i].next) b[e[i].to]=;
Solve(x,Sz-size[rt]+);
a[top=]=(data){rt,,dep[rt]};
for(int t=fa[rt];t!=fa[x];t=fa[t]) a[++top]=(data){t,,dep[t]};
for(int i=head[rt];i;i=e[i].next) Dfs(e[i].to);
sort(a+,a+top+,cmp);H.clear();
for(int i=;i<=top;++i)
if(a[i].kind) H.ins(a[i].x);
else if(H.tp) f[a[i].x]=min(f[a[i].x],GetAns(a[i].x,H.query(p[a[i].x])));
for(int i=head[rt];i;i=e[i].next) Solve(e[i].to,size[e[i].to]);
} int main()
{
n=read();read();
memset(f,,sizeof(f));
for(int i=;i<=n;++i)
{
fa[i]=read();ll w=read();
p[i]=read(),q[i]=read();l[i]=read();
ins(fa[i],i,w);
}
Pre();f[]=;Mx[]=INF;Solve(,size[]);
for(int i=;i<=n;++i) printf("%lld\n",f[i]);
return ;
}
[Noi2014]购票的更多相关文章
- [BZOJ3672][UOJ#7][NOI2014]购票
[BZOJ3672][UOJ#7][NOI2014]购票 试题描述 今年夏天,NOI在SZ市迎来了她30周岁的生日.来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. ...
- 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)). 链的情况大家都会做...就是用栈维护个下凸包, 插入时暴力弹栈, 查询时就在凸包上二分/ ...
- 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][ ...
- [NOI2014]购票 --- 斜率优化 + 树形DP + 数据结构
[NOI2014]购票 题目描述 今年夏天,NOI在SZ市迎来了她30周岁的生日. 来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. 全国的城市构成了一棵以SZ市为根的有根树,每 ...
- 【BZOJ 3672】 3672: [Noi2014]购票 (CDQ分治+点分治+斜率优化)**
3672: [Noi2014]购票 Description 今年夏天,NOI在SZ市迎来了她30周岁的生日.来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. 全国 ...
- 【BZOJ3672】[Noi2014]购票 树分治+斜率优化
[BZOJ3672][Noi2014]购票 Description 今年夏天,NOI在SZ市迎来了她30周岁的生日.来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. ...
- BZOJ_3672_ [Noi2014]购票_CDQ分治+斜率优化
BZOJ_3672_ [Noi2014]购票_CDQ分治+斜率优化 Description 今年夏天,NOI在SZ市迎来了她30周岁的生日.来自全国 n 个城市的OIer们都会从各地出发,到SZ市参 ...
- BZOJ 3672[NOI2014]购票(树链剖分+线段树维护凸包+斜率优化) + BZOJ 2402 陶陶的难题II (树链剖分+线段树维护凸包+分数规划+斜率优化)
前言 刚开始看着两道题感觉头皮发麻,后来看看题解,发现挺好理解,只是代码有点长. BZOJ 3672[NOI2014]购票 中文题面,题意略: BZOJ 3672[NOI2014]购票 设f(i)f( ...
随机推荐
- io多路复用(一)
sever端 1 import socket sk1 = socket.socket() sk1.bind(('127.0.0.1',8001,)) sk1.listen() sk2 = socket ...
- nyoj Color the fence
Color the fence 时间限制:1000 ms | 内存限制:65535 KB 难度:2 描述 Tom has fallen in love with Mary. Now Tom w ...
- MySQL InnoDB锁机制
概述: 锁机制在程序中是最常用的机制之一,当一个程序需要多线程并行访问同一资源时,为了避免一致性问题,通常采用锁机制来处理.在数据库的操作中也有相同的问题,当两个线程同时对一条数据进行操作,为了保证数 ...
- 《javascript设计模式与开发实践》阅读笔记(10)—— 组合模式
组合模式:一些子对象组成一个父对象,子对象本身也可能是由一些孙对象组成. 有点类似树形结构的意思,这里举一个包含命令模式的例子 var list=function(){ //创建接口对象的函数 ret ...
- margin-top导致父标签偏移问题
从一个大神博客中看到这句话: 这个问题发生的原因是根据规范,一个盒子如果没有上补白(padding-top)和上边框(border-top),那么这个盒子的上边距会和其内部文档流中的第一个子元素的上边 ...
- jq 滚轮监听事件
windowAddMouseWheel(); function windowAddMouseWheel() { var i = 0; var scrollFunc = function (e) { e ...
- 你考虑清楚了吗就决定用 Bootstrap ?
近年来,在前端项目中, Bootstrap 已经成为了一个非常受欢迎的工具. Bootstrap 的确有很多优点,然而,如果你的团队中恰好有一个专职的前端工程师.那我推荐你们不要使用 Bootstra ...
- emqtt 试用(六)系统主题
$SYS-系统主题 EMQ 消息服务器周期性发布自身运行状态.MQTT 协议统计.客户端上下线状态到 $SYS/ 开头系统主题. $SYS 主题路径以 "$SYS/brokers/{node ...
- Linux CentOS7.0 (02)修改主机名和ip地址
一.主机名修改 1.查看命令 在CentOS中,有三种定义的主机名:静态的(static),瞬态的(transient),和灵活的(pretty). "静态"主机名也称为内核主机名 ...
- python中两种方法实现二分法查找,细致分析二分法查找算法
之前分析了好多排序算法,可难理解了呢!!(泣不成声)这次我要把二分查找总结一下,这个算法不算难度特别大,欢迎大家参考借鉴我不喜欢太官方的定义,太晦涩的语言,让人看了就头晕.我希望加入我自己的理解,能帮 ...