bzoj2631 tree lct

http://www.elijahqi.win/2018/01/28/bzoj2631-tree/ ‎
Description

　一棵n个点的树，每个点的初始权值为1。对于这棵树有q个操作，每个操作为以下四种操作之一：
+ u v c：将u到v的路径上的点的权值都加上自然数c；
- u1 v1 u2 v2：将树中原有的边(u1,v1)删除，加入一条新边(u2,v2)，保证操作完之后仍然是一棵树；
* u v c：将u到v的路径上的点的权值都乘上自然数c；
/ u v：询问u到v的路径上的点的权值和，求出答案对于51061的余数。

Input
　　第一行两个整数n，q
接下来n-1行每行两个正整数u，v，描述这棵树
接下来q行，每行描述一个操作

Output
　　对于每个/对应的答案输出一行

Sample Input
3 2
1 2
2 3
* 1 3 4
/ 1 1

Sample Output

4

HINT

数据规模和约定

10%的数据保证，1<=n，q<=2000

另外15%的数据保证，1<=n，q<=5*10^4，没有-操作，并且初始树为一条链

另外35%的数据保证，1<=n，q<=5*10^4，没有-操作

100%的数据保证，1<=n，q<=10^5，0<=c<=10^4

正常的一个lct但是凭空多出了乘法操作就比较难过了 但是其实可以参考洛谷线段树模板2 和这个就非常相似了 处理好那个再来处理这个就比较轻松了 这题有个坑点这个模数的平方会爆int 那我比较懒 用惯了unsigned int就直接上了 据说开long long 会tle 那蒟蒻我 大常数的lct就不要轻易尝试了 一开始写的时候差一点犯错误 就是我想一直让它自然溢出 然后最后再取模 这样显然答案是不对的 具体可以参考noi2012随机数生成器 里面的一个过程同理 有了这个乘法操作怎么搞呢 我每次下放的时候先乘后加即可 对于下放的时候维护的东西也需要注意一下 比如我现在已经有个加法标记了 这时候我再来要给乘法 那么显然我为了避免麻烦我就 直接给这个加法标记乘相同的倍数 后面累加的时候便可以直接累加了 剩下还是先乘后加 搞定

#include<cstdio>
#include<algorithm>
#define mod 51061
#define N 110000
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=gc();}
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x*f;
}
struct node{
    int y,z,next;
}data[N<<1];
unsigned int sum[N],v[N],add[N],mul[N];
int q[N],top,fa[N],h[N],n,qq,num,size[N],c[N][2],rev[N];
inline void update(int x){
    size[x]=size[c[x][0]]+size[c[x][1]]+1;
    sum[x]=sum[c[x][0]]+v[x]+sum[c[x][1]];sum[x]%=mod;
}
inline void push(int x,unsigned int ad,int mu,int re){
    v[x]=(v[x]*mu+ad)%mod;
    sum[x]=(sum[x]*mu+ad*size[x])%mod;
    (mul[x]*=mu)%=mod;add[x]=(add[x]*mu+ad)%mod;
    if (re) rev[x]^=1,swap(c[x][0],c[x][1]);
}
inline void pushdown(int x){
    push(c[x][0],add[x],mul[x],rev[x]);push(c[x][1],add[x],mul[x],rev[x]);
    rev[x]=add[x]=0;mul[x]=1;
}
inline bool isroot(int x){return c[fa[x]][0]!=x&&c[fa[x]][1]!=x;}
inline void rotate(int x){
    int y=fa[x],z=fa[y];
    if (!isroot(y)) c[z][c[z][1]==y]=x;
    int l=c[y][1]==x,r=l^1;
    fa[c[x][r]]=y;fa[y]=x;fa[x]=z;
    c[y][l]=c[x][r];c[x][r]=y;update(y);update(x);
}
inline void splay(int x){
    top=1;q[top]=x;for (int i=x;!isroot(i);i=fa[i]) q[++top]=fa[i];
    while(top) pushdown(q[top--]);
    while(!isroot(x)){
        int y=fa[x],z=fa[y];
        if (!isroot(y)){
            if (c[y][0]==x^c[z][0]==y) rotate(x);else rotate(y);
        }rotate(x);
    }
}
inline void access(int x){
    for (int t=0;x;t=x,x=fa[x]) {splay(x);
        c[x][1]=t;
        update(x);}}
inline void makeroot(int x){access(x);splay(x);rev[x]^=1;swap(c[x][0],c[x][1]);}
inline void print(int x){
    pushdown(x);
    if (c[x][0]) print(c[x][0]);
    printf("%u %u %d %d %d %d\n",v[x],sum[x],size[x],x,c[x][0],c[x][1]);
    if (c[x][1]) print(c[x][1]);
}
int main(){
    freopen("bzoj2631.in","r",stdin);
    n=read();qq=read();
    for (int i=1;i<n;++i) {
        int x=read(),y=read();
        data[++num].y=y;data[num].next=h[x];h[x]=num;
        data[++num].y=x;data[num].next=h[y];h[y]=num;
    }int l=1,r=0;q[++r]=1;
    while(l<=r){
        int x=q[l++];
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y;if (fa[x]==y) continue;fa[y]=x;q[++r]=y;
        }
    }for (int i=1;i<=n;++i) size[i]=1,mul[i]=1,v[i]=1;
    for (int i=1;i<=qq;++i){
        char ch=gc();while(ch!='+'&&ch!='-'&&ch!='*'&&ch!='/') ch=gc();
        if (ch=='+'){
            int x=read(),y=read(),c=read();
            makeroot(x);access(y);splay(y);push(y,c,1,0);
        }
        if (ch=='-'){
            int u1=read(),v1=read(),u2=read(),v2=read();
            makeroot(u1);access(v1);splay(v1);c[v1][0]=fa[u1]=0;update(v1);
            makeroot(u2);fa[u2]=v2;
        }
        if (ch=='*'){
            int x=read(),y=read(),c=read();
            makeroot(x);access(y);splay(y);push(y,0,c,0);//print(x);puts("");
        }
        if (ch=='/'){
            int x=read(),y=read();makeroot(x);access(y);splay(y);printf("%u\n",sum[y]);
        }
    }
    return 0;
}