https://www.luogu.org/problemnew/show/P3994

设dp[i] 表示第i个城市到根节点的最小花费

dp[i]=min{ (dis[i]-dis[j])*P[i]+Q[i]+dp[j] }

这是O(n^2)的

这个式子可以斜率优化

dp[i]+dis[j]*P[i]=dis[i]*P[i]+Q[i]+dp[j]

就是一条斜率为P[i]的直线,截(dis[j],dp[j])的最小截距

在根往下走的过程中,斜率单调递增

这就体现了 为什么题目中说“i号城市是j号城市的某个祖先,那么一定存在Pi<=Pj”

我们按dfs序dp

现在唯一的问题就是如何得到 一个点到根节点路径上的单调队列

只需要考虑如何去除兄弟节点的子树对单调队列的影响

即在一个节点退出dfs时,将单调队列恢复为这个节点开始dfs的情况

头指针只是不断的+1,没有涉及到单调队列中元素的修改,所以记录下头指针在哪个位置即可

尾指针涉及到元素的替换,但是它只会替换一个元素,所以记录下尾指针的位置,以及被当前点替换的元素是谁

当节点退出dfs时,恢复记录的这三个值即可

这样的话,一个节点多次出队入队,时间复杂度就不是O(n)了

所以二分出队位置,时间复杂度为O(nlogn)

朴素的DP:

#include<cstdio>
#include<iostream>
#include<algorithm> using namespace std; #define N 1000001 typedef long long LL; int P[N],Q[N]; int front[N],to[N<<],nxt[N<<],val[N<<],tot; int fa[N]; LL dis[N]; int t;
LL mi[N]; void read(int &x)
{
x=; char c=getchar();
while(!isdigit(c)) c=getchar();
while(isdigit(c)) { x=x*+c-''; c=getchar(); }
} void add(int u,int v,int w)
{
to[++tot]=v; nxt[tot]=front[u]; front[u]=tot; val[tot]=w;
to[++tot]=u; nxt[tot]=front[v]; front[v]=tot; val[tot]=w;
} void dfs(int x,int f)
{
for(int i=front[x];i;i=nxt[i])
{
if(to[i]==f) continue;
dis[to[i]]=dis[x]+val[i];
mi[to[i]]=dis[to[i]]*P[to[i]]+Q[to[i]];
t=fa[to[i]];
while(t!=)
{
mi[to[i]]=min(mi[to[i]],(dis[to[i]]-dis[t])*P[to[i]]+Q[to[i]]+mi[t]);
t=fa[t];
}
dfs(to[i],x);
}
} int main()
{
int n,s;
read(n);
for(int i=;i<n;++i)
{
read(fa[i+]); read(s); read(P[i+]); read(Q[i+]);
add(fa[i+],i+,s);
}
dfs(,);
for(int i=;i<=n;++i) cout<<mi[i]<<'\n';
}

斜率优化,暴力出队:

#include<cstdio>
#include<iostream> using namespace std; #define N 1000001 typedef long long LL; int front[N],nxt[N],to[N],tot,val[N]; int P[N],Q[N]; int q[N],head,tail; LL dis[N];
LL dp[N]; void read(int &x)
{
x=; char c=getchar();
while(!isdigit(c)) c=getchar();
while(isdigit(c)) { x=x*+c-''; c=getchar(); }
} void add(int u,int v,int w)
{
to[++tot]=v; nxt[tot]=front[u]; front[u]=tot; val[tot]=w;
} inline double X(int i,int j) { return dis[j]-dis[i]; }
inline double Y(int i,int j) { return dp[j]-dp[i]; } void dfs(int x)
{
int now_h=head,now_t=tail;
while(head<tail- && Y(q[head],q[head+])<P[x]*X(q[head],q[head+])) head++;
int j=q[head];
dp[x]=(dis[x]-dis[j])*P[x]+dp[j]+Q[x];
while(head<tail- && Y(q[tail-],q[tail-])*X(q[tail-],x)>X(q[tail-],q[tail-])*Y(q[tail-],x)) tail--;
int rr=q[tail];
q[tail++]=x;
for(int i=front[x];i;i=nxt[i])
dis[to[i]]=dis[x]+val[i],dfs(to[i]);
head=now_h; q[tail-]=rr; tail=now_t;
} int main()
{
int n;
read(n);
int fa,d;
for(int i=;i<=n;++i)
{
read(fa); read(d);
add(fa,i,d);
read(P[i]); read(Q[i]);
}
for(int i=front[];i;i=nxt[i])
{
dis[to[i]]=val[i];
q[head=]=; tail=;
dfs(to[i]);
}
for(int i=;i<=n;++i) cout<<dp[i]<<'\n';
}

斜率优化,二分出队

#include<cstdio>
#include<iostream> using namespace std; #define N 1000001 typedef long long LL; int front[N],nxt[N],to[N],tot,val[N]; int P[N],Q[N]; int q[N],head,tail; LL dis[N];
LL dp[N]; void read(int &x)
{
x=; char c=getchar();
while(!isdigit(c)) c=getchar();
while(isdigit(c)) { x=x*+c-''; c=getchar(); }
} void add(int u,int v,int w)
{
to[++tot]=v; nxt[tot]=front[u]; front[u]=tot; val[tot]=w;
} inline double X(int i,int j) { return dis[j]-dis[i]; }
inline double Y(int i,int j) { return dp[j]-dp[i]; } void dfs(int x)
{
int now_h=head,now_t=tail;
int l=head,r=tail-,mid,tmp=-;
while(l<=r)
{
mid=l+r>>;
if(Y(q[mid],q[mid+])>=P[x]*X(q[mid],q[mid+])) tmp=mid,r=mid-;
else l=mid+;
}
if(tmp!=-) head=tmp;
else head=tail-;
int j=q[head];
dp[x]=(dis[x]-dis[j])*P[x]+dp[j]+Q[x];
l=head,r=tail-,tmp=-;
while(l<=r)
{
mid=l+r>>;
if(Y(q[mid],q[mid+])*X(q[mid+],x)<=X(q[mid],q[mid+])*Y(q[mid+],x)) tmp=mid,l=mid+;
else r=mid-;
}
if(tmp!=-) tail=tmp+;
else tail=head+;
int rr=q[tail];
q[tail++]=x;
for(int i=front[x];i;i=nxt[i])
dis[to[i]]=dis[x]+val[i],dfs(to[i]);
head=now_h; q[tail-]=rr; tail=now_t;
} int main()
{
int n;
read(n);
int fa,d;
for(int i=;i<=n;++i)
{
read(fa); read(d);
add(fa,i,d);
read(P[i]); read(Q[i]);
}
for(int i=front[];i;i=nxt[i])
{
dis[to[i]]=val[i];
q[head=]=; tail=;
dfs(to[i]);
}
for(int i=;i<=n;++i) cout<<dp[i]<<'\n';
}

洛谷 P3994 高速公路的更多相关文章

  1. 洛谷 P3994 高速公路(斜率优化)

    题目链接 题意:给出一棵树,\(1\) 号点为根,边上有边权. 每个点有两个参数 \(p_i,q_i\) 如果你想从 \(i\) 号点到与其距离为 \(d\) 的 \(j\) 号点,那么你需花费 \( ...

  2. 【洛谷p3994】Highway 二分+斜率优化DP

    题目大意:给你一颗$n$个点的有根树,相邻两个点之间有距离,我们可以从$x$乘车到$x$的祖先,费用为$dis\times P[x]+Q[x]$,问你除根以外每个点到根的最小花费. 数据范围:$n≤1 ...

  3. 洛谷P3994 Highway(树形DP+斜率优化+可持久化线段树/二分)

    有点类似NOI2014购票 首先有方程$f(i)=min\{f(j)+(dep_i-dep_j)*p_i+q_i\}$ 这个显然是可以斜率优化的... $\frac {f(j)-f(k)}{dep_j ...

  4. 洛谷P2221 高速公路【线段树】

    题目:https://www.luogu.org/problemnew/show/P2221 题意:有n个节点排成一条链,相邻节点之间有一条路. C u v val表示从u到v的路径上的每条边权值都加 ...

  5. 洛谷P2242 公路维修问题

    To 洛谷.2242 公路维修问题 题目描述 由于长期没有得到维修,A国的高速公路上出现了N个坑.为了尽快填补好这N个坑,A国决定对M处地段采取交通管制.为了求解方便,假设A国的高速公路只有一条,而且 ...

  6. [洛谷P3761] [TJOI2017]城市

    洛谷题目链接:[TJOI2017]城市 题目描述 从加里敦大学城市规划专业毕业的小明来到了一个地区城市规划局工作.这个地区一共有ri座城市,<-1条高速公路,保证了任意两运城市之间都可以通过高速 ...

  7. 洛谷 P5638 光骓者的荣耀

    洛谷 P5638 [CSGRound2]光骓者的荣耀 洛谷传送门 题目背景 小 K 又在做白日梦了.他进入到他的幻想中,发现他打下了一片江山. 题目描述 小 K 打下的江山一共有nn个城市,城市ii和 ...

  8. 洛谷 P6383 -『MdOI R2』Resurrection(DP)

    洛谷题面传送门 高速公路上正是补 blog 的时候,难道不是吗/doge,难不成逆在高速公路上写题/jy 首先形成的图显然是连通图并且有 \(n-1\) 条边.故形成的图是一棵树. 我们考虑什么样的树 ...

  9. 洛谷1640 bzoj1854游戏 匈牙利就是又短又快

    bzoj炸了,靠离线版题目做了两道(过过样例什么的还是轻松的)但是交不了,正巧洛谷有个"大牛分站",就转回洛谷做题了 水题先行,一道傻逼匈牙利 其实本来的思路是搜索然后发现写出来类 ...

随机推荐

  1. 重置Oracle配置

    经常被ORACLE坑,作为一个只需要开发时候连连ORACLE的程序员,在经历了一次又一次的折腾之后,决定还是把这些琐碎的事情写下来. 经常在虚拟机中使用ORACLE,ORACLE的网络配置有一些变化就 ...

  2. jira webhook 事件触发并程序代码调用jenkins接口触发构建操作

    要解决的问题 开发管理工具触发站点构建事件,事件处理中需要调用Jenkins接口开始构建动作. 我的应用场景: 使用jira作为管理工具,在jira中创建自定义的工作流来规定测试,上线,发布等流程,并 ...

  3. [算法总结] 6 道题搞定 BAT 面试——堆栈和队列

    本文首发于我的个人博客:尾尾部落 0. 基础概念 栈:后进先出(LIFO) 队列:先进先出(FIFO) 1. 栈的 java 实现 import java.util.Arrays; public cl ...

  4. 【容器魔方解读】AWS Re:Invent 2018大会

    每年云计算领域技术与商业风向标之一的AWS Re:Invent大会上周在美国拉斯维加斯召开,如往届一样,AWS密集发布了上百项的新产品或新技术.随着国内近两年云计算尤其是公有云的普及度越来越高,国内各 ...

  5. 企业落地Kubernetes的问题与对策

    在当今云计算领域,“容器技术”已经从三四年前的炒作期正式进入了产业落地期,而Kubernetes作为容器平台的标准已经得到了广泛应用. Kubernetes从2014年6月由Google宣布开源,到2 ...

  6. sql注入语句整理

    1.判断有无注入点; and 1=1 and 1=2 2.猜表一般的表的名称无非是admin adminuser user pass password 等..and 0<>(select ...

  7. 每天学一点---document.createDocumentFragment

    document.createDocumentFragment  用于创建文档对象,创建好的对象存在于内存中(不会引起回流,对元素位置和几何上的运算),不是附着在DOM树上,所以有更好的性能 可将该文 ...

  8. Python Pandas read_csv报错

    为实现文本去重(将前面采集的数据进行两两对比删除重复),写了以下代码. #-*- coding: utf-8 -*-import pandas as pd inputfile = 'e:/data/H ...

  9. Web网络服务介绍

    Web网络服务也叫WWW(World Wide Web),一般是指能够让用户通过浏览器访问到互联网中文档资源服务.目前提供WEB网络服务的程序有Apache .Nginx 和  IIS  等等,Web ...

  10. C#简述(三)

    详细请参考:http://www.runoob.com/csharp/csharp-string.html 1.C# 字符串(String) 在 C# 中,可以使用字符数组来表示字符串,但是,更常见的 ...