【BZOJ3514】Codechef MARCH14 GERALD07加强版

题目

这是一道好题。（只是我这种蒟蒻这样认为，求大牛勿喷）

先不考虑L和R的限制。一开始每个点都是一个联通块，则答案是n，我们发现如果加进来一条边，边的两边之前没有联通，那么答案就会减一（及贡献为-1）。如果边的两边已经联通了，则对答案没有贡献。

现在有了L和R的限制。相当于原来有一些对答案没有贡献的边可能会对答案有-1的贡献，那什么样的边会对答案有影响呢？我们对每条边求一个值po[i],即在加第i条边时，如果若L等于po[i]+1,则此边会对答案产生影响，怎么求呢？？若两端没有联通，则po[i]=0；若联通了，我就求最早的那个罪魁祸首的编号K（即若L=K+1的时候再加这条边，这条边就会对答案产生-1的贡献），则po[i]=K，并且要删除原来那条边（因为现在又是一个新的联通块辣~~）。

有了po数组我们就知道了对于每条边，L等于好多的时候此边会对答案产生影响。很明显po要用LCT求，化边为点，令此点权为这条边的标号，原来图中的点点权为+inf，维护一个minn，即是那个最小的罪魁祸首。

现在给出L和R，则答案为n-sum(sum为在L~R内且po值小于L的个数)（明显的主席树啊233333）

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define N 400005
#define inf 0x7fffffff/2
using namespace std;
struct eee{int x,y;}edge
;
int n,m,op,Q,mi
,val
,fz
,fa
,cnt,po
,c
[2],au
;

int rt
,sum[N*19],num=0,ls[N*19],rs[N*19];
void insert(int la,int &rt,int l,int r,int v){
int mid=(l+r)/2;
rt=++num;
sum[rt]=sum[la];sum[rt]++;
ls[rt]=ls[la];rs[rt]=rs[la];
if(l==r)return;
if(v<=mid)insert(ls[la],ls[rt],l,mid,v);
else insert(rs[la],rs[rt],mid+1,r,v);
}
int ask(int r1,int r2,int l,int r,int lim){
int mid=(l+r)/2;
if(l>=lim)return 0;
if(r<lim)return sum[r2]-sum[r1];
return ask(ls[r1],ls[r2],l,mid,lim)+ask(rs[r1],rs[r2],mid+1,r,lim);
}
void pushup(int x){
int l=c[x][0],r=c[x][1];
mi[x]=val[x];
if(l)mi[x]=min(mi[x],mi[l]);
if(r)mi[x]=min(mi[x],mi[r]);
}
void pushdown(int x){
int l=c[x][0],r=c[x][1];
if(fz[x]){
fz[x]=0;
fz[l]^=1;fz[r]^=1;
swap(c[x][0],c[x][1]);
}
}
void rot(int x){
int y=fa[x],z=fa[y];
int l=(x!=c[y][0]),r=l^1;
if(!au[y]){au[x]=0;au[y]=1;}
else {if(c[z][0]==y)c[z][0]=x;else c[z][1]=x;}
fa[x]=z;fa[y]=x;fa[c[x][r]]=y;
c[y][l]=c[x][r];c[x][r]=y;
pushup(y);pushup(x);
}
void splay(int x){
pushdown(x);
while(au[x]){
int y=fa[x],z=fa[y];
if(au[y])pushdown(z);pushdown(y);pushdown(x);
if(au[y]){if(c[z][0]==y^c[y][0]==x)rot(x);else rot(y);}
rot(x);
}
}
void access(int x){
int y=0;
while(x){
splay(x);
au[c[x][1]]=0;au[y]=1;
fa[y]=x;c[x][1]=y;
pushup(x);
y=x;x=fa[x];
}
}
void evert(int x){access(x);splay(x);fz[x]^=1;}
void link(int x,int y){evert(y);fa[y]=x;}
void cut(int x,int y){
evert(x);access(y);splay(y);
au[c[y][0]]=fa[c[y][0]]=0;c[y][0]=0;
}
int find_rt(int x){
access(x);splay(x);
pushdown(x);
while(c[x][0]){x=c[x][0];pushdown(x);}
return x;
}
void Pre_lct(){
int i,j,k;
for(i=1;i<=m;i++){
cnt++;//华点：之前放在else里面去了，T了一下午23333
if(edge[i].x==edge[i].y)po[i]=i;
else {
mi[cnt]=val[cnt]=cnt;
if(find_rt(edge[i].x)==find_rt(edge[i].y)){
evert(edge[i].x);access(edge[i].y);
splay(edge[i].y);
int u1=mi[edge[i].y],u2=u1-n;
cut(u1,edge[u2].x);
cut(u1,edge[u2].y);
po[i]=u2;
}
link(cnt,edge[i].x);
link(cnt,edge[i].y);
}
}
for(i=1;i<=m;i++)insert(rt[i-1],rt[i],0,m,po[i]);
}
void Que(){
int i,j,k,x,y,ans=0;
for(i=1;i<=Q;i++){
scanf("%d%d",&x,&y);
if(op){x^=ans;y^=ans;}
ans=n-ask(rt[x-1],rt[y],0,m,x);
printf("%d\n",ans);
}
}
int main(){
int i,j,k,x,y;
scanf("%d%d%d%d",&n,&m,&Q,&op);
for(i=1;i<=n;i++){mi[i]=val[i]=inf;}//华点：之前写成了m
mi[0]=val[0]=inf;
cnt=n;//华点：之前写成m了
for(i=1;i<=m;i++)
scanf("%d%d",&edge[i].x,&edge[i].y);
Pre_lct();
Que();
return 0;
}

注解： 华点：出自福尔摩斯：“盲生，我发现了华点。”