P2147 [SDOI2008]洞穴勘测(LCT维护树的连通性)

题目描述

辉辉热衷于洞穴勘测。

某天,他按照地图来到了一片被标记为JSZX的洞穴群地区。经过初步勘测,辉辉发现这片区域由n个洞穴(分别编号为1到n)以及若干通道组成,并且每条通道连接了恰好两个洞穴。假如两个洞穴可以通过一条或者多条通道按一定顺序连接起来,那么这两个洞穴就是连通的,按顺序连接在一起的这些通道则被称之为这两个洞穴之间的一条路径。 洞穴都十分坚固无法破坏,然而通道不太稳定,时常因为外界影响而发生改变,比如,根据有关仪器的监测结果,123号洞穴和127号洞穴之间有时会出现一条通道,有时这条通道又会因为某种稀奇古怪的原因被毁。

辉辉有一台监测仪器可以实时将通道的每一次改变状况在辉辉手边的终端机上显示:

如果监测到洞穴u和洞穴v之间出现了一条通道,终端机上会显示一条指令 Connect u v

如果监测到洞穴u和洞穴v之间的通道被毁,终端机上会显示一条指令 Destroy u v

经过长期的艰苦卓绝的手工推算,辉辉发现一个奇怪的现象:无论通道怎么改变,任意时刻任意两个洞穴之间至多只有一条路径。

因而,辉辉坚信这是由于某种本质规律的支配导致的。因而,辉辉更加夜以继日地坚守在终端机之前,试图通过通道的改变情况来研究这条本质规律。 然而,终于有一天,辉辉在堆积成山的演算纸中崩溃了……他把终端机往地面一砸(终端机也足够坚固无法破坏),转而求助于你,说道:“你老兄把这程序写写吧”。

辉辉希望能随时通过终端机发出指令 Query u v,向监测仪询问此时洞穴u和洞穴v是否连通。现在你要为他编写程序回答每一次询问。 已知在第一条指令显示之前,JSZX洞穴群中没有任何通道存在。

输入格式

第一行为两个正整数n和m,分别表示洞穴的个数和终端机上出现过的指令的个数。 以下m行,依次表示终端机上出现的各条指令。每行开头是一个表示指令种类的字符串s("Connect”、”Destroy”或者”Query”,区分大小写),之后有两个整数u和v (1≤u, v≤n) 分别表示两个洞穴的编号。

输出格式

对每个Query指令,输出洞穴u和洞穴v是否互相连通:是输出”Yes”,否则输出”No”。(不含双引号)

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e6+100;
const int mod=51061;
typedef long long ll;
ll f[maxn];
ll ch[maxn][2];
ll v[maxn];
ll size[maxn];
ll lm[maxn];
ll la[maxn];
ll s[maxn];
ll st[maxn];
bool r[maxn];
bool nroot (int x) {
    return ch[f[x]][0]==x||ch[f[x]][1]==x;
}
void pushup (int x) {
    s[x]=s[ch[x][0]]+s[ch[x][1]]+v[x];
    s[x]%=mod;
    
    size[x]=size[ch[x][0]]+size[ch[x][1]]+1;
}

void rev (int x) {
    //翻转
    int t=ch[x][0];
    ch[x][0]=ch[x][1];
    ch[x][1]=t;
    r[x]^=1; 
}

void push_mul (int x,int c) {
    //
    s[x]*=c;s[x]%=mod;
    v[x]*=c;v[x]%=mod;
    lm[x]*=c;lm[x]%=mod;
    la[x]*=c;la[x]%=mod;
}

void push_add (int x,int c) {
    //
    s[x]+=c*size[x];s[x]%=mod;
    v[x]+=c;v[x]%=mod;
    la[x]+=c;la[x]%=mod;
}

void pushdown (int x) {
    if (lm[x]!=1) push_mul(ch[x][0],lm[x]),push_mul(ch[x][1],lm[x]),lm[x]=1;
    if (la[x]) push_add(ch[x][0],la[x]),push_add(ch[x][1],la[x]),la[x]=0;
    if (r[x]) {
        if (ch[x][0]) rev(ch[x][0]);
        if (ch[x][1]) rev(ch[x][1]);
        r[x]=0;
    } 
}

void rotate (int x) {
    int y=f[x],z=f[y],k=ch[y][1]==x,w=ch[x][!k];
    if (nroot(y)) ch[z][ch[z][1]==y]=x;
    ch[x][!k]=y;
    ch[y][k]=w;
    if (w) f[w]=y;
    f[y]=x;f[x]=z;
    pushup(y);
}

void splay (int x) {
    int y=x,z=0;
    st[++z]=y;
    while (nroot(y)) st[++z]=y=f[y];
    while (z) pushdown(st[z--]);
    while (nroot(x)) {
        y=f[x];z=f[y];
        if (nroot(y)) rotate((ch[y][0]==x)^(ch[z][0]==y)?x:y);
        rotate(x);
    }
    pushup(x);
}

void access (int x) {
    for (int y=0;x;x=f[y=x]) splay(x),ch[x][1]=y,pushup(x);
}

void makeroot(int x) {
    access(x);
    splay(x);
    rev(x);
}

void split (int x,int y) {
    makeroot(x);
    access(y);
    splay(y);
}

void link (int x,int y) {
    makeroot(x);
    f[x]=y;
}

void cut (int x,int y) {
    split(x,y);
    f[x]=ch[y][0]=0;
}

int find (int x) {
    //找树根
    access(x);
    splay(x);
    while (ch[x][0]) x=ch[x][0];
    return x; 
}
int n,q;
int main () {
    scanf("%d%d",&n,&q);
    while (q--) {
        string op;
        cin>>op;
        if (op=="Connect") {
            int x,y;
            scanf("%d%d",&x,&y);
            link(x,y);
        }
        else if (op=="Destroy") {
            int x,y;
            scanf("%d%d",&x,&y);
            cut(x,y);
        }
        else if (op=="Query") {
            int x,y;
            scanf("%d%d",&x,&y);
            if (find(x)==find(y))
                printf("Yes
");
            else
                printf("No
");
        }
    }
}
原文地址:https://www.cnblogs.com/zhanglichen/p/13770376.html