【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[N];
int n,m,op,Q,mi[N],val[N],fz[N],fa[N],cnt,po[N],c[N][2],au[N];

int rt[N],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;
}
注解:    华点:出自福尔摩斯:“盲生,我发现了华点。”

你可能感兴趣的:(题解,主席树,LCT)