HDU 3078 LCA转RMQ

点击打开链接

题意：给一个树和每个点的权值，然后m询问a b c，若a为0，则将b点的权值改为c，否则问你的是从b走到c的所有权值中的第a大的数，若不足a个输出那个

思路：就是用LCA转RMQ就可以了，每次的询问就是直接找到它们的最近公共祖先，然后两个点都向上找到最近公共祖先，沿途的所有的点全部记录下来，最后排个序输出第K大就行了

#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int maxn=80010;
vector<int>G[maxn];
bool vis[maxn];
int L[maxn*2],E[maxn*2],H[maxn],dis[maxn],dp[2*maxn][20],num[maxn],pre[maxn];
//H为每个元素第一次出现的位置，E为dfs遍历的序号，L为深度，dis为每个点到根的距离
//dp[i][j]代表的是i到j深度最小的位置
int k,n;
void dfs(int t,int deep){
k++;E[k]=t;L[k]=deep;H[t]=k;//遍历时叶子在E中只会出现一次，而非叶子节点会出现多次
for(unsigned int i=0;i<G[t].size();i++){
int tt=G[t][i];
if(!vis[tt]){//出现过就不再继续，保证dis的大小和H第一次出现的位置
pre[tt]=t;
vis[tt]=1;
dfs(tt,deep+1);
k++;E[k]=t;L[k]=deep;
}
}
}
void RMQ_init(){//RMQ的ST算法
for(int i=1;i<=2*n-1;i++) dp[i][0]=i;
for(int i=1;(1<<i)<=2*n-1;i++){
for(int j=1;j+(1<<i)-1<=2*n-1;j++){
if(L[dp[j][i-1]]<L[dp[j+(1<<(i-1))][i-1]]) dp[j][i]=dp[j][i-1];
else dp[j][i]=dp[j+(1<<(i-1))][i-1];
}
}
}
int RMQ(int le,int ri){
le=H[le];ri=H[ri];//找到第一次出现的位置
if(le>ri) swap(le,ri);//le大于ri说明ri在le之前出现，交换
int kk=0;
while((1<<(kk+1))<=ri-le+1) kk++;
if(L[dp[le][kk]]<L[dp[ri-(1<<kk)+1][kk]]) return E[dp[le][kk]];
else return E[dp[ri-(1<<kk)+1][kk]];//返回的是dp位置的值
}
int tmp[maxn];
int main(){
int q,a,b,c;
while(scanf("%d%d",&n,&q)!=-1){
for(int i=1;i<=n;i++) scanf("%d",&num[i]);
for(int i=0;i<maxn;i++) G[i].clear();
memset(vis,0,sizeof(vis));
for(int i=0;i<n-1;i++){
scanf("%d%d",&a,&b);
G[a].push_back(b);
G[b].push_back(a);
}
k=0;vis[1]=1;
dfs(1,1);RMQ_init();
while(q--){
scanf("%d%d%d",&a,&b,&c);
if(a==0){
num[b]=c;continue;
}
int en=RMQ(b,c);
int sum=0;
while(b!=en){
tmp[sum++]=num[b];
b=pre[b];
}
while(c!=en){
tmp[sum++]=num[c];
c=pre[c];
}
tmp[sum++]=num[en];
if(sum<a) printf("invalid request!\n");
else{
sort(tmp,tmp+sum);
printf("%d\n",tmp[sum-a]);
}
}
}
return 0;
}