bzoj1180 [CROATIAN2009]OTOCI

http://www.elijahqi.win/2018/01/28/bzoj1180-croatian2009otoci/ ‎

Description

给出n个结点以及每个点初始时对应的权值wi。起始时点与点之间没有连边。有3类操作： 1、bridge A B：询问结点A与结点B是否连通。如果是则输出“no”。否则输出“yes”，并且在结点A和结点B之间连一条无向边。 2、penguins A X：将结点A对应的权值wA修改为X。 3、excursion A B：如果结点A和结点B不连通，则输出“impossible”。否则输出结点A到结点B的路径上的点对应的权值的和。给出q个操作，要求在线处理所有操作。数据范围：1<=n<=30000, 1<=q<=300000, 0<=wi<=1000。
Input

第一行包含一个整数n(1<=n<=30000)，表示节点的数目。第二行包含n个整数，第i个整数表示第i个节点初始时对应的权值。第三行包含一个整数q(1<=n<=300000)，表示操作的数目。以下q行，每行包含一个操作，操作的类别见题目描述。任意时刻每个节点对应的权值都是1到1000的整数。
Output

输出所有bridge操作和excursion操作对应的输出，每个一行。
Sample Input
5
4 2 4 5 6
10
excursion 1 1
excursion 1 2
bridge 1 2
excursion 1 2
bridge 3 4
bridge 3 5
excursion 4 5
bridge 1 3
excursion 2 4
excursion 2 5
Sample Output
4
impossible
yes
6
yes
yes
15
yes
15
16
这代码我调试了将近一整天 因为写代码的时候不严谨..
注意这题bridge的时候需要输出的是反着的 如果不连通则输出yes 然后把他们再连起来 其他就是正常的lct

#include<cstdio>
#include<algorithm>
#define N 33000
using namespace std;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=getchar();}
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int c[N][2],fa[N],v[N],n,m,top,q[N],sum[N];
bool rev[N];char op[10];
inline bool isroot(int x){return c[fa[x]][0]!=x&&c[fa[x]][1]!=x;}
inline void update(int x){sum[x]=sum[c[x][0]]+sum[c[x][1]]+v[x];}
inline void pushdown(int x){
    int l=c[x][0],r=c[x][1];if (!rev[x]) return;
    rev[l]^=1;rev[r]^=1;rev[x]=0;swap(c[x][0],c[x][1]);
}
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){
    q[top=1]=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;}
inline void print(int x){
    pushdown(x);
    if (c[x][0]) print(c[x][0]);
    printf("%d %d %d %d\n",v[x],x,c[x][0],c[x][1]);
    if (c[x][1]) print(c[x][1]);
}
inline int find(int x){access(x);splay(x);while(c[x][0])x=c[x][0];return x;}
int main(){
    freopen("bzoj1180.in","r",stdin);
    n=read();for (int i=1;i<=n;++i) v[i]=read();m=read();
    for (int i=1;i<=m;++i){
        scanf("%s",op);int x=read(),y=read();
        if(op[0]=='e'){

            if (find(x)!=find(y)) puts("impossible");
            else {
            //  print(y),puts("laji");
            //  for (int i=1;i<=n;++i) printf("%d ",fa[i]);puts("无语");
                makeroot(x);
                access(y);splay(y);printf("%d\n",sum[y]);
            }


        }
        if (op[0]=='b') if (find(x)!=find(y)) puts("yes"),makeroot(x),fa[x]=y;else puts("no");
        if (op[0]=='p') {
            //print(x);puts("zhazha");access(x);print(x);puts("标记");splay(x);print(x);puts("标记");rev[x]^=1; 
            makeroot(x);
            access(x);
            splay(x);v[x]=sum[x]=y;//puts("1");puts("3");print(x);print(x);puts("1");print(x);puts("2");
        }/*for (int i=1;i<=n;++i) printf("%d ",fa[i]);puts("无语");
        for (int i=1;i<=n;++i) printf("%d %d\n",c[i][0],c[i][1]);
        for (int i=1;i<=n;++i) printf("%d ",rev[i]);puts("\\");*/
    }
    return 0;
}