bzoj4399: 魔法少女LJJ
Description
在森林中见过会动的树，在沙漠中见过会动的仙人掌过后，魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了
LJJ感叹道“这里真是个迷人的绿色世界,空气清新、淡雅,到处散发着醉人的奶浆味；小猴在枝头悠来荡去,好不自在；各式各样的鲜花争相开放,各种树枝的枝头挂满沉甸甸的野果；鸟儿的歌声婉转动听,小河里飘着落下的花瓣真是人间仙境”
SHY觉得LJJ还是太naive，一天，SHY带着自己心爱的图找到LJJ，对LJJ说：“既然你已经见识过动态树，动态仙人掌了，那么今天就来见识一下动态图吧”
LJJ：“要支持什么操作？”
SHY：“
1.新建一个节点，权值为x。
2.连接两个节点。
3.将一个节点a所属于的联通快内权值小于x的所有节点权值变成x。
4.将一个节点a所属于的联通快内权值大于x的所有节点权值变成x。
5.询问一个节点a所属于的联通块内的第k小的权值是多少。
6.询问一个节点a所属联通快内所有节点权值之积与另一个节点b所属联通快内所有节点权值之积的大小。
7.询问a所在联通快内节点的数量
8.若两个节点a，b直接相连，将这条边断开。
9.若节点a存在，将这个点删去。
”
LJJ：“我可以离线吗？”
SHY：“可以，每次操作是不加密的，”
LJJ：“我可以暴力吗？”
SHY：“自重”
LJJ很郁闷，你能帮帮他吗
Input
第一行有一个正整数m，表示操作个数。
接下来m行，每行先给出1个正整数c。
若c=1，之后一个正整数x，表示新建一个权值为x的节点，并且节点编号为n+1（当前有n个节点）。
若c=2，之后两个正整数a，b，表示在a，b之间连接一条边。
若c=3，之后两个正整数a，x，表示a联通快内原本权值小于x的节点全部变成x。
若c=4，之后两个正整数a，x，表示a联通快内原本权值大于x的节点全部变成x。
若c=5，之后两个正整数a，k，表示询问a所属于的联通块内的第k小的权值是多少。
若c=6，之后两个正整数a，b，表示询问a所属联通快内所有节点权值之积与b所属联通快内所有节点权值之积的大小，
若a所属联通快内所有节点权值之积大于b所属联通快内所有节点权值之积，输出1，否则为0。
若c=7，之后一个正整数a，表示询问a所在联通块大小
若c=8，之后两个正整数a，b，表示断开a，b所连接的边。
若c=9，之后一个正整数a，表示断开a点的所有连边
具体输出格式见样例
Output
Sample Input
1 2
1 3
1 4
1 5
1 6
2 1 2
2 2 3
2 3 4
2 4 5
9 1
3 2 5
5 3 4
Sample Output
HINT
对100%的数据 0<=m<=400000,c<=7,所有出现的数均<=1000000000,所有出现的点保证存在
【HINT】请认真阅读题面
Source
啊4399 小游戏？
【HINT】请认真阅读题面 c<=7
于是我们只需支持合并。
考虑用权值线段树维护<=k的有几个，合并时就线段树合并。
修改是先查询size，再在x点打上标记。
这时线段树合并先pushdown()再合并。
还剩一个棘手的操作6。权值太大存不下，我们取ln(x),然后比较ln的大小。
注意注意注意：线段树合并不能用pushup(x)来维护，边界会炸。
调的欲哭无泪啊
#include<cstdio> #include<iostream> #include<cstdlib> #include<cstring> #include<algorithm> #include<cmath> #define maxn 400005 #define Max 1000000000 using namespace std; int m,op,t1,t2,t3,f[maxn],tot,cnt; int rt[maxn]; struct node{ int ls,rs,sz,bj;double v; }tr[maxn*10]; int getf(int k){return f[k]==k?k:f[k]=getf(f[k]);} void wh(int k){ int ls=tr[k].ls,rs=tr[k].rs; tr[k].sz=tr[ls].sz+tr[rs].sz; tr[k].v=tr[ls].v+tr[rs].v; } void upd(int k){ if(!k)return; tr[k].sz=0;tr[k].v=0; tr[k].bj=1; } void down(int k){ if(tr[k].bj){ upd(tr[k].ls);upd(tr[k].rs); tr[k].bj=0; } } void add(int &k,int l,int r,int pl,int v){ if(!k)k=++cnt; if(l==r){ tr[k].sz=v;tr[k].v=log(pl)*v; return; } down(k); int mid=l+r>>1; if(pl<=mid)add(tr[k].ls,l,mid,pl,v); else add(tr[k].rs,mid+1,r,pl,v); wh(k); } int merge(int x,int y){ if(!x||!y)return x+y; down(x);down(y); tr[x].sz+=tr[y].sz;tr[x].v+=tr[y].v; tr[x].ls=merge(tr[x].ls,tr[y].ls); tr[x].rs=merge(tr[x].rs,tr[y].rs); return x; } int lian(int k,int l,int r,int li,int ri){ if(!k)return 0; if(l>=li&&r<=ri){ int v=tr[k].sz; upd(k);return v; } down(k); int mid=l+r>>1;int sz=0; if(li<=mid)sz+=lian(tr[k].ls,l,mid,li,ri); if(ri>mid)sz+=lian(tr[k].rs,mid+1,r,li,ri); wh(k); return sz; } int Kth(int k,int l,int r,int kth){ if(l==r)return l; down(k); int mid=l+r>>1; int sl=tr[tr[k].ls].sz; if(sl>=kth)return Kth(tr[k].ls,l,mid,kth); else return Kth(tr[k].rs,mid+1,r,kth-sl); } int main() { cin>>m; for(int i=1;i<=m;i++)f[i]=i; for(int i=1;i<=m;i++){ scanf("%d",&op); if(op==1){ scanf("%d",&t3); tot++; add(rt[tot],1,Max,t3,1); } if(op==2){ scanf("%d%d",&t1,&t2); int f1=getf(t1),f2=getf(t2); if(f1!=f2){ rt[f1]=merge(rt[f1],rt[f2]); f[f2]=f1; } } if(op==3){ scanf("%d%d",&t1,&t2); t1=getf(t1); int siz=lian(rt[t1],1,Max,1,t2); add(rt[t1],1,Max,t2,siz); } if(op==4){ scanf("%d%d",&t1,&t2); t1=getf(t1); int siz=lian(rt[t1],1,Max,t2,Max); add(rt[t1],1,Max,t2,siz); } if(op==5){ scanf("%d%d",&t1,&t2); t1=getf(t1); printf("%d ",Kth(rt[t1],1,Max,t2)); } if(op==6){ scanf("%d%d",&t1,&t2); t1=getf(t1);t2=getf(t2); if(tr[rt[t1]].v>tr[rt[t2]].v)puts("1"); else puts("0"); } if(op==7){ scanf("%d",&t1); t1=getf(t1); printf("%d ",tr[rt[t1]].sz); } } return 0; }