[BZOJ4399]魔法少女LJJ----------线段树进阶
感谢线段树进阶,给了我重新做人的机会。---------------某不知名OIer,Keen_z
Description
题目描述
输入格式
HINT
对100%的数据 0<=m<=400000,c<=7,所有出现的数均<=1000000000,所有出现的点保证存在 【HINT】请认真阅读题面
Solution
给了8,9两个阴间删点删边操作之后告诉你c<=7。。
一个很喵的转化,可以将连乘转化为对数之和。
数据最大1e9,显然不能直接开线段树,需要离散化。
将每个联通块都看作一个线段树,在其中维护权值为i的点的个数与元素对数的和。连边时将两点所在的线段树合并。
过程中要用一波并察集操作,来及时更新合并后线段树的根。
主要还是看代码吧,应该挺好理解的(挺不好调的
附去掉9 1后的样例输入,输出应为5。
code:
#include<bits/stdc++.h> #define debug cout<<"wrong"<<endl using namespace std; const int NN=4e5+5; int hal,id,c[NN],rt[NN],m,dat[NN][3],has[NN],cnt,fa[NN],s; double logg[NN]; inline signed 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<<1)+(x<<3)+(ch^48); ch=getchar(); } return x*f; } void write(int x){ if(x<0) putchar('-'), x=-x; if(x>9) write(x/10); putchar(x%10+'0'); } int get(int x){ return fa[x]==x?x:fa[x]=get(fa[x]); } inline int cag(int x){ return lower_bound(has+1,has+1+hal,x)-has; } struct node{ int seg,ls[NN*20],rs[NN*20],caf[NN*20],siz[NN*20]; double sum[NN*20]; void pushup(int x){ sum[x]=sum[ls[x]]+sum[rs[x]]; siz[x]=siz[ls[x]]+siz[rs[x]]; } void insert(int &x,int l,int r,int pos,int w){ if(!x) x=++seg; if(l==r){ sum[x]+=logg[pos]*w; siz[x]+=w; return; } int mid=(l+r)>>1; if(pos<=mid) insert(ls[x],l,mid,pos,w); else insert(rs[x],mid+1,r,pos,w); pushup(x); } int query(int x,int l,int r,int pos){ if(l==r) return l; int mid=(l+r)>>1; if(siz[ls[x]]>=pos) return query(ls[x],l,mid,pos); else return query(rs[x],mid+1,r,pos-siz[ls[x]]); } void marge(int &x,int y,int l,int r){ if(!x||!y){ x=x+y; return; } if(l==r){ siz[x]+=siz[y]; sum[x]+=sum[y]; return; } int mid=(l+r)>>1; marge(ls[x],ls[y],l,mid); marge(rs[x],rs[y],mid+1,r); pushup(x); } void delet(int x,int l,int r,int opl,int opr){ if(!siz[x]) return; if(l==r){ s+=siz[x]; sum[x]=siz[x]=0; return; } int mid=(l+r)>>1; if(opl<=mid) delet(ls[x],l,mid,opl,opr); if(opr>mid) delet(rs[x],mid+1,r,opl,opr); pushup(x); } }segt; int main(){ m=read(); for(int i=1;i<=m;i++){ c[i]=read(); dat[i][1]=read(); if(c[i]!=1&&c[i]!=7) dat[i][2]=read(); if(c[i]==1) has[++cnt]=dat[i][1], dat[i][2]=++id; if(c[i]==3||c[i]==4) has[++cnt]=dat[i][2]; } sort(has+1,has+cnt+1); hal=unique(has+1,has+cnt+1)-(has+1); for(int i=1;i<=m;i++) switch(c[i]){ case 1:{ int x=cag(dat[i][1]); logg[x]=log(dat[i][1]); segt.insert(rt[dat[i][2]],1,hal,x,1); fa[dat[i][2]]=dat[i][2]; break; } case 2:{ int rx=get(dat[i][1]),ry=get(dat[i][2]); if(rx!=ry){ fa[ry]=rx; segt.marge(rt[rx],rt[ry],1,hal); } break; } case 3:{ int x=cag(dat[i][2]),ro=get(dat[i][1]); s=0; logg[x]=log(dat[i][2]); segt.delet(rt[ro],1,hal,1,max(1,x-1)); if(s) segt.insert(rt[ro],1,hal,x,s); break; } case 4:{ int x=cag(dat[i][2]),ro=get(dat[i][1]); s=0; logg[x]=log(dat[i][2]); segt.delet(rt[ro],1,hal,min(hal,x+1),hal); if(s) segt.insert(rt[ro],1,hal,x,s); break; } case 5:{ int ro=get(dat[i][1]); int ans=segt.query(rt[ro],1,hal,dat[i][2]); write(has[ans]); putchar('\n'); break; } case 6:{ int r1=get(dat[i][1]),r2=get(dat[i][2]); if(segt.sum[rt[r1]]>segt.sum[rt[r2]]) puts("1"); else puts("0"); break; } case 7:{ int ro=get(dat[i][1]); write(segt.siz[rt[ro]]); putchar('\n'); break; } } } /* 11 1 2 1 3 1 4 1 5 1 6 2 1 2 2 2 3 2 3 4 2 4 5 3 2 5 5 3 4 */
1 #include<bits/stdc++.h> 2 #define debug cout<<"lbwnb"<<endl 3 using namespace std; 4 const int NN=4e5+5; 5 int hal,id,c[NN],rt[NN],m,dat[NN][3],has[NN],cnt,fa[NN],s; 6 double logg[NN]; 7 inline signed read(){ 8 int x=0,f=1; 9 char ch=getchar(); 10 while(ch<'0'||ch>'9') 11 { 12 if(ch=='-') f=-1; 13 ch=getchar(); 14 } 15 while(ch<='9'&&ch>='0') 16 { 17 x=(x<<1)+(x<<3)+(ch^48); 18 ch=getchar(); 19 } 20 return x*f; 21 } 22 void write(int x){ 23 if(x<0) putchar('-'), x=-x; 24 if(x>9) write(x/10); 25 putchar(x%10+'0'); 26 } 27 int get(int x){ 28 return fa[x]==x?x:fa[x]=get(fa[x]); 29 } 30 inline int cag(int x){ 31 return lower_bound(has+1,has+1+hal,x)-has; 32 } 33 struct node{ 34 int seg,ls[NN*20],rs[NN*20],caf[NN*20],siz[NN*20]; 35 double sum[NN*20]; 36 void pushup(int x){ 37 sum[x]=sum[ls[x]]+sum[rs[x]]; 38 siz[x]=siz[ls[x]]+siz[rs[x]]; 39 } 40 void insert(int &x,int l,int r,int pos,int w){ 41 if(!x) x=++seg; 42 if(l==r){ 43 sum[x]+=logg[pos]*w; 44 siz[x]+=w; 45 return; 46 } 47 int mid=(l+r)>>1; 48 if(pos<=mid) insert(ls[x],l,mid,pos,w); 49 else insert(rs[x],mid+1,r,pos,w); 50 pushup(x); 51 } 52 int query(int x,int l,int r,int pos){ 53 if(l==r) return l; 54 int mid=(l+r)>>1; 55 if(siz[ls[x]]>=pos) return query(ls[x],l,mid,pos); 56 else return query(rs[x],mid+1,r,pos-siz[ls[x]]); 57 } 58 void marge(int &x,int y,int l,int r){ 59 if(!x||!y){ 60 x=x+y; 61 return; 62 } 63 if(l==r){ 64 siz[x]+=siz[y]; 65 sum[x]+=sum[y]; 66 return; 67 } 68 int mid=(l+r)>>1; 69 marge(ls[x],ls[y],l,mid); 70 marge(rs[x],rs[y],mid+1,r); 71 pushup(x); 72 } 73 void delet(int x,int l,int r,int opl,int opr){ 74 if(!siz[x]) return; 75 if(l==r){ 76 s+=siz[x]; 77 sum[x]=siz[x]=0; 78 return; 79 } 80 int mid=(l+r)>>1; 81 if(opl<=mid) delet(ls[x],l,mid,opl,opr); 82 if(opr>mid) delet(rs[x],mid+1,r,opl,opr); 83 pushup(x); 84 } 85 }segt; 86 int main(){ 87 m=read(); 88 for(int i=1;i<=m;i++){ 89 c[i]=read(); dat[i][1]=read(); 90 if(c[i]!=1&&c[i]!=7) dat[i][2]=read(); 91 if(c[i]==1) has[++cnt]=dat[i][1], dat[i][2]=++id; 92 if(c[i]==3||c[i]==4) has[++cnt]=dat[i][2]; 93 } 94 sort(has+1,has+cnt+1); 95 hal=unique(has+1,has+cnt+1)-(has+1); 96 for(int i=1;i<=m;i++) 97 switch(c[i]){ 98 case 1:{ 99 int x=cag(dat[i][1]); 100 logg[x]=log(dat[i][1]); 101 segt.insert(rt[dat[i][2]],1,hal,x,1); 102 fa[dat[i][2]]=dat[i][2]; 103 break; 104 } 105 case 2:{ 106 int rx=get(dat[i][1]),ry=get(dat[i][2]); 107 if(rx!=ry){ 108 fa[ry]=rx; 109 segt.marge(rt[rx],rt[ry],1,hal); 110 } 111 break; 112 } 113 case 3:{ 114 int x=cag(dat[i][2]),ro=get(dat[i][1]); 115 s=0; logg[x]=log(dat[i][2]); 116 segt.delet(rt[ro],1,hal,1,max(1,x-1)); 117 if(s) segt.insert(rt[ro],1,hal,x,s); 118 break; 119 } 120 case 4:{ 121 int x=cag(dat[i][2]),ro=get(dat[i][1]); 122 s=0; logg[x]=log(dat[i][2]); 123 segt.delet(rt[ro],1,hal,min(hal,x+1),hal); 124 if(s) segt.insert(rt[ro],1,hal,x,s); 125 break; 126 } 127 case 5:{ 128 int ro=get(dat[i][1]); 129 int ans=segt.query(rt[ro],1,hal,dat[i][2]); 130 write(has[ans]); putchar('\n'); 131 break; 132 } 133 case 6:{ 134 int r1=get(dat[i][1]),r2=get(dat[i][2]); 135 if(segt.sum[rt[r1]]>segt.sum[rt[r2]]) puts("1"); 136 else puts("0"); 137 break; 138 } 139 case 7:{ 140 int ro=get(dat[i][1]); 141 write(segt.siz[rt[ro]]); putchar('\n'); 142 break; 143 } 144 } 145 } 146 /* 147 11 148 1 2 149 1 3 150 1 4 151 1 5 152 1 6 153 2 1 2 154 2 2 3 155 2 3 4 156 2 4 5 157 3 2 5 158 5 3 4