【HDU】6110 路径交(2017百度之星) 线段树+RMQ-LCA+树链的交
【题意】给定n个点的带边权树,m条编号1~m的路径,Q次询问编号区间[L,R]所有链的交集的长度。n<=500000。
【算法】线段树+RMQ-LCA+树链的交
【题解】树链的交:记一条链为(a1,b1),LCA为c1。另一条链为(a2,b2),LCA为c2。记a1a2,a1b2,b1a2,b1b2的LCA为d1,d2,d3,d4,按深度排序后得deep[d1]<=deep[d2]<=deep[d3]<=deep[d4],deep[c1]<=deep[c2]。
两条链有交集当且仅当deep[c1]<=deep[d1]&&deep[c2]<=deep[d3],此时树链(d3,d4)是两条链的交。(判断是否有交还有另一种方法:一条链的LCA在另一条链上时有交)
基于此,用线段树维护m条链的交集,回答询问。用RMQ-LCA预处理后可以降低复杂度(注意欧拉序开2倍)。
总复杂度O(n log n)。
#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
#define ll long long
using namespace std;
int read(){
char c;int s=,t=;
while(!isdigit(c=getchar()))if(c=='-')t=-;
do{s=s*+c-'';}while(isdigit(c=getchar()));
return s*t;
}
const int maxn=;
int dfn[maxn*],p[maxn],first[maxn],deep[maxn],Q[maxn*][],c[],d[],tot=,num=;
ll dis[maxn];
int logs[maxn*],n,m;
struct edge{int v,w,from;}e[maxn*];
struct cyc{int x,y;}O[maxn];
struct tree{int l,r;cyc o;}t[maxn*]; void insert(int u,int v,int w){tot++;e[tot].v=v;e[tot].w=w;e[tot].from=first[u];first[u]=tot;}
void dfs(int x,int fa){
dfn[++num]=x;
if(!p[x])p[x]=num;
for(int i=first[x];i;i=e[i].from)if(e[i].v!=fa){
deep[e[i].v]=deep[x]+;
dis[e[i].v]=dis[x]+e[i].w;
dfs(e[i].v,x);
dfn[++num]=x;
}
}
void prepare(){
logs[]=-;for(int i=;i<=num;i++)logs[i]=logs[i>>]+;
for(int i=;i<=num;i++)Q[i][]=dfn[i];
for(int j=;(<<j)<=num;j++){
for(int i=;i+(<<j)-<=num;i++){
if(deep[Q[i][j-]]<deep[Q[i+(<<(j-))][j-]])Q[i][j]=Q[i][j-];
else Q[i][j]=Q[i+(<<(j-))][j-];
}
}
}
int lca(int x,int y){
if(p[x]>p[y])swap(x,y);
int k=logs[p[y]-p[x]+];
return deep[Q[p[x]][k]]<deep[Q[p[y]-(<<k)+][k]]?Q[p[x]][k]:Q[p[y]-(<<k)+][k];
}
bool cmp(int a,int b){return deep[a]<deep[b];}
cyc merge(cyc o1,cyc o2){
if(!o1.x||!o2.x)return (cyc){,};
d[]=lca(o1.x,o2.x);d[]=lca(o1.x,o2.y);d[]=lca(o1.y,o2.x);d[]=lca(o1.y,o2.y);
c[]=lca(o1.x,o1.y);c[]=lca(o2.x,o2.y);
sort(d+,d++,cmp);if(deep[c[]]>deep[c[]])swap(c[],c[]);
if(deep[c[]]<=deep[d[]]&&deep[c[]]<=deep[d[]])return (cyc){d[],d[]};else return (cyc){,};//
}
void build(int k,int l,int r){
t[k].l=l;t[k].r=r;
if(l==r){t[k].o=O[l];return;}
int mid=(l+r)>>;
build(k<<,l,mid);build(k<<|,mid+,r);
t[k].o=merge(t[k<<].o,t[k<<|].o);
}
cyc query(int k,int l,int r){
if(l<=t[k].l&&t[k].r<=r)return t[k].o;
int mid=(t[k].l+t[k].r)>>;cyc x=(cyc){-,};
if(l<=mid)x=query(k<<,l,r);
if(r>mid){if(~x.x)x=merge(x,query(k<<|,l,r));else x=query(k<<|,l,r);}
return x;
}
int main(){
n=read();
for(int i=;i<n;i++){
int u=read(),v=read(),w=read();
insert(u,v,w);insert(v,u,w);
}
dfs(,);
prepare();
m=read();
for(int i=;i<=m;i++)O[i].x=read(),O[i].y=read();
build(,,m);
int T=read();
while(T--){
int l=read(),r=read();cyc x;
x=query(,l,r);
if(!x.x)printf("0\n");else
printf("%lld\n",dis[x.x]+dis[x.y]-*dis[lca(x.x,x.y)]);
}
return ;
}
这题WA之后对拍不出错,发现是生成的数据每个点的父亲编号严格小于它,导致深度关系和编号等价。
然后我在判断c1<=d1&&c2<=d3的时候比较了编号而非深度,于是GG。
【HDU】6110 路径交(2017百度之星) 线段树+RMQ-LCA+树链的交的更多相关文章
- 2017"百度之星"程序设计大赛 - 初赛(A) [ hdu 6108 小C的倍数问题 ] [ hdu 6109 数据分割 ] [ hdu 6110 路径交 ] [ hdu 6112 今夕何夕 ] [ hdu 6113 度度熊的01世界 ]
这套题体验极差. PROBLEM 1001 - 小C的倍数问题 题 OvO http://acm.hdu.edu.cn/showproblem.php?pid=6108 (2017"百度之星 ...
- HDU 6109 数据分割 【并查集+set】 (2017"百度之星"程序设计大赛 - 初赛(A))
数据分割 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submis ...
- HDU 6108 小C的倍数问题 【数学】 (2017"百度之星"程序设计大赛 - 初赛(A))
小C的倍数问题 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Sub ...
- HDU 6122 今夕何夕 【数学公式】 (2017"百度之星"程序设计大赛 - 初赛(A))
今夕何夕 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submis ...
- HDU 6113 度度熊的01世界 【DFS】(2017"百度之星"程序设计大赛 - 初赛(A))
度度熊的01世界 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Su ...
- 2017"百度之星"程序设计大赛 - 复赛1005&&HDU 6148 Valley Numer【数位dp】
Valley Numer Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Tota ...
- 2017"百度之星"程序设计大赛 - 复赛1003&&HDU 6146 Pokémon GO【数学,递推,dp】
Pokémon GO Time Limit: 3000/1500 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total ...
- 2017"百度之星"程序设计大赛 - 复赛1001&&HDU 6144 Arithmetic of Bomb【java大模拟】
Arithmetic of Bomb Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Other ...
- HDU 6118 度度熊的交易计划 【最小费用最大流】 (2017"百度之星"程序设计大赛 - 初赛(B))
度度熊的交易计划 Time Limit: 12000/6000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total S ...
随机推荐
- lintocde-247-线段树的查询 II
247-线段树的查询 II 对于一个数组,我们可以对其建立一棵 线段树, 每个结点存储一个额外的值 count 来代表这个结点所指代的数组区间内的元素个数. (数组中并不一定每个位置上都有元素) 实现 ...
- Java Map 在put值时value值不被覆盖
一.问题描述 最近在代码开发中遇到一个问题,在往Map中put文件路径值然后把Map放到List中去时,遇到问题是Map的后一个值总是把前一个值覆盖,导致最后Map中只有一个值. 二.解决办法(有如下 ...
- 使用nginx反向代理时,如何正确获取到用户的真实ip
在记录日志的的时候,获取用户的信息,比如用户的ip,浏览器等等信息是十分重要的. 但是在使用nginx反向代理的时候,可能经过转发无法获取到用户的真实的ip, 在此情况下需要配置nginx,让其在转发 ...
- Contest 6
A:容易发现这要求所有子集中元素的最高位1的位置相同,并且满足这个条件也是一定合法的.统计一下即可. #include<iostream> #include<cstdio> # ...
- oracle job定时执行存储过程
JOB定时跑插入语句1.建插入数据的存储过程create or replace procedure report_web asV_START_DATE DATE;V_END_DATE DATE;b ...
- 【转】NHibernate 各种数据库配置
转载自:http://terrylee.cnblogs.com/archive/2006/04/05/367381.html 自己只用了Oracle 可用! 三.常见的配置示例 Castle网站为我们 ...
- NOIP2016愤怒的小鸟 题解报告 【状压DP】
题目什么大家都清楚 题解 我们知道,三点确定一条抛物线,现在这条抛物线过原点,所以任意两只猪确定一条抛物线.通过运算的出对于两头猪(x1,y1),(x2,y2),他们所在抛物线a=(y1*x2-y2* ...
- [学习笔记]快速幂&&快速乘
本质:二进制拆分(你说倍增我也没脾气).然后是一个配凑. 合起来就是边二进制拆分,边配凑. 快速乘(其实龟速):把乘数二进制拆分.利用乘法分配率. 用途:防止爆long long 代码: ll qk( ...
- python 常用 time, datetime处理
python 中 time 有三种格式: float, struct tuple(time.struct_time 或 datetime.datetime), str 常用的: float --> ...
- Git之安装及使用
学习使用Git来管理平时自己写的demo代码和阅读的一些源码,因为一直在windows中操作所以开始学习用Git Bash操作在github上的代码.git命令和svn命令是很相似的,我觉得没有必要把 ...