【bzoj3083】 遥远的国度

比较难想的一道链剖（其实就是不那么裸的一道链剖），看来我这个zz只会刷裸题
op2就是树链剖分的基本应用，关键是op3操作
我们可以对op3进行讨论
我们可以将结果分为3种情况
1.点x就是root，这种情况需要查询整个树
2.点x是在root的子树上（即点x与root的LCA是root），这样直接查询x子树的就行
3.root在点x的子树上（即点x与root的LCA是x），此时我们需要求出处于x与root这条链上的x的下面的一个点，需要查询在树中root到这个点的补集，然后这道题就可以直接树剖做了.
平时多思考一下，不要光想着看题解
注意对于dfs序建的线段树顺序是点的dfs序，查询的时候一点要写num[x]，因为这个错误查了好久.

#include<iostream>#include<cstdio>#include<cstring>#include<string>#include<algorithm>//树剖的时候没写num using namespace std;const int N=100010,inf=0x7fffffff;int n,m,sz,root,te,qx;int bin[20],fa[N][20],h[N],tp[N],val[N];int son[N],num[N],tree[N],size[N],head[N];struct edge{    int v,next;}e[200010];struct seg{    int l,r,tag,mn;}tr[800010];inline int F(){    register int aa,bb;register char ch;    while(ch=getchar(),(ch<'0'||ch>'9')&&ch!='-');ch=='-'?aa=bb=0:(aa=ch-'0',bb=1);    while(ch=getchar(),ch<='9'&&ch>='0')aa=(aa<<3)+(aa<<1)+ch-'0';return bb?aa:-aa;}void clear(){    memset(head,0,sizeof(head));    memset(son,0,sizeof(son));    memset(fa,0,sizeof(fa));    size[0]=0,fa[1][0]=h[1]=bin[0]=1;    for (int i=1;i<=18;++i)    bin[i]=bin[i-1]<<1;}void add(int u,int v){    e[++te].v=v;    e[te].next=head[u];    head[u]=te;}void dfs1(int x){    size[x]=1;    for (int i=1;i<=16;++i)    if (h[x]>bin[i])fa[x][i]=fa[fa[x][i-1]][i-1];    else break;    for (int i=head[x];i;i=e[i].next)    {        int v=e[i].v;        if (v==fa[x][0])continue;        h[v]=h[x]+1;        fa[v][0]=x;        dfs1(v);        size[x]+=size[v];        if (size[v]>size[son[x]])son[x]=v;    }} void dfs2(int x,int chain){    tp[x]=chain;    num[x]=++sz;    if (!son[x])return ;    dfs2(son[x],chain);    for (int i=head[x];i;i=e[i].next)    {        int v=e[i].v;        if (v==fa[x][0]||v==son[x])continue;        dfs2(v,v);      } }int lca(int x,int y){    if (h[x]<h[y])swap(x,y);    int t=h[x]-h[y];    for (int i=0;i<=16;i++)    if(t&bin[i])x=fa[x][i];    for (int i=16;i>=0;--i)    if (fa[x][i]!=fa[y][i])    x=fa[x][i],y=fa[y][i];    if (x==y)return x;    return fa[x][0];}int lca1(int x,int y){    if (h[x]<h[y])swap(x,y);    int t=h[x]-h[y]-1;    for (int i=0;i<=16;++i)    if (bin[i]&t)x=fa[x][i];    return x;}void updata(int k){tr[k].mn=min(tr[k<<1].mn,tr[k<<1|1].mn);}void pushdown(int k){    if (!tr[k].tag||tr[k].l==tr[k].r)return;    int t=tr[k].tag;    tr[k<<1].mn=tr[k<<1|1].mn=t;    tr[k<<1].tag=tr[k<<1|1].tag=t;    tr[k].tag=0;}void build(int k,int l,int r){    tr[k].l=l;tr[k].r=r;tr[k].tag=0;    if (l==r)    {        tr[k].mn=val[tree[l]];        return ;    }    int mid=(l+r)>>1;    build(k<<1,l,mid);    build(k<<1|1,mid+1,r);    updata(k);}void query(int k,int x,int y){    pushdown(k);    int l=tr[k].l,r=tr[k].r;    if(x<=l&&r<=y)    {        qx=min(qx,tr[k].mn);        return ;    }    int mid=(l+r)>>1;    if (x<=mid)query(k<<1,x,y);    if (y>mid)query(k<<1|1,x,y);    updata(k);}void change(int k,int x,int y,int z){    pushdown(k);    int l=tr[k].l,r=tr[k].r;    if (x<=l&&r<=y)    {        tr[k].mn=tr[k].tag=z;        return;    }    int mid=(l+r)>>1;    if (x<=mid)change(k<<1,x,y,z);    if (y>mid)change(k<<1|1,x,y,z);    updata(k);}void solvechange(int x,int y,int z){    while(tp[x]!=tp[y])    {        if(h[x]<h[y])swap(x,y);        change(1,num[tp[x]],num[x],z);        x=fa[tp[x]][0];    }    if (h[x]>h[y])swap(x,y);    change(1,num[x],num[y],z);}int solvequery(int x){    int z;    qx=inf;    int f=lca(x,root);    if (x==root)qx=tr[1].mn;    else if (f!=x)query(1,num[x],num[x]+size[x]-1);    else if (f==x)    {        z=lca1(x,root);        int l=num[z]-1,r=num[z]+size[z];        if (l>=1)query(1,1,l);        if (r<=n)query(1,r,n);//      cout<<x<<' '<<root<<' '<<z<<endl;    }    return qx;}int main(){//  freopen("std.in","r",stdin);//  freopen("dp.out","w",stdout);    clear();    int u,v,opx,x,y,z;    n=F(),m=F();    for(int i=1;i<n;++i)    u=F(),v=F(),add(u,v),add(v,u);    for(int i=1;i<=n;++i)    val[i]=F();    root=F();    dfs1(1);    dfs2(1,1);    for (int i=1;i<=n;++i)    tree[num[i]]=i;    build(1,1,n);//  for (int i=1;i<=20;++i)//  cout<<tr[i].l<<' '<<tr[i].r<<' '<<tr[i].mn<<' '<<tr[i].tag<<endl;    for (int i=1;i<=m;++i)    {        opx=F();        if (opx==1)root=F();        else if(opx==2)x=F(),y=F(),z=F(),solvechange(x,y,z);        else if(opx==3)x=F(),printf("%d\n",solvequery(x));    }    return 0;}