●BZOJ 3551 [ONTAK2010]Peaks(在线)
题链:
http://www.lydsy.com/JudgeOnline/problem.php?id=3551
题解:
最小生成树 Kruskal,主席树,在线
这个做法挺巧妙的。。。
以Kruskal算法为基础,如果在用边 e(u,v,w) 合并 u 和 v 所在的联通块时,
我们新加一个节点 x(同时给它一个权值 w,即边 e 的权值),
使得 u的联通块和 v的联通块通过这个节点 x 来合并为一个联通块。
那么当Kruskal算法完成时,那么也就生成了一颗二叉树。
不难发现,这个二叉树非常棒啊:
1).叶子代表原图中的点,共有 N 个叶子
2).对于一个节点 x 来说,如果其点权为 w,
则表明 x 的子树的叶子节点所代表的的那些原图中的点可以通过边权不超过 w 的边互相通达。
所以按照如上方法生成了一颗二叉树后,
dfs一遍计算出每个节点所包含叶子节点的范围
(l[u],r[u]表示 u 这个节点的所包含叶子节点的范围是第 l[u]个叶子到第 R[u]个叶子)
然后按照叶子节点的顺序对原图的点的高度建立主席树。
每次查询v,x,k时,
就在二叉树中从该叶子向上倍增,
找到一个最大的子树使得该子树的根的权值不超过 x
然后得到范围 l[x] r[x],并在主席树 rt[l[x]-1]~rt[r[x]] 中查询第 k 大就好了。
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 105000
using namespace std;
int H[MAXN],tmp[MAXN],Ord[MAXN];
int N,M,Q,tnt,Ont;
struct Edge{
int u,v,w;
bool operator <(const Edge &rtm)const{
return w<rtm.w;
}
}E[MAXN*5];
struct CMT{
int rt[MAXN],ls[MAXN*20],rs[MAXN*20],cnt[MAXN*20],sz;
void Insert(int &u,int l,int r,int p){
sz++; cnt[sz]=cnt[u]; ls[sz]=ls[u]; rs[sz]=rs[u];
u=sz; cnt[u]++;
if(l==r) return;
int mid=(l+r)>>1;
if(p<=mid) Insert(ls[u],l,mid,p);
else Insert(rs[u],mid+1,r,p);
}
int Query(int v,int u,int l,int r,int K){
if(l==r) return l;
int mid=(l+r)>>1,rcnt=cnt[rs[u]]-cnt[rs[v]];
if(K<=rcnt) return Query(rs[v],rs[u],mid+1,r,K);
else return Query(ls[v],ls[u],l,mid,K-rcnt);
}
void Build(){
for(int i=1;i<=N;i++){
rt[i]=rt[i-1];
Insert(rt[i],1,tnt,H[Ord[i]]);
}
}
}DT;//维护成重构的二叉树的底层叶子节点信息的主席树
struct BT{
int bel[MAXN*2],ls[MAXN*2],rs[MAXN*2],l[MAXN*2],r[MAXN*2],val[MAXN*2],fa[MAXN*2][20];
int sz,rt;
void Newnode(int u,int lson,int rson,int w){
val[u]=w; bel[lson]=u; bel[rson]=u;
ls[u]=lson; rs[u]=rson;
}
int Find(int x){
return x==bel[x]?x:bel[x]=Find(bel[x]);
}
void Dfs(int u,int f){
fa[u][0]=f;
for(int k=1;k<20;k++){
if(!fa[u][k-1]||!fa[fa[u][k-1]][k-1]) break;
fa[u][k]=fa[fa[u][k-1]][k-1];
}
if(!ls[u]&&!rs[u]){
Ord[++Ont]=u; l[u]=r[u]=Ont;
return;
}
Dfs(ls[u],u); Dfs(rs[u],u);
l[u]=l[ls[u]]; r[u]=r[rs[u]];
}
void Build(){
sz=N;
for(int i=1;i<=2*N;i++) bel[i]=i;
for(int i=1,fu,fv;i<=M;i++){
fu=Find(E[i].u);
fv=Find(E[i].v);
if(fu==fv) continue;
Newnode(++sz,fu,fv,E[i].w);
}
rt=sz; Dfs(rt,0);
}
int Query(int u,int w,int K){
for(int k=19;k>=0;k--){
if(!fa[u][k]||val[fa[u][k]]>w) continue;
u=fa[u][k];
}
int lrt=DT.rt[l[u]-1],rrt=DT.rt[r[u]];
if(DT.cnt[rrt]-DT.cnt[lrt]<K) return -1;
return tmp[DT.Query(lrt,rrt,1,tnt,K)];
}
}GT;//根据Kruskal重构的二叉树
void read(int &x){
static int f; static char ch;
x=0; f=1; ch=getchar();
while(ch<'0'||'9'<ch){if(ch=='-') f=-1; ch=getchar();}
while('0'<=ch&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}
x=x*f;
}
void readin(){
read(N); read(M); read(Q);
for(int i=1;i<=N;i++) read(H[i]),tmp[i]=H[i];
sort(tmp+1,tmp+N+1);
tnt=unique(tmp+1,tmp+N+1)-tmp-1;
for(int i=1;i<=N;i++)
H[i]=lower_bound(tmp+1,tmp+tnt+1,H[i])-tmp;
for(int i=1;i<=M;i++)
read(E[i].u),read(E[i].v),read(E[i].w);
sort(E+1,E+M+1);
}
void answer(){
static int V,X,K,lastANS;
for(int i=1;i<=Q;i++){
scanf("%d%d%d",&V,&X,&K);
if(lastANS!=-1) V^=lastANS,X^=lastANS,K^=lastANS;
lastANS=GT.Query(V,X,K);
printf("%d\n",lastANS);
}
}
int main(){
readin();
GT.Build();
DT.Build();
answer();
return 0;
}
●BZOJ 3551 [ONTAK2010]Peaks(在线)的更多相关文章
- BZOJ 3551: [ONTAK2010]Peaks加强版 [Kruskal重构树 dfs序 主席树]
3551: [ONTAK2010]Peaks加强版 题意:带权图,多组询问与一个点通过边权\(\le lim\)的边连通的点中点权k大值,强制在线 PoPoQQQ大爷题解传送门 说一下感受: 容易发现 ...
- bzoj 3551: [ONTAK2010]Peaks加强版
Description [题目描述]同3545 Input 第一行三个数N,M,Q. 第二行N个数,第i个数为h_i 接下来M行,每行3个数a b c,表示从a到b有一条困难值为c的双向路径. 接下来 ...
- BZOJ.3551.[ONTAK2010]Peaks加强版(Kruskal重构树 主席树)
题目链接 \(Description\) 有n个座山,其高度为hi.有m条带权双向边连接某些山.多次询问,每次询问从v出发 只经过边权<=x的边 所能到达的山中,第K高的是多少. 强制在线. \ ...
- bzoj 3551 [ONTAK2010]Peaks加强版(kruskal,主席树,dfs序)
Description [题目描述]同3545 Input 第一行三个数N,M,Q. 第二行N个数,第i个数为h_i 接下来M行,每行3个数a b c,表示从a到b有一条困难值为c的双向路径. 接下来 ...
- 【刷题】BZOJ 3551 [ONTAK2010]Peaks加强版
Description [题目描述]同3545 Input 第一行三个数N,M,Q. 第二行N个数,第i个数为h_i 接下来M行,每行3个数a b c,表示从a到b有一条困难值为c的双向路径. 接下来 ...
- BZOJ 3551: [ONTAK2010]Peaks加强版 Kruskal重构树+dfs序+主席树+倍增
建出来 $Kruskal$ 重构树. 将询问点向上跳到深度最小,且合法的节点上. 那么,得益于重构树优美的性质,这个最终跳到的点为根的所有子节点都可以与询问点互达. 对于子树中求点权第 $k$ 大的问 ...
- 3551: [ONTAK2010]Peaks加强版
3551: [ONTAK2010]Peaks加强版 https://www.lydsy.com/JudgeOnline/problem.php?id=3551 分析: kruskal重构树 + 倍增 ...
- bzoj 3545&&3551: [ONTAK2010]Peaks &&加强版 平衡树&&并查集合并树&&主席树
3545: [ONTAK2010]Peaks Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 635 Solved: 177[Submit][Stat ...
- bzoj 3545/3551: [ONTAK2010]Peaks -- 主席树,最小生成树,倍增
3545: [ONTAK2010]Peaks Time Limit: 10 Sec Memory Limit: 128 MB Description 在Bytemountains有N座山峰,每座山峰 ...
随机推荐
- itchat 微信的使用
#coding=utf8 import itchat # 自动回复 # 封装好的装饰器,当接收到的消息是Text,即文字消息 @itchat.msg_register('Text') def text ...
- 关于tomcat部署应用的三种方式
关于tomcat部署应用虽然不是一个经常的操作,因为一旦选择了一种部署方式,我们其他的应用就会不经大脑的使用这种既定模式, 如果不使用这种部署方式,但是对于其他的部署方式不是很清楚的话,很容易抓瞎,所 ...
- var 和 let 的异同?
相同点 声明后未赋值表现一致 不同点 1.使用未声明的变量表现不同 2.变量作用范围不同 3.var可以声明多次 let只能声明一次 let的好处就是当我们在写代码的时候可以避免在不知道的情况下重复声 ...
- 我的PCB电路设计(一)
我的制板规则 过孔大小:14/24mil-(12/22-28/50) 一般过孔没必要太大.如果电流较大可以适当增大过孔,或者多加几个过孔 线宽大小:小信号线8mil,大电流线不等按1A电流30mil ...
- 使用Python3爬虫抓取网页来下载小说
很多时候想看小说但是在网页上找不到资源,即使找到了资源也没有提供下载,小说当然是下载下来用手机看才爽快啦! 于是程序员的思维出来了,不能下载我就直接用爬虫把各个章节爬下来,存入一个txt文件中,这样, ...
- istio入门(03)istio的helloworld-场景说明
一.原生应用 四个微服务: python微服务:一个deployment(deployment含有一个pod,pod内含有一个容器) java微服务:三个deployment(deployment含有 ...
- 新概念英语(1-1)Excuse me!
Excuse me!Whose handbag is it? A:Excuse me! B:Yes? A:Is this your handbag? B:Pardon? A:Is this your ...
- Python学习之再议row_input
再议raw_input birth = raw_input('birth: ') if birth < 2000: print '00前' else: print '00后' 运行结果: bir ...
- Mysql 5.1的坑
1.数据库表是区分大小写的 之前程序在5.7数据库没问题,测试环境上数据库是5.1的,就提示表找不到. 2.同样的sql,在5.1上会提示事务获取锁失败,超时返回.而5.7上正常.原因暂未找到.
- 初学Java Web(4)——Servlet学习总结
经过一段时间的学习,对于Servlet有了新的不一样的见解,在这里做一下总结,将近来学习到的知识总结一下. Servlet 的请求流程 浏览器发出请求:http://localhost:80/xxx1 ...