bzoj 3551: [ONTAK2010]Peaks加强版

题意：

同3545，但强制在线。

题解：

这里要用到一个有趣的东西：kruskal重构树。
就是在做kruskal的时候，将两个联通块合并时，增加一个新的点，且将两个联通块接向它，点权为当前边的边权。
下转popoqqq：
这样生成的树有一些十分优美的性质：
1.二叉树(好吧这题意义不大)
2.原树与新树两点间路径上边权(点权)的最大值相等
3.子节点的边权小于等于父亲节点(大根堆)
4.原树中两点之间路径上边权的最大值等于新树上两点的LCA的点权
于是对于每个询问 我们从这个询问向上倍增寻找深度最小的点权小于等于x的点 易证这个节点的子树就是v所能到达的所有点。

然后dfs序+主席树。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int inf=1e9;
int n,m,Q,h[1000010];
struct node{
    int x,y,c,next;
}e[500010],a[200010];int len=0,last[200010];
struct trnode{
    int dep,fa[17],C[17];
}tr[200010];int tot=0,root[200010];
struct TRnode{
    int lc,rc,c;
}TR[6100010];
int fa[200010],num,w[200010],ys[200010],z=0,la[200010],A[200010];
int findfa(int x) {return fa[x]==x?x:fa[x]=findfa(fa[x]);}
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
bool cmp(node a,node b) {return a.c<b.c;}
void update(int &x,int fr,int l,int r,int k)
{
    x=++tot;TR[x]=TR[fr];
    if(l==r) {TR[x].c++;return;}
    int mid=(l+r)/2;
    if(k<=mid) update(TR[x].lc,TR[fr].lc,l,mid,k);
    else update(TR[x].rc,TR[fr].rc,mid+1,r,k);
    TR[x].c=TR[TR[x].lc].c+TR[TR[x].rc].c;
}
int findans(int x1,int x2,int l,int r,int k)
{
    if(TR[x2].c-TR[x1].c<k) return -1;
    if(l==r) return l;
    int mid=(l+r)/2,rcc=TR[TR[x2].rc].c-TR[TR[x1].rc].c;
    if(rcc>=k) return findans(TR[x1].rc,TR[x2].rc,mid+1,r,k);
    else return findans(TR[x1].lc,TR[x2].lc,l,mid,k-rcc);
}
void ins(int x,int y)
{
    a[++len].y=y;
    a[len].next=last[x];last[x]=len;
}
void dfs(int x,int fa)
{
    tr[x].dep=tr[fa].dep+1;tr[x].fa[0]=fa;tr[x].C[0]=w[fa];
    for(int i=1;(1<<i)<=tr[x].dep;i++)
    {
        tr[x].fa[i]=tr[tr[x].fa[i-1]].fa[i-1];
        tr[x].C[i]=max(tr[x].C[i-1],tr[tr[x].fa[i-1]].C[i-1]);
    }
    ys[x]=la[x]=++z;
    for(int i=last[x];i;i=a[i].next)
    {
        int y=a[i].y;
        dfs(y,x);la[x]=la[y];
    }
}
int find(int x,int v)
{
    for(int i=16;i>=0;i--)
        if((1<<i)<=tr[x].dep&&tr[x].C[i]<=v) x=tr[x].fa[i];
    return x;
}
int main()
{
    n=read();m=read();Q=read();num=n;
    for(int i=1;i<=n;i++) h[i]=read(),fa[i]=i;
    for(int i=1;i<=m;i++) e[i].x=read(),e[i].y=read(),e[i].c=read();
    sort(e+1,e+m+1,cmp);
    for(int i=1;i<=m;i++)
    {
        int x=findfa(e[i].x),y=findfa(e[i].y);
        if(x!=y)
        {
            w[++num]=e[i].c;ins(num,x);ins(num,y);
            fa[num]=fa[x]=fa[y]=num;
        }
    }
    tr[0].dep=-1;dfs(num,0);
    for(int i=1;i<=n;i++) A[ys[i]]=h[i];
    for(int i=1;i<=num;i++)
    {
        if(A[i]==0) root[i]=root[i-1];
        else update(root[i],root[i-1],1,inf,A[i]);
    }
    int ans=0;
    while (Q--)
    {
        int v=read(),x=read(),k=read();
        if (ans>-1) v^=ans,x^=ans,k^=ans;
        int p=find(v,x);ans=findans(root[ys[p]-1],root[la[p]],1,inf,k);
        printf("%d\n",ans);
    }
}