题意

给出一个图,边有边权,点有点权,每次询问一个点 \(x\) 只走边权小于等于 \(d\) 的边能到达的点中点权第 \(k\) 大。

强制在线,\(n\le 10^5,m,q\le 5\times 10^5\)

分析

如果可以离线的话,我们可以用一个并查集(路径压缩)维护连通性,在并查集的每个点上存一个权值线段树,每次merge的时候就把下面的线段树合并到上面。做kruskal最小生成树,从小到大每次加完一种权值的边后处理这里的询问。这样的复杂度为 \(O(n\log n\alpha(n))\) 。

然而现在是强制在线,引入一种新的数据结构——kruskal重构树(我也不知道为什么是这个名字)。它的思想非常简单。kruskal的时候我们是从小到大添加边,那么我们还是用类似并查集的方法,每次合并的时候新建一个点,把需要合并的两个点都挂在这个新点的下面,新点的权值为这条边的边权。这样我们会得到一棵树,它满足新点的结构组成一个大根堆。那么如果我们要查询从一个点出发,边权小于等于 \(d\) 的边能到达的点,其实就是从这个点不断往上跳,保证经过的每个新点 \(p\) 的权值都小于等于 \(d\) 。最后走到的这个新点的整个子树就是我们能够走到的点。

由于这棵树的高度没有保证,我们把这棵树建出来,做一个可持久化线段树的合并,倍增跳到合法的 \(p\) ,查询子树即可。复杂度为 \(O(n\log n)\),非常优秀。

其实这题我还有一个想法。我们每次不新建点,而是按秩合并并查集,在每个点上开一个map,记录一下每个权值对应的线段树根,每次暴力往上跳到合法位置,在map上lower_bound以下,直接查询即可。复杂度也为 \(O(n\log n)\) 。

代码

这不是kruskal重构树,是我那个方法。

#include<bits/stdc++.h>
using namespace std;
inline char nchar() {
static const int bufl=1<<20;
static char buf[bufl],*a,*b;
return a==b && (b=(a=buf)+fread(buf,1,bufl,stdin),a==b)?EOF:*a++;
}
inline int read() {
int x=0,f=1;
char c=nchar();
for (;!isdigit(c);c=nchar()) if (c=='-') f=-1;
for (;isdigit(c);c=nchar()) x=x*10+c-'0';
return x*f;
}
inline void write(int x) {
if (x==0) puts("0"); else if (x==-1) puts("-1"); else {
static char s[20];
int tot=0;
for (;x;x/=10) s[++tot]=x%10+'0';
while (tot) putchar(s[tot--]);
puts("");
}
}
const int maxn=1e5+1;
const int maxm=5e5+1;
const int inf=1e9+7;
int h[maxn],c[maxn],ls=0,n,m,q;
pair<int,pair<int,int> > e[maxm];
namespace sgt {
const int maxp=5e6+1;
struct node {
int l,r,size;
} t[maxp];
int tot=0;
void inc(int &x,int l,int r,int p) {
if (!x) x=++tot;
++t[x].size;
if (l==r) return;
int mid=(l+r)>>1;
p<=mid?inc(t[x].l,l,mid,p):inc(t[x].r,mid+1,r,p);
}
inline void inc(int &x,int p) {inc(x,1,ls,p);}
int merge(int x,int y,int l,int r) {
if (!x) return y;
if (!y) return x;
int p=++tot,mid=(l+r)>>1;
if (l==r) {
t[p].size=t[x].size+t[y].size;
return p;
}
t[p].l=merge(t[x].l,t[y].l,l,mid);
t[p].r=merge(t[x].r,t[y].r,mid+1,r);
t[p].size=t[t[p].l].size+t[t[p].r].size;
return p;
}
int query(int x,int l,int r,int k) {
if (l==r) return l;
int mid=(l+r)>>1;
return k<=t[t[x].l].size?query(t[x].l,l,mid,k):query(t[x].r,mid+1,r,k-t[t[x].l].size);
}
inline int query(int x,int k) {
if (k==0 || k>t[x].size) return -1;
return query(x,1,ls,t[x].size-k+1);
}
}
namespace st {
int rk[maxn],f[maxn],w[maxn],last[maxn];
map<int,int> mp[maxn];
inline void init(int n) {for (int i=1;i<=n;++i) rk[i]=0,w[i]=-1,f[i]=i;}
int find(int x,int d=inf) {
for (;f[x]!=x && w[x]<=d;x=f[x]);
return x;
}
int uni(int x,int y,int v) {
if (rk[x]>rk[y]) swap(x,y);
f[x]=y,w[x]=v,rk[y]+=(rk[x]==rk[y]);
int tmp=mp[y].rbegin()->second;
int &rt=mp[y][v];
rt=sgt::merge(rt?rt:tmp,mp[x].rbegin()->second,1,ls);
}
int query(int x,int val,int k) {
int fx=find(x,val);
map<int,int>::iterator it=mp[fx].upper_bound(val);
int rt=(--it)->second;
return sgt::query(rt,k);
}
}
inline bool cmp(int x,int y) {return st::rk[x]<st::rk[y];}
int main() {
#ifndef ONLINE_JUDGE
freopen("test.in","r",stdin);
#endif
n=read(),m=read(),q=read();
st::init(n);
for (int i=1;i<=n;++i) c[++ls]=h[i]=read();
sort(c+1,c+ls+1),ls=unique(c+1,c+ls+1)-c-1;
for (int i=1;i<=n;++i) h[i]=lower_bound(c+1,c+ls+1,h[i])-c;
for (int i=1;i<=m;++i) {
int x=read(),y=read(),w=read();
e[i]=make_pair(w,make_pair(x,y));
}
sort(e+1,e+m+1);
for (int i=1;i<=n;++i) {
sgt::inc(st::mp[i][-1],h[i]);
}
for (int i=1;i<=m;++i) {
int &w=e[i].first,&x=e[i].second.first,&y=e[i].second.second;
int fx=st::find(x),fy=st::find(y);
if (fx!=fy) st::uni(fx,fy,w);
}
int ans=0;
while (q--) {
int x=read()^ans,val=read()^ans,k=read()^ans;
ans=st::query(x,val,k);
if (ans!=-1) ans=c[ans];
write(ans);
if (ans==-1) ans=0;
}
return 0;
}

bzoj3545-bzoj3551-Peaks的更多相关文章

  1. bzoj3545/bzoj3551 [ONTAK2010]Peaks/Peaks加强版

    bzoj3545/bzoj3551 [ONTAK2010]Peaks/Peaks加强版 传送门:bzoj  bzoj wdnmd为什么加强版不是权限题原题却是啊 3545: [ONTAK2010]Pe ...

  2. [您有新的未分配科技点][BZOJ3545&BZOJ3551]克鲁斯卡尔重构树

    这次我们来搞一个很新奇的知识点:克鲁斯卡尔重构树.它也是一种图,是克鲁斯卡尔算法求最小生成树的升级版首先看下面一个问题:BZOJ3545 Peaks. 在Bytemountains有N座山峰,每座山峰 ...

  3. bzoj3545: [ONTAK2010]Peaks 重构树 主席树

    题目链接 bzoj3545: [ONTAK2010]Peaks 题解 套路重构树上主席树 代码 #include<cstdio> #include<algorithm> #de ...

  4. [BZOJ3551]Peaks

    [BZOJ3551]Peaks BZOJ luogu 建Kruskal重构树,点权为边权 按dfn序建出主席树 倍增找到能跳到的最浅的祖先 主席树查询一下 #include<bits/stdc+ ...

  5. 【bzoj3545/bzoj3551】[ONTAK2010]Peaks/加强版 Kruskal+树上倍增+Dfs序+主席树

    bzoj3545 题目描述 在Bytemountains有N座山峰,每座山峰有他的高度h_i.有些山峰之间有双向道路相连,共M条路径,每条路径有一个困难值,这个值越大表示越难走,现在有Q组询问,每组询 ...

  6. bzoj3545 [ONTAK2010]Peaks、bzoj3551 [ONTAK2010]Peaks加强版

    题目描述: bzoj3545,luogu bzoj3551 题解: 重构树+线段树合并. 可以算是板子了吧. 代码(非强制在线): #include<cstdio> #include< ...

  7. 【BZOJ3545&BZOJ3551】Peaks(kruskal重构树,主席树,dfs序)

    题意:在Bytemountains有N座山峰,每座山峰有他的高度h_i. 有些山峰之间有双向道路相连,共M条路径,每条路径有一个困难值,这个值越大表示越难走, 现在有Q组询问,每组询问询问从点v开始只 ...

  8. BZOJ3545 [ONTAK2010]Peaks kruskal 并查集 主席树 dfs序

    欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ3545 题意概括 Description 在Bytemountains有N座山峰,每座山峰有他的高度 ...

  9. bzoj3551 Peaks加强版

    这个题--感觉离线和在线的代码难度差不多(pb_ds不要说话). 离线的话,就是把所有询问按照w排个序,然后一边Kruskal+平衡树启发式合并一边回答询问就好了. 在线也不难写.首先Kruskal重 ...

  10. bzoj3545: [ONTAK2010]Peaks

    Description 在Bytemountains有N座山峰,每座山峰有他的高度h_i.有些山峰之间有双向道路相连,共M条路径,每条路径有一个困难值,这个值越大表示越难走,现在有Q组询问,每组询问询 ...

随机推荐

  1. 20155338 2016-2017-2 《JAVA程序设计》课程总结

    ---恢复内容开始--- 20155338 2016-2017-2 <JAVA程序设计>课程总结 每周作业汇总: 预备作业一 预备作业二 预备作业三 <JAVA程序设计>第一周 ...

  2. Azkaban系统的安装和分析。

    Azkaban系统是一个数据处理的很好用的工具,可以用来运行hadoop任务,管理hdfs,可以进行schedule任务调度,总体来说功能还是很强大的. 研究了一下azkaban,做了以下总结性的东西 ...

  3. django项目的配置文件settings.py详解

    我们创建好了一个Python项目(mysite/)之后,需要在项目中添加模块应用(polls/),在模块应用中添加处理功能逻辑,如添加模块中的视图处理函数(polls.views.index()),这 ...

  4. python全栈开发-前方高能-函数

    python_day_9 一.今日主要内容 函数: 函数:对功能的封装 语法: def 函数名(形参): 函数体 函数名(实参) 函数名:命名规则和变量一样 函数的返回值: return, 函数执行完 ...

  5. 第五章Web应用与应用层协议

    Web应用与应用层协议 本篇博文中的主要参考文献是<计算机网络高级教程>,分别是吴功宜老先生和吴英教授合著.这部教程是我研究生老师所推荐的网络必读科目,由于该教程讲解的基础知识详细,但内容 ...

  6. 天马行空DevOps-Dev平台建设概述

    概述 DevOps(Development和Operations的组合词)是一组过程.方法与系统的统称,用于促进开发(应用程序/软件工程).技术运营和质量保障(QA)部门之间的沟通.协作与整合.它是一 ...

  7. 从零开始的Python学习Episode 11——装饰器

    装饰器 装饰器是用来处理其他函数的函数,主要作用是在不修改原有函数的情况下添加新的功能,装饰器的返回值也是一个函数对象. 简单的装饰器 import time def show_time(f): de ...

  8. gets函数的完美替代

    众所周知 在C语言中scanf用来读取一行字符串时遇到空格或回车会停止 而若要读入一行带空格的字符串时 有些人会用gets来代替 然而,gets的最大问题在于:会读取超过数组长度上限个字符,而超出长度 ...

  9. loadrunner socket协议问题归纳(3)

    摘要:通过实例讲解loadrunner中的socket协议性能测试的一种测试方法,如何不依赖loadrunner既定规则,自行控制收发数据包 关键词:Loadrunner,socket,自行控制,收发 ...

  10. Echarts数据可视化全解

    点击进入 Echarts数据可视化全解