数据结构(树链剖分):COGS 2109. [NOIP2015] 运输计划
2109. [NOIP2015] 运输计划
★★★ 输入文件:transport.in
输出文件:transport.out
简单对比
时间限制:1 s 内存限制:256 MB
【题目描述】
公元 2044 年,人类进入了宇宙纪元。
L 国有 n 个星球,还有 n-1 条双向航道,每条航道建立在两个星球之间,这 n-1 条航道连通了 L 国的所有星球。
小 P 掌管一家物流公司,该公司有很多个运输计划,每个运输计划形如:有一艘物
流飞船需要从 ui 号星球沿最快的宇航路径飞行到 vi 号星球去。显然,飞船驶过一条航道 是需要时间的,对于航道 j,任意飞船驶过它所花费的时间为 tj,并且任意两艘飞船之 间不会产生任何干扰。
为了鼓励科技创新,L 国国王同意小 P 的物流公司参与 L 国的航道建设,即允许小 P 把某一条航道改造成虫洞,飞船驶过虫洞不消耗时间。
在虫洞的建设完成前小 P 的物流公司就预接了 m 个运输计划。在虫洞建设完成后, 这 m 个运输计划会同时开始,所有飞船一起出发。当这 m 个运输计划都完成时,小 P 的 物流公司的阶段性工作就完成了。
如果小 P 可以自由选择将哪一条航道改造成虫洞,试求出小 P 的物流公司完成阶段 性工作所需要的最短时间是多少?
【输入格式】
第一行包括两个正整数 n、m,表示 L 国中星球的数量及小 P 公司预接的运输计划的数量,星球从 1 到 n 编号。
接下来 n-1 行描述航道的建设情况,其中第 i 行包含三个整数 ai, bi 和 ti,表示第i 条双向航道修建在 ai 与 bi 两个星球之间,任意飞船驶过它所花费的时间为 ti。 接下来 m 行描述运输计划的情况,其中第 j 行包含两个正整数 uj 和 vj,表示第 j 个运输计划是从 uj 号星球飞往 vj 号星球。
【输出格式】
共 1 行,包含 1 个整数,表示小 P 的物流公司完成阶段性工作所需要的最短时间。
【样例输入】
6 3
1 2 3
1 6 4
3 1 7
4 3 6
3 5 5
3 6
2 5
4 5
【样例输出】
11
【提示】
所有测试数据的范围和特点如下表所示
请注意常数因子带来的程序效率上的影响。
哎,拼命优化还是95分。
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int maxn=; int n,m;
int E[maxn][],cnt,fir[maxn],to[maxn<<],nxt[maxn<<],val[maxn<<];
void addedge(int a,int b,int v){
nxt[++cnt]=fir[a];to[cnt]=b;fir[a]=cnt;val[cnt]=v;
} int dep[maxn],son[maxn],sz[maxn],fa[maxn],dis[maxn]; void DFS(int x){
sz[x]=;
for(int i=fir[x];i;i=nxt[i])
if(to[i]!=fa[x]){
dep[to[i]]=dep[x]+;
dis[to[i]]=dis[x]+val[i];
fa[to[i]]=x;
DFS(to[i]);
sz[x]+=sz[to[i]];
if(sz[son[x]]<sz[to[i]])
son[x]=to[i];
}
} int top[maxn],ID[maxn],tot;
void DFS(int x,int tp){
ID[x]=++tot;top[x]=tp;
if(son[x])DFS(son[x],tp);
for(int i=fir[x];i;i=nxt[i])
if(!ID[to[i]])
DFS(to[i],to[i]);
} int Lca(int x,int y){
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]])swap(x,y);
x=fa[top[x]];
}
return dep[x]>dep[y]?y:x;
} int Mx1[maxn*],Mx2[maxn*];
struct Node{
int l,r;
Node(int L=,int R=){
l=L;r=R;
}
bool operator <(const Node &b)const{
return l<b.l;
}
}st[maxn]; int V[maxn];
void Update1(int x,int l,int r,int a,int b,int g){
if(l>=a&&r<=b){
Mx1[x]=max(Mx1[x],g);
return;
}
int mid=(l+r)>>;
if(mid>=a)Update1(x<<,l,mid,a,b,g);
if(mid<b)Update1(x<<|,mid+,r,a,b,g);
} void Update2(int x,int l,int r,int a,int b,int g){
if(l>=a&&r<=b){
Mx2[x]=max(Mx2[x],g);
return;
}
int mid=(l+r)>>;
if(mid>=a)Update2(x<<,l,mid,a,b,g);
if(mid<b)Update2(x<<|,mid+,r,a,b,g);
} int Query1(int x,int l,int r,int g){
if(l==r)
return Mx1[x];
int mid=(l+r)>>;
if(mid>=g)return max(Query1(x<<,l,mid,g),Mx1[x]);
else return max(Query1(x<<|,mid+,r,g),Mx1[x]);
} int Query2(int x,int l,int r,int g){
if(l==r)
return Mx2[x];
int mid=(l+r)>>;
if(mid>=g)return max(Query2(x<<,l,mid,g),Mx2[x]);
else return max(Query2(x<<|,mid+,r,g),Mx2[x]);
}
void Solve(int x,int y,int sum){
int tp=;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]])swap(x,y);
st[++tp]=Node(ID[top[x]],ID[x]);
x=fa[top[x]];
}
if(dep[x]<dep[y])swap(x,y);
if(x!=y)st[++tp]=Node(ID[son[y]],ID[x]);
sort(st,st+tp+);
for(int i=;i<=tp;i++)
Update1(,,n,st[i].l,st[i].r,sum); int L=;
for(int i=;i<=tp;i++){
if(L<=st[i].l-)
Update2(,,n,L,st[i].l-,sum);
L=st[i].r+;
}
if(L<=n)Update2(,,n,L,n,sum);
} int main(){
#ifndef ONLINE_JUDGE
freopen("transport.in","r",stdin);
freopen("transport.out","w",stdout);
#endif
scanf("%d%d",&n,&m);
for(int i=;i<n;i++){
scanf("%d%d%d",&E[i][],&E[i][],&E[i][]);
addedge(E[i][],E[i][],E[i][]);
addedge(E[i][],E[i][],E[i][]);
} DFS();
DFS(,); memset(Mx1,0x80,sizeof(Mx1));
memset(Mx2,0x80,sizeof(Mx2)); for(int i=;i<n;i++){
if(dep[E[i][]]<dep[E[i][]])
swap(E[i][],E[i][]);
V[E[i][]]=E[i][];
} int Max=,pa,pb;
while(m--){
int a,b,lca;
scanf("%d%d",&a,&b);
lca=Lca(a,b);
if(dis[a]+dis[b]-*dis[lca]>Max){
pa=a,pb=b;
Max=dis[a]+dis[b]-*dis[lca];
}
Solve(a,b,dis[a]+dis[b]-*dis[lca]);
} int ans=;
while(pa!=pb){
if(dep[pa]<dep[pb])swap(pa,pb);
int ret=max(Query1(,,n,ID[pa])-V[pa],Query2(,,n,ID[pa]));
ans=min(ans,ret);
pa=fa[pa];
} printf("%d\n",ans);
return ;
}
数据结构(树链剖分):COGS 2109. [NOIP2015] 运输计划的更多相关文章
- 数据结构(树链剖分):BZOJ 4034: [HAOI2015]T2
Description 有一棵点数为 N 的树,以点 1 为根,且树点有边权.然后有 M 个 操作,分为三种: 操作 1 :把某个节点 x 的点权增加 a . 操作 2 :把某个节点 x 为根的子树中 ...
- 数据结构--树链剖分准备之LCA
有关LCA的模板题 传送门 题目描述 如题,给定一棵有根多叉树,请求出指定两个点直接最近的公共祖先. 输入输出格式 输入格式: 第一行包含三个正整数N.M.S,分别表示树的结点个数.询问的个数和 ...
- 数据结构(树链剖分,堆):HNOI 2016 network
2215. [HNOI2016]网络 ★★★☆ 输入文件:network_tenderRun.in 输出文件:network_tenderRun.out 简单对比时间限制:2 s 内存 ...
- cogs 2109. [NOIP 2015] 运输计划 提高组Day2T3 树链剖分求LCA 二分答案 差分
2109. [NOIP 2015] 运输计划 ★★★☆ 输入文件:transport.in 输出文件:transport.out 简单对比时间限制:3 s 内存限制:256 MB [题 ...
- BZOJ_4326_[NOIP2015]_运输计划_(二分+LCA_树链剖分/Tarjan+差分)
描述 http://www.lydsy.com/JudgeOnline/problem.php?id=4326 给出一棵带有边权的树,以及一系列任务,任务是从树上的u点走到v点,代价为u到v路径上的权 ...
- 树链剖分-Hello!链剖-[NOIP2015]运输计划-[填坑]
This article is made by Jason-Cow.Welcome to reprint.But please post the writer's address. http://ww ...
- bzoj 4326: NOIP2015 运输计划(二分+树链剖分)
传送门 题解: 树链剖分快速求解任意两点间的路径的权值和: 然后,二分答案: 此题的难点是如何快速求解重合路径? 差分数组可以否??? 在此之前先介绍一下相关变量: int fa[maxn]; int ...
- JZYZOJ1454 NOIP2015 D2T3_运输计划 二分 差分数组 lca tarjan 树链剖分
http://172.20.6.3/Problem_Show.asp?id=1454 从这道题我充分认识到我的脑子里好多水orz. 如果知道了这个要用二分和差分写,就没什么思考上的难点了(屁咧你写了一 ...
- NOIP2015 运输计划 - 二分 + 树链剖分 / (倍增 + 差分)
BZOJ CodeVS Uoj 题目大意: 给一个n个点的边带权树,给定m条链,你可以选择树中的任意一条边,将它置为0,使得最长的链长最短. 题目分析: 最小化最大值,二分. 二分最短长度mid,将图 ...
随机推荐
- MD5 密码破解 碰撞 网站
MD5反向查询网站 http://www.cmd5.com/ 文件MD5值查询网站 http://www.atool.org/file_hash.php 个人对密码破解的理解 1.使用MD5对密码加密 ...
- Orchard路由随记(一)
对于Orchard来说,个人以为要真正理解Orchard,必须理解其路由工作方式. 一.Orchard的自定义路由由三种类型组成 1.分发类: HubRoute:其功能是按租户筛选出当前访问租户的路由 ...
- Struts2 中拦截器和Action的调用关系(写的很好)
http://blog.csdn.net/hackerain/article/details/6991082
- PropertyPlaceholderConfigurer的用法(使用spring提供的类读取数据库配置信息.properties)
http://www.cnblogs.com/wanggd/archive/2013/07/04/3172042.html(写的很好)
- CI 笔记5 (CI3.0 默认控制器,多目录)
在ci3.x中,不支持多级子目录的默认控制器设置, 解决方法如下: 在index.php中,添加 $routing['directory'] = 'admin';然后在默认的router.php的默 ...
- iOS菜鸟之FMDB的二次封装简单易用
闲来无事写点东西,希望大家多多指正! 大家先去git下载FMDB,然后将其中source文件夹中的fmdb文件夹拖入自己的项目中.最后就可以引用下面的代码对fmdb进行一次简单的封装. 这样可以更直观 ...
- linux防火墙解封某端口
首先,使用netstat –tunlp查看是否23端口被防火墙封掉了: 再使用iptables修改设置, # iptables -I INPUT -p tcp --dport 23 –jACCEPT ...
- PAT - 基础 - 龟兔赛跑
题目: 乌龟与兔子进行赛跑,跑场是一个矩型跑道,跑道边可以随地进行休息.乌龟每分钟可以前进3米,兔子每分钟前进9米:兔子嫌乌龟跑得慢,觉得肯定能跑赢乌龟,于是,每跑10分钟回头看一下乌龟,若发现自己超 ...
- php里 \r\n换行问题
<?php echo "hello"; echo "\r\n"; echo "world"; ?> 在浏览器输出的是hello ...
- DedeCMS时间格式
时间格式 {dede:field name='pubdate' function='strftime("%Y年%m月%d日 %H:%M:%S","@me")' ...