hdu3726 Graph and Queries

hdu3726
题意呢就是说，给你一个n个点m条边的..东西？然后有三种操作，第一个删除第x条边，第二个查询与x联通的第k大的数，第三个把第x点的值改为v
看到第二个和第三个..能知道用并查集维护会很开心，然后并查集并不能支持删除（我不会..不过好像也可以？？？），就又又又不是很开心了。不过我们再看数据，发现只有删除而没有插入，所以有一个很妙的方法——离散做！倒着做！
ps：预处理好麻烦诶，写到没脾气 而且…变量太多..导致..数组开小了再见再见，改了好久
哦这题还很坑…int会爆…

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 22000
int T=0,n,m,num,cnt,q,nn,n1,h[N],h1[N],now[N],ch[N][2],fa[N],size[N],v[N],f[N],root[N],data[N];
double ans=0;
bool flag[3*N];
char str[5];
struct node1{int x,next;}mp[11*N];
struct node2{int x,y;}a[3*N];
struct node3{int x,y;char ch;}query1[20*N];
int find(int x){
    if(f[x]!=x) f[x]=find(f[x]);
    return f[x];
}
void update(int x){size[x]=size[ch[x][0]]+size[ch[x][1]]+1;}
void rotate(int &x,int rt){
    int y=fa[x],z=fa[y],t=ch[y][0]==x;
    if(z) ch[z][ch[z][1]==y]=x;
    else root[rt]=x;
    fa[y]=x;fa[x]=z;fa[ch[x][t]]=y;
    ch[y][t^1]=ch[x][t];ch[x][t]=y;
    update(y);update(x);
}
void splay(int x,int &rt,int rt1){
    while(x!=rt){
        int y=fa[x];
        if(y==rt){rotate(x,rt1);return;}
        if(ch[y][0]==x^ch[fa[y]][0]==y) rotate(x,rt1);
        else rotate(y,rt1);rotate(x,rt1);
    }
}
void insert(int &p,int x,int father,int rt){
    if(!p){
        p=x;size[p]=1;ch[p][0]=ch[p][1]=0;fa[p]=father;
        splay(p,root[rt],rt);return;
    }if(v[x]>=v[p]) insert(ch[p][0],x,p,rt);
    else insert(ch[p][1],x,p,rt);
}
void merge(int xx,int yy){
    int x=find(xx),y=find(yy);
    if(x==y) return;
    if(size[root[x]]<size[root[y]]) swap(x,y);//y往x中插 
    int op=0,cl=1;data[1]=root[y];
    while(op<cl){
        int p=data[++op];
        if(ch[p][0]) data[++cl]=ch[p][0];
        if(ch[p][1]) data[++cl]=ch[p][1];
        insert(root[x],p,0,x);
    }f[y]=x;
}
void change(int x,int y){   //把x的值变成y 
    int xx=find(x);
    if(size[root[xx]]==1){v[x]=y;return;}
    while(ch[x][0]*ch[x][1]) splay(ch[x][0],x,xx);
    while(ch[x][0]) splay(ch[x][0],x,xx);
    while(ch[x][1]) splay(ch[x][1],x,xx);
    ch[fa[x]][ch[fa[x]][1]==x]=0;
    size[fa[x]]--;splay(fa[x],root[xx],xx);
    v[x]=y;insert(root[xx],x,0,xx);
}
int query(int p,int x){
    if(x<=0 || size[p]<x) return 0;
    if(x<=size[ch[p][0]]) return query(ch[p][0],x);
    if(x==size[ch[p][0]]+1) return v[p];
    return query(ch[p][1],x-size[ch[p][0]]-1);
}
int main(){
    freopen("hdu3726.in","r",stdin);
    while(1){
        scanf("%d%d",&n,&m);
        if(!n && !m) break;T++;
        num=0;cnt=n;q=0;nn=0;ans=0;n1=0;
        memset(h1,0,sizeof(h1));
        memset(h,0,sizeof(h));
        for(int i=1;i<=n;i++){
            int x;scanf("%d",&x);
            mp[i].x=x;mp[i].next=h[i];h[i]=i;now[i]=i;
        }
        memset(flag,true,sizeof(flag));
        for(int i=1;i<=m;i++) scanf("%d%d",&a[i].x,&a[i].y);
        while(1){
            scanf("%s",str);if(str[0]=='E') break;query1[++q].ch=str[0];
            if(str[0]=='Q') scanf("%d%d",&query1[q].x,&query1[q].y);
            if(str[0]=='D'){scanf("%d",&query1[q].x);flag[query1[q].x]=false;}
            if(str[0]=='C'){
                int x,y;scanf("%d%d",&x,&y);
                query1[q].x=x;query1[q].y=y;
                mp[++cnt].x=y;mp[cnt].next=h[x];h[x]=cnt;now[x]=cnt;
            }
        }
        for(int i=1;i<=n;i++) f[i]=root[i]=i,size[i]=1,v[i]=mp[now[i]].x,ch[i][0]=ch[i][1]=fa[i]=0;
        for(int i=1;i<=m;i++)
            if(flag[i]) merge(a[i].x,a[i].y);
        for(int i=q;i>0;i--){
            int x=query1[i].x,y=query1[i].y;
            if(query1[i].ch=='D') merge(a[x].x,a[x].y);
            else if(query1[i].ch=='C'){
                now[x]=mp[now[x]].next;
                change(x,mp[now[x]].x);
            }else if(query1[i].ch=='Q') n1++,ans+=(double)query(root[find(x)],y);
        }
        if(n1) printf("Case %d: %.6lf\n",T,ans/n1);
        else printf("Case %d: 0\n",T); //这句话似乎...不重要【捂脸】
    }
    return 0;
}