BZOJ

思路

我觉得这题可持久化线段树合并也可以做

我觉得这题建出最小生成树之后动态点分治+线段树也可以做

还是学习一下Kruskal重构树吧……

Kruskal重构树,就是在做最小生成树的时候,如果一条边\(e\)被选中了,就让那两个连通块的根都连向它,变成新的根。显然,最后会做出一个二叉树,其中叶子是点、非叶子节点是边,且边权形成了一个堆的形式。

分析一下性质,发现“用不超过某个权值的边构成的连通块”其实就是从一个点往上跳,最后的某一棵子树。

于是倍增+主席树区间k大,这题就做完了。

Kruskal重构树的例题还有[NOI2018]归程。只要知道这个东西那题似乎就一sb题……

代码

#include<bits/stdc++.h>
//clock_t t=clock();
namespace my_std{
using namespace std;
// #define pii pair<int,int>
// #define fir first
// #define sec second
// #define MP make_pair
#define rep(i,x,y) for (int i=(x);i<=(y);i++)
#define drep(i,x,y) for (int i=(x);i>=(y);i--)
#define go(x) for (int i=head[x];i;i=edge[i].nxt)
// #define templ template<typename T>
#define sz 606060
typedef long long ll;
typedef double db;
// mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
// templ inline T rnd(T l,T r) {return uniform_int_distribution<T>(l,r)(rng);}
// templ inline bool chkmax(T &x,T y){return x<y?x=y,1:0;}
// templ inline bool chkmin(T &x,T y){return x>y?x=y,1:0;}
inline void read(int& t)
{
t=0;char f=0,ch=getchar();double d=0.1;
while(ch>'9'||ch<'0') f|=(ch=='-'),ch=getchar();
while(ch<='9'&&ch>='0') t=t*10+ch-48,ch=getchar();
t=(f?-t:t);
}
// template<typename T,typename... Args>inline void read(T& t,Args&... args){read(t); read(args...);}
// char __sr[1<<21],__z[20];int __C=-1,__zz=0;
// inline void Ot(){fwrite(__sr,1,__C+1,stdout),__C=-1;}
// inline void print(register int x)
// {
// if(__C>1<<20)Ot();if(x<0)__sr[++__C]='-',x=-x;
// while(__z[++__zz]=x%10+48,x/=10);
// while(__sr[++__C]=__z[__zz],--__zz);__sr[++__C]='\n';
// }
// void file()
// {
// #ifdef NTFOrz
// freopen("a.in","r",stdin);
// #endif
// }
// inline void chktime()
// {
// #ifndef ONLINE_JUDGE
// cout<<(clock()-t)/1000.0<<'\n';
// #endif
// }
// #ifdef mod
// ll ksm(ll x,int y){ll ret=1;for (;y;y>>=1,x=x*x%mod) if (y&1) ret=ret*x%mod;return ret;}
// ll inv(ll x){return ksm(x,mod-2);}
// #else
// ll ksm(ll x,int y){ll ret=1;for (;y;y>>=1,x=x*x) if (y&1) ret=ret*x;return ret;}
// #endif
// inline ll mul(ll a,ll b){ll d=(ll)(a*(double)b/mod+0.5);ll ret=a*b-d*mod;if (ret<0) ret+=mod;return ret;}
}
using namespace my_std; int n,m,Q;
int val[sz],V[sz];
struct hhh{int u,v,w;const bool operator < (const hhh &x) const {return w<x.w;}}e[sz];
struct hh{int t,nxt;}edge[sz<<1];
int head[sz],ecnt;
void make_edge(int f,int t)
{
edge[++ecnt]=(hh){t,head[f]};
head[f]=ecnt;
edge[++ecnt]=(hh){f,head[t]};
head[t]=ecnt;
} int rt;
int fa[sz][25],dfn[sz],low[sz],id[sz],T;
#define v edge[i].t
void dfs(int x,int f)
{
id[dfn[x]=++T]=x;
fa[x][0]=f;
rep(i,1,20) fa[x][i]=fa[fa[x][i-1]][i-1];
go(x) if (v!=f) dfs(v,x);
low[x]=T;
}
#undef v int f[sz];
int getfa(int x){return x==f[x]?x:f[x]=getfa(f[x]);}
void merge(int x,int y){f[x]=y;make_edge(x,y);} int root[sz];
#define Tree sz*4
int size[Tree],ls[Tree],rs[Tree],cnt;
#undef Tree
#define lson ls[k],l,mid
#define rson rs[k],mid+1,r
void insert(int &k,int l,int r,int x,int pre)
{
k=++cnt;size[k]=size[pre]+1;ls[k]=ls[pre],rs[k]=rs[pre];
if (l==r) return;
int mid=(l+r)>>1;
if (x<=mid) insert(lson,x,ls[pre]);
else insert(rson,x,rs[pre]);
}
int query(int k1,int k2,int l,int r,int K)
{
if (l==r) return l;
int mid=(l+r)>>1,S=size[rs[k2]]-size[rs[k1]];
if (S>=K) return query(rs[k1],rs[k2],mid+1,r,K);
return query(ls[k1],ls[k2],l,mid,K-S);
} int main()
{
// file();
read(n),read(m),read(Q);
rep(i,1,n) read(val[i]),V[i]=val[i];
sort(V+1,V+n+1);int _n=unique(V+1,V+n+1)-V-1;
rep(i,1,n) val[i]=lower_bound(V+1,V+_n+1,val[i])-V;
int x,y,z;
rep(i,1,m) read(x),read(y),read(z),e[i]=(hhh){x,y,z};
sort(e+1,e+m+1);
rep(i,1,n+m) f[i]=i;
rep(i,1,m)
{
x=e[i].u,y=e[i].v;
x=getfa(x),y=getfa(y);
if (x==y) continue;
merge(x,i+n),merge(y,i+n);
}
rt=getfa(1);
dfs(rt,0);
rep(i,1,n+n-1) if (id[i]<=n) insert(root[i],1,n,val[id[i]],root[i-1]); else root[i]=root[i-1];
int ans=-1;
while (Q--)
{
read(x),read(y),read(z);
if (ans!=-1) x^=ans,y^=ans,z^=ans;
drep(i,20,0) if (fa[x][i]&&e[fa[x][i]-n].w<=y) x=fa[x][i];
if (size[root[low[x]]]-size[root[dfn[x]-1]]<z) ans=-1;
else ans=V[query(root[dfn[x]-1],root[low[x]],1,n,z)];
printf("%d\n",ans);
}
return 0;
}

BZOJ3551 Peaks加强版 [Kruskal重构树,主席树]的更多相关文章

  1. [BZOJ3551][ONTAK2010]Peaks(加强版)(Kruskal重构树,主席树)

    3551: [ONTAK2010]Peaks加强版 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 2438  Solved: 763[Submit][ ...

  2. 【BZOJ-3545&3551】Peaks&加强版 Kruskal重构树 + 主席树 + DFS序 + 倍增

    3545: [ONTAK2010]Peaks Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1202  Solved: 321[Submit][Sta ...

  3. BZOJ.3551.[ONTAK2010]Peaks加强版(Kruskal重构树 主席树)

    题目链接 \(Description\) 有n个座山,其高度为hi.有m条带权双向边连接某些山.多次询问,每次询问从v出发 只经过边权<=x的边 所能到达的山中,第K高的是多少. 强制在线. \ ...

  4. luoguP4197:Peaks(Kruskal重构树+主席树)或者(点分树+离线)

    题意:有N座山,M条道路.山有山高,路有困难值(即点权和边权).现在Q次询问,每次给出(v,p),让求从v出发,只能结果边权<=p的边,问能够到达的山中,第K高的高度(从大到小排序). 思路:显 ...

  5. BZOJ 3551: [ONTAK2010]Peaks加强版 [Kruskal重构树 dfs序 主席树]

    3551: [ONTAK2010]Peaks加强版 题意:带权图,多组询问与一个点通过边权\(\le lim\)的边连通的点中点权k大值,强制在线 PoPoQQQ大爷题解传送门 说一下感受: 容易发现 ...

  6. 【BZOJ 3551】[ONTAK2010] Peaks加强版 Kruskal重构树+树上倍增+主席树

    这题真刺激...... I.关于Kruskal重构树,我只能开门了,不过补充一下那玩意还是一棵满二叉树.(看一下内容之前请先进门坐一坐) II.原来只是用树上倍增求Lca,但其实树上倍增是一种方法,L ...

  7. BZOJ 3551: [ONTAK2010]Peaks加强版 Kruskal重构树+dfs序+主席树+倍增

    建出来 $Kruskal$ 重构树. 将询问点向上跳到深度最小,且合法的节点上. 那么,得益于重构树优美的性质,这个最终跳到的点为根的所有子节点都可以与询问点互达. 对于子树中求点权第 $k$ 大的问 ...

  8. BZOJ3545&3551[ONTAK2010]Peaks——kruskal重构树+主席树+dfs序+树上倍增

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

  9. luogu4197 Peaks (kruskal重构树+主席树)

    按照边权排序建出kruskal重构树,每次就变成了先找一个权值<=x的最远的祖先,然后看这个子树的第k小.离散化一下,在dfs序上做主席树即可 而且只需要建叶节点的主席树 注意输出的是第k小点的 ...

随机推荐

  1. idea 中激活 JRebel

    JRebel介绍: JRebel是一款JVM插件,它使得Java代码修改后不用重启系统,立即生效.IDEA上原生是不支持热部署的,一般更新了 Java 文件后要手动重启 Tomcat 服务器,修改才能 ...

  2. java程序员必须熟悉的一些操作

    1.mysql数据库服务启动命令 /etc/init.d/mysqld start     --启动命令      mysql数据库安装方法参考           http://www.blogja ...

  3. android使用http3

    http3的github地址: https://github.com/cloudflare/quiche

  4. python pip 切换阿里云镜像

    示例: pip install xxx -i http://mirrors.aliyun.com/pypi/simple/ --trusted-host mirrors.aliyun.com

  5. 最简单 无返回值 无参数 sql server存储过程

    CREATE proc Upadte_stateas update Table_1 set [state]=2 where id in (select id from Table_1 where st ...

  6. Matlab图像处理基础知识

    Matlab图像处理基础知识 Matlab的图片以矩阵的形式存储,矩阵的行列值为图片的行列的色彩值. 1图像表达方式: 像素索引 图像被视为离散单元.如使用I(2,2)可以获取第二行第二列的像素值 空 ...

  7. 父元素设置min-height子元素设置100%问题

    问题:父元素设置min-height子元素高度设置100%取不到值,这是因为子元素 div设置 height:100%: 只有当父级元素满足min-height:1000px;设置的条件才触发: 浏览 ...

  8. oracle 查询 10题

    说明:表数据来自oracle 初始用户之一scott里面的三个初始表:emp,dept,salgrade --1.查询员工表中工资最高的雇员的员工号.员工姓名.工资和部门号. select empno ...

  9. AutoMapper 初次使用心得

    本例以asp.net webform为例: 结构: 主要代码:AutoMapperConfig 类 public class AutoMapperConfig { public static void ...

  10. DELL R730 做raid10

    1.服务器开机,在出现下图提示时,同时按着<ctrl >+ < R >键,即可进入配置界面 2.会进入下图 3.按上下键到第一项PERC H730P MINI ,按F2,选择c ...