[bzoj3083]遥远的国度

3083: 遥远的国度

Time Limit: 10 Sec Memory Limit: 1280 MB
Submit: 1733 Solved: 429
[Submit][Status][Discuss]
Description

描述
zcwwzdjn在追杀十分sb的zhx，而zhx逃入了一个遥远的国度。当zcwwzdjn准备进入遥远的国度继续追杀时，守护神RapiD阻拦了zcwwzdjn的去路，他需要zcwwzdjn完成任务后才能进入遥远的国度继续追杀。

问题是这样的：遥远的国度有n个城市，这些城市之间由一些路连接且这些城市构成了一颗树。这个国度有一个首都，我们可以把这个首都看做整棵树的根，但遥远的国度比较奇怪，首都是随时有可能变为另外一个城市的。遥远的国度的每个城市有一个防御值，有些时候RapiD会使得某两个城市之间的路径上的所有城市的防御值都变为某个值。RapiD想知道在某个时候，如果把首都看做整棵树的根的话，那么以某个城市为根的子树的所有城市的防御值最小是多少。由于RapiD无法解决这个问题，所以他拦住了zcwwzdjn希望他能帮忙。但zcwwzdjn还要追杀sb的zhx，所以这个重大的问题就被转交到了你的手上。

Input

第1行两个整数n m，代表城市个数和操作数。
第2行至第n行，每行两个整数 u v，代表城市u和城市v之间有一条路。
第n+1行，有n个整数，代表所有点的初始防御值。
第n+2行一个整数 id，代表初始的首都为id。
第n+3行至第n+m+2行，首先有一个整数opt，如果opt=1，接下来有一个整数id，代表把首都修改为id；如果opt=2，接下来有三个整数p1 p2 v，代表将p1 p2路径上的所有城市的防御值修改为v；如果opt=3，接下来有一个整数 id，代表询问以城市id为根的子树中的最小防御值。

Output

对于每个opt=3的操作，输出一行代表对应子树的最小点权值。

Sample Input

3 7

1 2

1 3

1 2 3

1

3 1

2 1 1 6

3 1

2 2 2 5

3 1

2 3 3 4

3 1

Sample Output

1

2

3

4

提示

对于20%的数据，n<=1000 m<=1000。

对于另外10%的数据，n<=100000，m<=100000，保证修改为单点修改。

对于另外10%的数据，n<=100000，m<=100000，保证树为一条链。

对于另外10%的数据，n<=100000，m<=100000，没有修改首都的操作。

对于100%的数据，n<=100000，m<=100000，0<所有权值<=2^31。

题解：这个题主要就是有个换根操作。
假设当前的根式root，要查询的是以x为子树的最小值，那么可以分成三种情况：
①：x==root 那么就是查询整棵树就行了。
②：LCA(x,root)!=x 那么查询x在原来的树中的子树就行了，换根对答案不会有影响。
③：LCA(x,root)==x 这种情况就是在查询的时候扣掉root所在的x的那棵子树就行了。

其他的就是裸的链剖了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define inf 0x7fffffff
const int N=100010;
struct S{int st,en;}aa[N*2];
int n,m,tot,point[N],next[N*2],a[N],root,fa[N][20];
int deep[N],pos[N],num,belong[N],s[N],tr[N*4],de[N*4],siz[N];
inline int in(){
    int x=0; char ch=getchar();
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x;
}
inline void add(int x,int y){
    tot+=1;next[tot]=point[x];point[x]=tot;
    aa[tot].st=x;aa[tot].en=y;
    tot+=1;next[tot]=point[y];point[y]=tot;
    aa[tot].st=y;aa[tot].en=x;
}
inline void dfs_1(int x,int last){
    int i;
    siz[x]=1;
    for(i=1;i<=16;++i){
        if(deep[x]<(1<<i)) break;
        fa[x][i]=fa[fa[x][i-1]][i-1];
    }
    for(i=point[x];i;i=next[i])
      if(aa[i].en!=last){
        deep[aa[i].en]=deep[x]+1;
        fa[aa[i].en][0]=x;
        dfs_1(aa[i].en,x);
        siz[x]+=siz[aa[i].en];
      }
}
inline void dfs_2(int x,int y){
    int i,k=0;
    num+=1;
    belong[x]=y;
    pos[x]=s[x]=num;
    for(i=point[x];i;i=next[i])
      if(deep[aa[i].en]>deep[x]&&siz[k]<siz[aa[i].en])
        k=aa[i].en;
    if(!k) return ;
    dfs_2(k,y);
    s[x]=max(s[x],s[k]);
    for(i=point[x];i;i=next[i])
      if(deep[aa[i].en]>deep[x]&&k!=aa[i].en){
        dfs_2(aa[i].en,aa[i].en);
        s[x]=max(s[x],s[aa[i].en]);
      }
}
inline int LCA(int x,int y){
    if(deep[x]<deep[y]) swap(x,y);
    int t=deep[x]-deep[y],i;
    for(i=0;i<=16;++i)
      if(t&(1<<i)) x=fa[x][i];
    for(i=16;~i;--i)
      if(fa[x][i]!=fa[y][i])
        x=fa[x][i],y=fa[y][i];
    return x==y?x:fa[x][0];
}
#define mid (l+r)/2
#define L k<<1,l,mid
#define R k<<1|1,mid+1,r
inline void paint(int k,int z){tr[k]=de[k]=z;}
inline void pushdown(int k){
    paint(k<<1,de[k]);
    paint(k<<1|1,de[k]);
    de[k]=0;
}
inline void insert(int k,int l,int r,int x,int y,int z){
    if(x<=l&&y>=r){
        paint(k,z);
        return ;
    }
    if(de[k]) pushdown(k);
    if(x<=mid) insert(L,x,y,z);
    if(y>mid) insert(R,x,y,z);
    tr[k]=min(tr[k<<1],tr[k<<1|1]);
}
inline int query(int k,int l,int r,int x,int y){
    int minn=inf;
    if(x<=l&&y>=r) return tr[k];
    if(de[k]) pushdown(k);
    if(x<=mid) minn=min(minn,query(L,x,y));
    if(y>mid) minn=min(minn,query(R,x,y));
    return minn;
}
inline void change(int x,int y,int z){
    while(belong[x]!=belong[y]){
        insert(1,1,n,pos[belong[x]],pos[x],z);
        x=fa[belong[x]][0];
    }
    insert(1,1,n,pos[y],pos[x],z);
    return ;
}
inline int get_point(int x,int y){
    int t=deep[x]-deep[y]-1,i;
    for(i=0;i<=16;++i)
      if(t&(1<<i)) x=fa[x][i];
    return x; 
}
int main(){
    int i,j,x,y,z,t;
    n=in();m=in();
    for(i=1;i<n;++i){
        x=in();y=in();
        add(x,y);
    }
    for(i=1;i<=n;++i) a[i]=in();
    root=in();
    dfs_1(1,0);
    dfs_2(1,1);
    memset(tr,127,sizeof(tr));
    for(i=1;i<=n;++i) insert(1,1,n,pos[i],pos[i],a[i]);
    while(m--){
        t=in();
        if(t==1) root=in();
        if(t==2){
            x=in();y=in();z=in();
            int lca=LCA(x,y);
            change(x,lca,z); change(y,lca,z);
        }
        if(t==3){
            x=in();
            int lca=LCA(x,root);
            if(x==root) printf("%d\n",tr[1]);
            else{
                if(x!=lca) printf("%d\n",query(1,1,n,pos[x],s[x]));
                else{
                    int r=get_point(root,x);
                    printf("%d\n",min(query(1,1,n,1,pos[r]-1),query(1,1,n,s[r]+1,n)));
                }
            }
        }
    }
}