香港记者跑的比谁都快是众所周知的常识。

现在,香港记者站在一颗有 n 个点的树的根结点上(即1号点),编号为 i 的点拥有权值 a[i] ,数据保证每个点的编号都小于它任意孩子结点的别号。

我们假定这棵树的每个叶子结点都在发生一个大新闻,香港记者要用最少的耗时去报道其中的任意一个。

若香港记者目前处于第 i 号点上,那么它可以移动至以 i 为根的子树上的任意一点 j ,耗时 a[i]+(j−i)p ,p为给定常数。

请问这位香港记者搞哪个大新闻的耗时最短?所耗时间是多少?

解题报告:

用时:3h,2WA4TLE

这题开始写没写出来,然后想着优化\(O(n^2)\)的暴力,于是开了个以\(val[x]+f[x]\)为判断条件的单调栈,然后还是TLE65

最后写的决策单调性:

对于一个点x,存在一个管辖的范围,在这个范围内x的转移是最优的,然后我们就维护这个东西,每次新加入一个点,我们拿这个点和栈顶的点分别和管辖的那个点计算然后作比较,维护单调性,然后每一次转移,我们二分到第一个能管辖到x的点就是最优转移,求一个点能管辖的范围,那么就二分一个\((x'-x)^p\)中的\(x'\),直到它刚好比栈顶元素优为止,其后就是它能管辖的范围

还有对于树上的栈维护,还要消除子树间的影响,我们进入一颗子树前,还要把在其他子树弹出的元素都放回栈中.

复杂度\(O(nlogn)\)

#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <queue>
#include <cstdio>
#include <cmath>
#define RG register
#define il inline
#define iter iterator
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
using namespace std;
const int N=1e5+5;
typedef long long ll;
const ll inf=1e18;
int head[N],to[N<<1],n,nxt[N<<1],num=0,val[N],p;
void link(int x,int y){nxt[++num]=head[x];to[num]=y;head[x]=num;}
ll ans=inf,f[N],mul[N];
ll qm(ll x,ll k){
ll sum=1;
while(k){
if(k&1)sum*=x;
x*=x;k>>=1;
}
return sum;
}
bool vis[N];
struct node{
int x,last,fo;
node(){}
node(int _x,int _last,int _fo){x=_x;last=_last;fo=_fo;}
}q[N],st[N];
int top=0,stop=0,du[N];
void dfs(int x,int last){
int u;
if(x!=1){
int l=1,r=top,mid,ret=l,last;
while(l<=r){
mid=(l+r)>>1;
if(q[mid].last<=x)ret=mid,l=mid+1;
else r=mid-1;
}
f[x]=f[q[ret].x]+val[q[ret].x]+mul[x-q[ret].x];
while(top && q[top].last>=x){
u=q[top].x;last=q[top].last;
if(f[x]+val[x]+mul[last-x]<=f[u]+val[u]+mul[last-u])
st[++stop]=q[top--],st[stop].fo=x;
else break;
}
l=Max(q[top].last+1,x+1);r=n;ret=0;u=q[top].x;
while(l<=r){
mid=(l+r)>>1;
if(f[x]+val[x]+mul[mid-x]<=f[u]+val[u]+mul[mid-u])
ret=mid,r=mid-1;
else l=mid+1;
}
if(ret)q[++top]=node(x,ret,0);
}
else f[x]=0,q[++top]=node(1,1,0);
if(!du[x])ans=Min(ans,f[x]);
for(int i=head[x];i;i=nxt[i]){
u=to[i];if(u==last)continue;
dfs(u,x);
}
if(q[top].x==x)top--;
while(stop && st[stop].fo==x)q[++top]=st[stop--];
}
void work()
{
int x;
scanf("%d%d",&n,&p);
for(int i=1;i<=n;i++){
scanf("%d%d",&val[i],&x);
if(x)link(x,i),link(i,x),du[x]++;
}
for(int i=0;i<=n;i++){
mul[i]=qm(i,p);
if(mul[i]>=inf){
for(int j=i;j<=n;j++)
mul[j]=inf;
break;
}
}
dfs(1,1);
printf("%lld\n",ans);
} int main()
{
work();
return 0;
}

51Nod 1781 跑的比谁都快的更多相关文章

  1. 跑的比谁都快 51Nod - 1789

    香港记者跑的比谁都快是众所周知的常识.   现在,香港记者站在一颗有  nn 个点的树的根结点上(即1号点),编号为  ii 的点拥有权值  a[i]a[i] ,数据保证每个点的编号都小于它任意孩子结 ...

  2. BZOJ:4219: 跑得比谁都快 3007: 拯救小云公主

    4219: 跑得比谁都快 3007: 拯救小云公主 三角剖分的解释可以看这里:http://www.cnblogs.com/Enceladus/p/6706444.html 后一道是前一道的弱化版. ...

  3. [51nod1789] 跑得比谁都快

    题面 题解 设\(f[i]\)为根节点到\(i\)的最小耗时 设\(S\)为\(i\)的祖先集合, 可以得到 \[ f[i] = min(f[j] + (i - j)^p),j \in S \] 对于 ...

  4. iOS 1 到 iOS 10 ,我都快老了

    iOS 1:iPhone诞生 虽然很难想像,但初代iPhone在问世时在功能方面其实远远落后于那时的竞争对手,比如Windows Mobile.Palm OS.塞班.甚至是黑莓.它不支持3G.多任务. ...

  5. jdk1.8新特性,还不知道的朋友还不看看,1.9都快出来了

    一.接口的默认方法 Java 8允许我们给接口添加一个非抽象的方法实现,只需要使用 default关键字即可,这个特征又叫做扩展方法,示例如下:代码如下:interface Formula {     ...

  6. 入坑IT都快十年了

    一起帮的开发直播已经告一段落:一是主体的功能差不多都实现了,二是用到的架构技术都展示得差不多了.以后就算继续开发,也应该都是一些“技术上”重复的工作而已.整个直播过程耗时近半年,SVN提交1062次, ...

  7. 今天我看了一个H5游戏EUI的例子,我都快分不清我到底是在用什么语言编译了代码了,作为刚刚学习H5游戏开发的菜鸟只能默默的收集知识

    今天看了一个EUI的demo,也是接触H5游戏开发的第五天了,我想看看我能不能做点什么出来,哎,自己写果然还是有问题的.在看EUI哪一个demo的时候就遇见了一些摇摆不定的问题,我觉得提出来 1.to ...

  8. Java9都快发布了,Java8的十大新特性你了解多少呢?

    Java 9预计将于今年9月份发布,这是否会是一次里程碑式的版本,我们拭目以待.今天,我们先来复习一下2014年发布的Java 8的十大新特性. Java 8可谓是自Java 5以来最具革命性的版本了 ...

  9. 51nod 1781 Pinball(线段树)

    题面 Pinball的游戏界面由m+2行.n列组成.第一行在顶端.一个球会从第一行的某一列出发,开始垂直下落,界面上有一些漏斗,一共有m个漏斗分别放在第2~m+1行,第i个漏斗的作用是把经过第i+1行 ...

随机推荐

  1. 第一周-JAVA基本概念

    1. 本周学习总结 本周学习内容: 1.JAVA的发展 2.JDK,JVM,JRE, 3.掌握JAVA的组成结构 4.掌握使用简单的编译器写javac与java命令, 关键概念之间的联系: JVM:将 ...

  2. Struts2之配置文件中Action的详细配置

    在Struts2之配置一文中,我们知道一个struts配置文件可以分为三部分:常量配置    包含其他配置文件的配置    Action配置  . 这其中 常量配置  和 包含其他配置文件的配置  二 ...

  3. linux下面的智能解压脚本smart解压

    file命令可以辨别出一个给定文件的文件类型,如:file lf.gz,其输出结果为: lf.gz: gzip compressed data, deflated, original filename ...

  4. Python内置函数(61)——eval

    英文文档: eval(expression, globals=None, locals=None) The arguments are a string and optional globals an ...

  5. You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'like '%逸%'' at line 1

    <?xml version="1.0" encoding="UTF-8"?> <!DOCTYPE mapper PUBLIC "-/ ...

  6. 其他—cooki和session

    cookie Cookie的由来 大家都知道HTTP协议是无状态的. 无状态的意思是每次请求都是独立的,它的执行情况和结果与前面的请求和之后的请求都无直接关系,它不会受前面的请求响应情况直接影响,也不 ...

  7. 新概念英语(1-9)How is Ema?

    A:Hello Helen. B:Hi Steven. A:How are you today? B:I'm very well, thank you. And you? A:I'm fine tha ...

  8. SQL Server 利用触发器对多表视图进行更新

    其步骤就是:利用update操作触发器产生的2个虚拟表[inserted]用来存储修改的数据信息和[deleted]表,然后将对应的数据更新到对应数据表中的字段信息中: 1.首先创建3个表: a.信息 ...

  9. Linux系统的基本使用

    曾经在网上看到一个一篇文章,说到了Linux学习的入门与学习技巧,也就是:初学者可以自己安装虚拟机,然后把linux常用命令例如cd.ls.chmod.useradd.vi等等多练习几十遍,把自己敲打 ...

  10. python列表基础操作

    Python列表基本操作 记住一句话,叫做顾首不顾尾 首先我们来定义一个列表 name = ["jixuege","dajiba","boduoye& ...