BZOJ4530:[BJOI2014]大融合(LCT)

时间:2023-03-09 03:04:42
BZOJ4530:[BJOI2014]大融合(LCT)

Description

小强要在N个孤立的星球上建立起一套通信系统。这套通信系统就是连接N个点的一个树。
这个树的边是一条一条添加上去的。在某个时刻,一条边的负载就是它所在的当前能够
联通的树上路过它的简单路径的数量。
BZOJ4530:[BJOI2014]大融合(LCT)
例如,在上图中,现在一共有了5条边。其中,(3,8)这条边的负载是6,因
为有六条简单路径2-3-8,2-3-8-7,3-8,3-8-7,4-3-8,4-3-8-7路过了(3,8)。
现在,你的任务就是随着边的添加,动态的回答小强对于某些边的负载的
询问。

Input

第一行包含两个整数N,Q,表示星球的数量和操作的数量。星球从1开始编号。
接下来的Q行,每行是如下两种格式之一:
A x y 表示在x和y之间连一条边。保证之前x和y是不联通的。
Q x y 表示询问(x,y)这条边上的负载。保证x和y之间有一条边。
1≤N,Q≤100000

Output

对每个查询操作,输出被查询的边的负载。

Sample Input

8 6
A 2 3
A 3 4
A 3 8
A 8 7
A 6 5
Q 3 8

Sample Output

6

Solution

一个讲lct维护子树讲的很好的blog
感觉很好理解,就是$pushup$的时候把虚树的也合并一下
就是很鸡肋,因为这样就不能维护边了
$Size$表示子树的全部信息,$Si$表示虚子树的信息
询问的时候就是边两个端点的子树和的乘积
此题$link$的时候要将$y~makeroot$不然信息就会错误(在这挂了好久

Code

 #include<iostream>
#include<cstring>
#include<cstdio>
#define N (100000+100)
using namespace std; int Father[N],Son[N][],Size[N],Si[N],Rev[N];
int n,m,x,y;
char opt[]; int Get(int x) {return Son[Father[x]][]==x;}
void Update(int x) {Size[x]=Si[x]+Size[Son[x][]]+Size[Son[x][]]+;}
int Is_root(int x) {return Son[Father[x]][]!=x && Son[Father[x]][]!=x;} void Rotate(int x)
{
int wh=Get(x);
int fa=Father[x],fafa=Father[fa];
if (!Is_root(fa)) Son[fafa][Son[fafa][]==fa]=x;
Father[fa]=x; Son[fa][wh]=Son[x][wh^];
if (Son[fa][wh]) Father[Son[fa][wh]]=fa;
Father[x]=fafa; Son[x][wh^]=fa;
Update(fa); Update(x);
} void Pushdown(int x)
{
if (Rev[x] && x)
{
if (Son[x][]) Rev[Son[x][]]^=;
if (Son[x][]) Rev[Son[x][]]^=;
swap(Son[x][],Son[x][]);
Rev[x]=;
}
} void Push(int x) {if (!Is_root(x)) Push(Father[x]); Pushdown(x);}
void Splay(int x)
{
for (int fa;!Is_root(x);Rotate(x))
if (!Is_root(fa=Father[x]))
Rotate(Get(fa)==Get(x)?fa:x);
} void Access(int x) {for (int y=;x;y=x,x=Father[x]) Splay(x), Si[x]+=Size[Son[x][]], Si[x]-=Size[y], Son[x][]=y, Update(x);}
void Make_root(int x) {Access(x); Splay(x); Rev[x]^=;}
int Find_root(int x) {Access(x); Splay(x); while (Son[x][]) x=Son[x][]; return x;}
void Link(int x,int y) {Make_root(x); Make_root(y); Father[x]=y; Si[y]+=Size[x]; Update(y);}
int Query(int x,int y)
{
Make_root(x);
Access(y);
Splay(y);
int sum=Size[y];
int sizex=Size[x];
int sizey=sum-sizex;
return sizex*sizey;
} int main()
{
scanf("%d%d",&n,&m);
for(int i=;i<=n;++i) Size[i]=;
for (int i=;i<=m;++i)
{
scanf("%s%d%d",opt,&x,&y);
if (opt[]=='A') Link(x,y);
else printf("%d\n",Query(x,y));
}
}