在此为LCT开一个永久的坑

其实我连splay都还不怎么会。

今天先抄了黄学长的bzoj2049，以后一定要把它理解了。

写LCT怎么能不%数据结构大神yeweining呢？%%%chrysanthemums  %%%切掉大森林的全ZJ唯一一人

#include<cstdio>
#include<algorithm>
#define N 1000005
using namespace std;
int n,m;
],c[][],f[][],st[];];
inline bool isroot(int x)
{
    ]!=x&&c[fa[x]][]!=x;
}
void push_down(int k)
{
    ],r=c[k][];
    if(rev[k])
    {
      rev[k]^=;rev[l]^=;rev[r]^=;
      swap(c[k][],c[k][]);
    }
}
void rotate(int x)
{
   int y=fa[x],z=fa[y],l,r;
   ]==x)l=;;
   r=l^;
   )
   {
       ]==y)c[z][]=x;]=x;
   }
   fa[x]=z;fa[y]=x;fa[c[x][r]]=y;
   c[y][l]=c[x][r];c[x][r]=y;
}
void splay(int x)
{
   ;st[++top]=x;
   for(int i=x;!isroot(i);i=fa[i])
   {
     st[++top]=fa[i];
   }
   for(int i=top;i;i--)push_down(st[i]);
   while(!isroot(x))
   {
      int y=fa[x],z=fa[y];
      if(!isroot(y))
      {
         ]==x^c[z][]==y)rotate(x);
         else rotate(y);
      }
      rotate(x);
   }
}
void access(int x)
{
  ;
  while(x)
  {
     splay(x);c[x][]=t;t=x;x=fa[x];
  }
}
void rever(int x)
{
    access(x);splay(x);rev[x]^=;
}
void link(int x,int y)
{
    rever(x);fa[x]=y;splay(x);
}
void cut(int x,int y)
{
    rever(x);access(y);splay(y);c[y][]=fa[x]=;
}
int find(int x)
{
    access(x);splay(x);int y=x;
    ])y=c[y][];
    return y;
}
int main()
{
    ];
    scanf("%d%d",&n,&m);
    ;i<=m;i++)
    {
        scanf("%s",s);int x,y;
        scanf("%d%d",&x,&y);
        ]=='C')link(x,y);
        ]=='D')cut(x,y);
        else
        {
            if(find(x)==find(y))printf("Yes\n");else printf("No\n");
        }
    }
 } 

洞穴勘测

http://wenku.baidu.com/link?url=x6FCcjtWt-mcj9-do9MFL7lsXB4o2G1T0aG3p-XFYkuNyLKW0Kg-o-u_DUWnI_L_yIVJfHcrE_ThkrfoBlmWqJM-TburN1KlZO9wOvKJmA3

---------这是一篇很详尽的讲稿