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
,h1
,now
,ch
[2],fa
,size
,v
,f
,root
,data
;
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;
}