$NOI2014$ 购票(斜率优化 点分治)
\(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]);
}
随机推荐
- 20180803-Java 流(Stream)、文件(File)和IO
Java 流(Stream).文件(File)和IO 下面的程序示范了用read()方法从控制台不断读取字符直到用户输入"q". // 使用BufferedReader 在控制台读 ...
- 调试 & 常数优化:我有特别的 debug 技巧
rxz 的调试技巧(https://www.zhihu.com/question/60719584/answer/179363450): #define DEBUG printf("Pass ...
- linux用setup命令来更改ip配置
在有安装系统桌面情况下,可以使用图形化形式来配置ip地址, 在命令行下,输入“setup”调出网卡.防火墙等配置界面: 2 选择“network configuration“,回车: 选择“devic ...
- MySQL innodb的组合索引各个列中的长度不能超过767,
MySQL索引的索引长度问题 MySQL的每个单表中所创建的索引长度是有限制的,且对不同存储引擎下的表有不同的限制. 在MyISAM表中,创建组合索引时,创建的索引长度不能超过1000,注意这里索 ...
- mysql_DML_select_聚合join
聚合函数: select avg(salary)//平均值 from wsb; select sum(salary)//总和 from wsb; select max(salary)//最大 from ...
- jQuery插件3种类型
1.封装对象方法的插件 这种插件是将对象方法封装起来,用于对通过选择器获取的jQuery对象进行操作,是最常见的一种插件. 此类插件可以发挥出jQuery选择器的强大优势,有相当一部分的jQuery方 ...
- 腾讯两大开源项目Tars、TSeer
6月25日,在LC3(LinuxCon + ContainerCon + CloudOpen)中国2018大会上,腾讯宣布其两大开源项目——RPC开发框架Tars.轻量化名字服务方案TSeer,加入L ...
- SEC3 - MySQL常见命令
1.查看当前所有的数据库 show databases; 2. 打开指定的库名 use 库名称: 3.查看当前库中所有的表 show tables; 4. 查看其他库的所有表 show tables ...
- 边界安全 - CDN/DMZ/网络协议
CDN 工具 - LuManager CDN DMZ 网络协议 - DNS Win7下搭建DNS服务器 - BIND 根域 顶级域(即相关国家域名管理机构的数据库,如中国的CNNIC) com n ...
- mysql使用localhost可以访问,使用ip地址无法访问
本地安装的mysql服务,使用localhost可以连接,ip地址发无法连接: 解决办法: 进入mysql命令界面,输入select host,user from mysql.user; host字段 ...