bzoj3083 遥远的国度 树链剖分+树上lca

由于懒得重打一遍直接贴上一题的剖分模板，结果因为上一题是单点查询所以一直T，T了5发才发现我是有多sb。。

明显不能直接换根，不然我全部操作都是换根你不就炸了。
假设开始的时候root，现在的根是newroot，查询的是x,那么明显有:
1.如果x=newroot，查询全部
2.如果newroot在x外，不受影响，直接查。
3.newroot在x内，查询x子树以外的，即1–pos[x]-1和mxpos[x]+1–n

#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define inf 0x7fffffff  
using namespace std;
const int N=1e5+5;
int n,m;
int head[N<<1],next[N<<1],go[N<<1];
int siz[N],dep[N],fa[N][20];
int pos[N],top[N],a[N],mxpos[N];
int tot=0,sz;
struct node
{
    int l,r,mn,lazy;
}t[1000000];
inline void add(int x,int y)
{
    go[++tot]=y;
    next[tot]=head[x];
    head[x]=tot;
}

inline void build(int x,int l,int r)
{
    t[x].l=l,t[x].r=r;
    t[x].mn=inf;
    if (l==r)return;
    int mid=(l+r)>>1;
    build(x<<1,l,mid);
    build(x<<1|1,mid+1,r);
}
inline void pushdown(int x)
{
    if (t[x].l==t[x].r)return;
    if (t[x].lazy)
    {
        t[x<<1].lazy=t[x<<1|1].lazy=t[x<<1].mn=t[x<<1|1].mn=t[x].lazy;
        t[x].lazy=0;
    }
}
inline void pushup(int x)
{
    t[x].mn=min(t[x<<1].mn,t[x<<1|1].mn);
}
inline void change(int x,int l,int r,int z)
{
    pushdown(x);
    if (t[x].l==l&&t[x].r==r)
    {
        t[x].lazy=t[x].mn=z;
        return;
    }
    int mid=(t[x].l+t[x].r)>>1;
    if (r<=mid)change(x<<1,l,r,z);
    else if (l>mid)change(x<<1|1,l,r,z);
    else
    {
        change(x<<1,l,mid,z);
        change(x<<1|1,mid+1,r,z);
    }
    pushup(x);
}
inline void dfs1(int x)
{
    fo(i,1,16)
    if (fa[x][i-1])fa[x][i]=fa[fa[x][i-1]][i-1];
    siz[x]=1;
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        if (v!=fa[x][0])
        {
            dep[v]=dep[x]+1;
            fa[v][0]=x;
            dfs1(v);
            siz[x]+=siz[v];
        }
    }
}
inline void dfs2(int x,int chain)
{
    int k=0;
    sz++;
    pos[x]=sz;
    top[x]=chain;
    mxpos[x]=sz;
    change(1,sz,sz,a[x]);
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        if (v!=fa[x][0]&&siz[v]>siz[k])
            k=v;
    }
    if (!k)return;
    dfs2(k,chain);
    mxpos[x]=max(mxpos[x],mxpos[k]);
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        if (v!=fa[x][0]&&k!=v)
        {
            dfs2(v,v);
            mxpos[x]=max(mxpos[x],mxpos[v]);
        }
    }
}
inline int query(int x,int l,int r)
{
    pushdown(x);
    if (t[x].l==l&&t[x].r==r)return t[x].mn;
    int mid=(t[x].l+t[x].r)>>1;
    if (r<=mid)return query(x<<1,l,r);
    else if (l>mid)return query(x<<1|1,l,r);
    else return min(query(x<<1,l,mid),query(x<<1|1,mid+1,r));
}
inline void solve(int x,int y,int z)
{
    while(top[x]!=top[y])
    {
        if (dep[top[x]]<dep[top[y]])swap(x,y);
        change(1,pos[top[x]],pos[x],z);
        x=fa[top[x]][0];
    }
    if (dep[x]<dep[y])swap(x,y);
    change(1,pos[y],pos[x],z);
}
inline int getlca(int x,int y)
{
    if (dep[x]<dep[y])swap(x,y);
    fd(i,16,0)
        if (dep[fa[x][i]]>=dep[y])x=fa[x][i];
    if (x==y)return x;
    fd(i,16,0)
    if (fa[x][i]!=fa[y][i])
    {
        x=fa[x][i],y=fa[y][i];
    }
    return fa[x][0];
}
int getv(int x,int y)
{
    fd(i,16,0)
    if (dep[fa[x][i]]>dep[y])x=fa[x][i];
    return x;
}
inline int solveq(int x,int root)
{
    if (x==root)return query(1,1,n);
    if (getlca(x,root)!=x)return query(1,pos[x],mxpos[x]);
    int v=getv(root,x),ans=inf;
    if (pos[v]>1)ans=min(ans,query(1,1,pos[v]-1));
    if (mxpos[v]<n)ans=min(ans,query(1,mxpos[v]+1,n));
    return ans;
}
int main()
{
    scanf("%d%d",&n,&m);
    fo(i,1,n-1)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        add(x,y);
        add(y,x);
    }
    build(1,1,n);
    fo(i,1,n)scanf("%d",&a[i]);
    int root;
    scanf("%d",&root);
    dep[root]=1;
    dfs1(root);
    dfs2(root,root);
    int nowroot=root;
    fo(i,1,m)
    {
        int id;
        scanf("%d",&id);
        if (id==1)
        {
            int x;
            scanf("%d",&x);
            nowroot=x;
        }
        else if (id==2)
        {
            int x1,x2,v;
            scanf("%d%d%d",&x1,&x2,&v);
            solve(x1,x2,v);
        }
        else {
            int x;
            scanf("%d",&x);
            printf("%d\n",solveq(x,nowroot));
        }
    }
}