\(NOI2014\)购票

哇终于可以碰电脑了赶快切些火题找找感觉。

拿到这道题的时候发现简单的斜率优化推一推可以秒掉平方做法,然后一条链也可以做。

然后呢。。。

卧槽这个在一棵树上怎么办啊。

大力\(YY\)了一个数据结构维护区间凸壳的东西,然而我感觉这东西好恶心啊我不会啊。。。

明明就是因为懒

然后翻了一波题解发现可以用点分治,然而都讲得非常不清楚。。。。

看了好久才看懂怎么做吧。。。我这里详细地记录一下。。

详细是不存在的这辈子的不存在的

总体思想是,每次点分治的时候有两个关键点,一个是分治重心,一个是当前子树的根节点。

每次先将含有当前子树根节点的那个分治子树加上分治重心作为一个全新的子树先递归处理

然后把分治中心到当前根节点的路径上的点抠出来,用这些点来更新其他不包含根节点的分治子树。

这个更新是可以统一处理的,具体操作就是把所有点抠出来按照可到达的最浅点从深到浅排个序

也就是斜率计算中的\(x\)

然后从分治重心一步步向上跳父亲加入单调栈,得到用于更新答案的单调栈

因为斜率不单调所以我们要在单调栈上二分

然后再继续递归不包含根节点的分治子树就完成啦。

这样做相当于从浅到深依次分治,依次把\(DP\)值确定好了。

方便理解的话注意一下这棵树是不需要建反边的。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<set>
#include<queue>
#include<map>
#include<vector> using namespace std; #define N 201000
#define EPS 1e-7
#define ll long long
#define RG register
#define inf 1e18+1 inline ll read(){
RG ll x=0,o=1; RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
if(ch=='-') o=-1,ch=getchar();
while(ch>='0'&&ch<='9') x=((x<<3)+(x<<1))+ch-'0',ch=getchar();
return x*o;
} int n,first[N],top,siz[N],vis[N],tmp[N],tot,Max,Sz,Rt,fa[N],Q[N],cnt;
ll dis[N],lim[N],p[N],q[N],f[N];
struct mona { int nxt,en; } s[N<<1];
inline void Insert(int x,int y) { s[++top]=(mona) { first[x],y },first[x]=top; }
inline void Getroot(int k){
siz[k]=1; int Maxn=0;
for(RG int i=first[k];i;i=s[i].nxt){
int en=s[i].en; if(vis[en]) continue ;
Getroot(en),Maxn=max(Maxn,siz[en]),siz[k]+=siz[en];
} Maxn=max(Maxn,Sz-siz[k]); if(Max>=Maxn) Max=Maxn,Rt=k;
}
inline void Getdis(int k){
tmp[++tot]=k;
for(RG int i=first[k];i;i=s[i].nxt)
if(!vis[s[i].en]) Getdis(s[i].en);
} inline bool cmp(const int &x,const int &y) { return dis[x]-lim[x]>dis[y]-lim[y]; }
inline double Slope(const int &i,const int &j){
if(dis[i]==dis[j]) return f[i]<f[j]?inf:-inf;
return (1.0*f[i]-f[j])/(1.0*dis[i]-dis[j]);
}
inline void Add(int i){
while(cnt>1&&Slope(Q[cnt],Q[cnt-1])<=Slope(i,Q[cnt-1])) --cnt;
Q[++cnt]=i;
}
inline int Binary(int cst){
int l=2,r=cnt,ans=1;
while(l<=r){
int mid=l+r>>1;
if(Slope(Q[mid],Q[mid-1])>=cst) ans=mid,l=mid+1;
else r=mid-1;
} return Q[ans];
} inline void Divide(int k,int Siz){
if(Siz==1) return ; Sz=Siz,Max=1e9,Getroot(k); int Cr=Rt;
for(RG int i=first[Cr];i;i=s[i].nxt) vis[s[i].en]=1;
Divide(k,Siz-siz[Cr]+1),tot=0;
for(RG int i=first[Cr];i;i=s[i].nxt) Getdis(s[i].en);
sort(tmp+1,tmp+1+tot,cmp); int now=Cr; cnt=0,Add(Cr);
for(RG int i=1;i<=tot;++i){ int t=tmp[i];
if(dis[now]<dis[t]-lim[t]) continue ;
while(now!=k&&dis[fa[now]]>=dis[t]-lim[t])
now=fa[now],Add(now);
int Bes=Binary(p[t]);
f[t]=min(f[t],f[Bes]+(dis[t]-dis[Bes])*p[t]+q[t]);
} for(RG int i=first[Cr];i;i=s[i].nxt) Divide(s[i].en,siz[s[i].en]);
} inline void Dfs(int k){
for(RG int i=first[k];i;i=s[i].nxt)
dis[s[i].en]+=dis[k],Dfs(s[i].en);
}
int main(){
n=read(),read();
for(RG int i=2;i<=n;++i)
fa[i]=read(),dis[i]=read(),Insert(fa[i],i),
p[i]=read(),q[i]=read(),lim[i]=read(),f[i]=inf;
Dfs(1),Sz=n,Divide(1,n);
for(RG int i=2;i<=n;++i) printf("%lld\n",f[i]);
}

随机推荐

  1. 4412 使用usb摄像头拍照YUYV格式

    一.内核设置 Linux内核中已经带有很完善的USB摄像头驱动,支持几乎所有的USB摄像头,我们只需要配置内核,选择上相应的Sensor型号即可. 配置内核,支持USB摄像头: Device Driv ...

  2. luoguP4721 【模板】分治 FFT (分治NTT)

    给定 $g[1....n-1]$,求 $f[0],f[1],...,f[n-1]$,其中   $f[i]=\sum_{j=1}^{i}f[i-j]g[j]$    变界为 $f[0]=1$ 答案模 9 ...

  3. js+html实现玫瑰花绽放

    玫瑰花示例demo: <!DOCTYPE html> <html> <head> <title>玫瑰花</title> <meta c ...

  4. spash和selenium浅析

    Splash是什么: Splash是一个Javascript渲染服务.它是一个实现了HTTP API的轻量级浏览器,Splash是用Python实现的,同时使用Twisted和QT.Twisted(Q ...

  5. RequestBody 请求格式

    @RequestMapping("/CreateComment")@ResponseBodypublic String createQuestionComments(@Reques ...

  6. SHADER 用 step替代 if-else

    今天聊起这个问题,百度发现了这个优化方式: https://blog.csdn.net/liu_if_else/article/details/77455639

  7. 【CDN+】Kafka 的初步认识与入门

    前言 项目中用到了Kafka 这种分布式消息队列来处理日志,本文将对Kafka的基本概念和原理做一些简要阐释 Kafka 的基本概念 官网解释: Kafka是最初由Linkedin公司开发,是一个分布 ...

  8. CSS-格式化上下文(Formatting Context)

    Formatting Context:指页面中的一个渲染区域,并且拥有一套渲染规则,他决定了其子元素如何定位,以及与其他元素的相互关系和作用. BFC 块级格式化上下文,它是指一个独立的块级渲染区域, ...

  9. Chrome开发小技巧--浏览器控制台现写并运行js代码--snippets

    想简单等运行一段js代码,以前可能会新建一个html 里面包含script标签,或者引入一个js,然后chrome浏览器打开.这样很麻烦. 想再console控制台写,也不方便,换行处理麻烦. 基于在 ...

  10. 拒绝从入门到放弃_《Python 核心编程 (第二版)》必读目录

    目录 目录 关于这本书 必看知识点 最后 关于这本书 <Python 核心编程 (第二版)>是一本 Python 编程的入门书,分为 Python 核心(其实并不核心,应该叫基础) 和 高 ...