这个题目名字好大众啊。
从昨天下午开始码,刚调出来。
由于LCT常数很大,一定要丧心病狂的优化常数,非常管用。
#include<iostream>
#include<cstdio>
#define LL unsigned int
#define P 51061
#define maxn 200005
using namespace std;
int n,q;
int fa[maxn],ch[maxn][2],sz[maxn];
LL v[maxn],sum[maxn],add[maxn],mul[maxn];
bool rev[maxn];
bool dir(int x)
{return x==ch[fa[x]][1];}
bool isroot(int x)
{return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x;}
void up(int x)
{
sz[x]=sz[ch[x][0]]+sz[ch[x][1]]+1;
sum[x]=(sum[ch[x][0]]+sum[ch[x][1]]+v[x])%P;
}
void down(int x)
{
if(rev[x])
{
rev[ch[x][0]]^=1;
rev[ch[x][1]]^=1;
rev[x]=0;
swap(ch[x][0],ch[x][1]);
}
if(mul[x]==1&&add[x]==0) return ;
if(ch[x][0])
{
v[ch[x][0]]=(v[ch[x][0]]*mul[x]+add[x])%P;
sum[ch[x][0]]=(sum[ch[x][0]]*mul[x]+add[x]*sz[ch[x][0]])%P;
mul[ch[x][0]]=(mul[x]*mul[ch[x][0]])%P;
add[ch[x][0]]=(add[ch[x][0]]*mul[x]+add[x])%P;
}
if(ch[x][1])
{
v[ch[x][1]]=(v[ch[x][1]]*mul[x]+add[x])%P;
sum[ch[x][1]]=(sum[ch[x][1]]*mul[x]+add[x]*sz[ch[x][1]])%P;
mul[ch[x][1]]=(mul[x]*mul[ch[x][1]])%P;
add[ch[x][1]]=(add[ch[x][1]]*mul[x]+add[x])%P;
}
mul[x]=1;add[x]=0;
}
void rotate(int x)
{
int y=fa[x],z=fa[y];
bool b=dir(x);
int a=ch[x][!b];
if(!isroot(y)) ch[z][dir(y)]=x;
fa[x]=z;ch[x][!b]=y;
fa[y]=x;ch[y][b]=a;
if(a) fa[a]=y;
up(y);up(x);
}
int S[maxn],tot;
void splay(int x)
{
S[++tot]=x;
for(int i=x;!isroot(i);i=fa[i]) //这里迷之错误,忘了加'!' 调了一天
S[++tot]=fa[i];
while(tot)down(S[tot--]);
while(!isroot(x))
{
int y=fa[x];int z=fa[y];
if(isroot(y)) rotate(x);
else
{
bool b=dir(x),c=dir(y);
if(b^c) {rotate(y);rotate(x);}
else {rotate(x);rotate(x);}
}
}
}
void access(int x)
{
for(int t=0;x;t=x,x=fa[x])
splay(x),ch[x][1]=t,up(x);
}
void make_root(int x)
{access(x);splay(x);rev[x]^=1;}
void link(int x,int y)
{make_root(x);fa[x]=y;}
void cut(int x,int y)
{
make_root(x);
access(y);splay(y);
ch[y][0]=fa[x]=0;
up(y);
}
void change(int x,int y,LL a,LL m)
{
make_root(x);
access(y);splay(y);
v[y]=(v[y]*m+a)%P;
sum[y]=(sum[y]*m+sz[y]*a)%P;
mul[y]*=m;mul[y]%=P;
add[y]=(add[y]*m+a)%P;
}
int ask(int x,int y)
{
make_root(x);
access(y);splay(y);
return sum[y]%P;
}
void read(int &a)
{
a=0;char c=getchar();
while(c<'0'||c>'9')
c=getchar();
while(c>='0'&&c<='9')
{
a=a*10+c-'0';
c=getchar();
}
}
char C[3];
int main()
{
read(n);read(q);
for(int i=1;i<=n;i++)
v[i]=sum[i]=mul[i]=sz[i]=1;
int u,v,u1,v1,c;
for(int i=1;i<n;i++)
{
read(u);read(v);
link(v,u);
}
for(int i=1;i<=q;i++)
{
scanf("%s",C);
if(C[0]=='+')
{
read(u);read(v);read(c);
change(u,v,c,1);
}
if(C[0]=='-')
{
read(u);read(v);read(u1);read(v1);
cut(u,v);link(u1,v1);
}
if(C[0]=='*')
{
read(u);read(v);read(c);
change(u,v,0,c);
}
if(C[0]=='/')
{
read(u);read(v);
printf("%d\n",ask(u,v));
}
}
return 0;
}