城市旅行

一、题目

点此看题

二、解法

考虑用$\text{ lct }$维护这个期望，我们先把柿子展开（分子考虑每个点的贡献，分母用等差数列求和）：
$$exp=\frac{a_1\times 1\times n+a_2\times 2\times (n-1)+....+a_n\times n\times 1}{\frac{n\times (n+1)}{2}}$$我们维护两个值，$ls=a_1+2a_2+...+n{a}_n,rs=n{a}_1+(n-1)a_2+...+a_n$，那么根节点期望也可以维护：
$$ex{p}_x=ex{p}_l+ex{p}_r+va{l}_x\times (si{z}_l+1)\times (si{z}_r+1)+l{s}_l\times (si{z}_r+1)+r{s}_r\times (si{z}_l+1)$$$ls,rs$也比较容易维护：
$$ls=l{s}_l+l{s}_r+su{m}_r\times (si{z}_l+1)+va{l}_x\times (si{z}_l+1)$$$$rs=r{s}_l+r{s}_r+su{m}_l\times (si{z}_r+1)+va{l}_x\times (si{z}_r+1)$$那么我们就完成了一半的工作（$push\_up$）

---

考虑修改，$val,sum,lazy$都很好维护，$ls$的维护也不难，用等差数列（$rs$也同理）：
$$ls=ls+p\times \frac{si{z}_x(si{z}_x+1)}{2}$$$exp$就有点**了，我们设$n=si{z}_x$，一开始有这个柿子：
$$exp=exp+p\times (1\times n+2\times (n-1)+...+n\times 1)$$问题变成了后面柿子的化简，网上大多是推导，我讲一种无脑方法，就是用拉格朗日插值（我的博客，看公式就行了），先盲猜它通项是一个三次多项式，我们选取四个点插它：$(0,0),(1,1),(2,4),(3,10)$，然后带入公式里化简，过程略，给出结果（因式分解后）：
$$\frac{n(n+1)(n+2)}{6}$$然后$exp$也可以维护了，$\text{lct}$ 打标记即可解决这个问题，总结一些坑点：

• $flip$时，不仅左右儿子要交换，$ls,rs$也要交换。
• 操作三跳过是要把第三个数个读进来（虽然用不到）。

#include <cstdio>
#include <iostream>
#define lc ch[x][0]
#define rc ch[x][1]
#define int long long
using namespace std;
const int M = 50005;
int read()
{
	int x=0,flag=1;char c;
	while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
	while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return x*flag;
}
int n,m,par[M],siz[M],ch[M][2],val[M],exp[M];
int fl[M],st[M],ls[M],rs[M],la[M],sum[M];
int nrt(int x)
{
	return ch[par[x]][0]==x || ch[par[x]][1]==x;
}
int chk(int x)
{
	return ch[par[x]][1]==x;
}
void push_up(int x)
{
	if(!x) return ;
	siz[x]=siz[ch[x][0]]+siz[ch[x][1]]+1;
	sum[x]=sum[ch[x][0]]+sum[ch[x][1]]+val[x];
	ls[x]=ls[lc]+ls[rc]+sum[rc]*(siz[lc]+1)+val[x]*(siz[lc]+1);
	rs[x]=rs[lc]+rs[rc]+sum[lc]*(siz[rc]+1)+val[x]*(siz[rc]+1);
	exp[x]=exp[lc]+exp[rc]+val[x]*(siz[lc]+1)*(siz[rc]+1)+ls[lc]*(siz[rc]+1)+rs[rc]*(siz[lc]+1);
}
void flip(int x)
{
	if(!x) return ;
	swap(ch[x][0],ch[x][1]);
	swap(ls[x],rs[x]);
	fl[x]^=1;
}
void add(int x,int p)
{
	if(!x) return ;
	val[x]+=p;la[x]+=p;
	sum[x]+=p*siz[x];
	ls[x]+=p*siz[x]*(siz[x]+1)/2;
	rs[x]+=p*siz[x]*(siz[x]+1)/2;
	exp[x]+=p*siz[x]*(siz[x]+1)*(siz[x]+2)/6;
}
void push_down(int x)
{
	if(!x) return ;
	if(fl[x])
	{
		flip(ch[x][0]);flip(ch[x][1]);
		fl[x]=0;
	}
	if(la[x])
	{
		add(ch[x][0],la[x]);
		add(ch[x][1],la[x]);
		la[x]=0;
	}
}
void rotate(int x)
{
	int y=par[x],z=par[y],k=chk(x),w=ch[x][k^1];
	ch[y][k]=w;par[w]=y;
	if(nrt(y)) ch[z][chk(y)]=x;par[x]=z;
	ch[x][k^1]=y;par[y]=x;
	push_up(y);push_up(x);
}
void splay(int x)
{
	int y=x,z=0;
	st[++z]=y;
	while(nrt(y)) st[++z]=y=par[y];
	while(z) push_down(st[z--]);
	while(nrt(x))
	{
		int y=par[x],z=par[y];
		if(nrt(y))
		{
			if(chk(y)==chk(x)) rotate(y);
			else rotate(x);
		}
		rotate(x);
	}
}
void access(int x)
{
	for(int y=0;x;x=par[y=x])
		splay(x),ch[x][1]=y,push_up(x); 
}
void makeroot(int x)
{
	access(x);splay(x);
	flip(x);
}
int findroot(int x)
{
	access(x);splay(x);
	while(ch[x][0]) push_down(x),x=ch[x][0];
	splay(x);
	return x;
}
void split(int x,int y)
{
	makeroot(x);
	access(y);splay(y);
}
void link(int x,int y)
{
	makeroot(x);
	if(findroot(y)!=x) par[x]=y;
}
void cut(int x,int y)
{
	makeroot(x);
	if(findroot(y)==x && par[y]==x && !ch[y][0])
	{
		par[y]=ch[x][1]=0;
		push_up(x); 
	}
}
int gcd(int a,int b)
{
	return !b?a:gcd(b,a%b);
}
signed main()
{
	n=read();m=read();
	for(int i=1;i<=n;i++)
	{
		val[i]=read();
		push_up(i);
	}
	for(int i=1;i<n;i++)
		link(read(),read());
	while(m--)
	{
		int op=read(),u=read(),v=read();
		if(op==1)
			cut(u,v);
		if(op==2)
			link(u,v);
		if(op==3)
		{
			if(findroot(u)^findroot(v))
			{
				read();//一定要读。。 
				continue ;
			}
			split(u,v);
			add(v,read());
		}
		if(op==4)
		{
			if(findroot(u)^findroot(v))
			{
				puts("-1");
				continue ;
			}
			split(u,v);
			int a=exp[v],b=siz[v]*(siz[v]+1)/2,t=gcd(a,b);
			a/=t;b/=t;
			printf("%lld/%lld\n",a,b);
		}
	}
}