bzoj3514 Codechef MARCH14 GERALD07加强版 LCT&&主席树

wulala

葱娘说这是一个很巧妙的题。。
有一个比较猎奇的做法：首先把边依次加到图中，若当前这条边与图中的边形成了环，那么把这个环中最早加进来的边弹出去
并将每条边把哪条边弹了出去记录下来：ntr[i] = j，特别地，要是没有弹出边，ntr[i] = 0;
这个显然是可以用LCT来弄的对吧。
然后对于每个询问，我们的答案就是对l~r中ntr小于l的边求和，并用n减去这个值
正确性可以YY一下：
如果一条边的ntr >= l,那么显然他可以与从l ~ r中的边形成环，那么它对答案没有贡献
反之如果一条边的ntr < l那么它与从l ~ r中的边是不能形成环的，那么他对答案的贡献为-1
对于查询从l ~ r中有多少边的ntr小于l，我反正是用的函数式线段树

       我觉得这个已经讲的很清楚了。

       然后发现自己最近做题自带2倍常数。。（被STL卡了以后就一直这样）。

AC代码如下（以前LCT的模板就是学(chao)黄学长的。。主席树也是学(chao)黄学长的。。然后就看起来和黄学长一模一样了o(╯□╰)o。。。23333）：

#include<iostream>
#include<cstdio>
#include<cstring>
#define inf 1000000000
#define N 400005
using namespace std;

int n,m,trtot,cnt,tp,ex[N],ey[N],stk[N],fa[N],c[N][2],pre[N],p[N],val[N];
int ls[N*10],rs[N*10],sum[N*10],rt[N];
bool rev[N];
int read(){
	int x=0; char ch=getchar();
	while (ch<'0' || ch>'9') ch=getchar();
	while (ch>='0' && ch<='9'){ x=x*10+ch-'0'; ch=getchar(); }
	return x;
}
bool isrt(int x){
	return c[fa[x]][0]^x && c[fa[x]][1]^x;
}
int les(int x,int y){
	if (val[p[x]]<val[p[y]]) return p[x];
		else return p[y];
}
void maintain(int k){
	p[k]=k;
	p[k]=les(k,c[k][0]); p[k]=les(k,c[k][1]);
}
void rotate(int x){
	int y=fa[x],z=fa[y],l=(c[y][0]==x)?0:1,r=l^1;
	if (!isrt(y)){
		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;
	maintain(y); maintain(x);
}
void pushdown(int x){
	if (rev[x]){
		swap(c[x][0],c[x][1]);
		rev[x]=0; rev[c[x][0]]^=1; rev[c[x][1]]^=1;
	}
}
void splay(int x){
	stk[tp=1]=x; int i,y,z;
	for (i=x; !isrt(i); i=fa[i]) stk[++tp]=fa[i];
	while (tp) pushdown(stk[tp--]);
	while (!isrt(x)){
		y=fa[x]; z=fa[y];
		if (!isrt(y)){
			if (c[y][0]==x ^ c[z][0]==y) rotate(x);
				else rotate(y);
		}
		rotate(x);
	}
}
void access(int x){
	int y=0;
	for (; x; y=x,x=fa[x]){
		splay(x); c[x][1]=y; maintain(x);
	}
}
void makert(int x){
	access(x); splay(x); rev[x]^=1;
}
void link(int x,int y){
	makert(x); fa[x]=y;
}
void cut(int x,int y){
	makert(x); access(y); splay(y);
	c[y][0]=fa[x]=0;
}
int getanc(int x){
	access(x); splay(x);
	while (c[x][0]) x=c[x][0]; return x;
}
int getmin(int x,int y){
	makert(x); access(y); splay(y);
	return p[y];
}
void ins(int l,int r,int x,int &y,int v){
	y=++trtot; sum[y]=sum[x]+1;
	if (l==r) return; int mid=(l+r)>>1;
	if (v<=mid){ rs[y]=rs[x]; ins(l,mid,ls[x],ls[y],v); }
	else{ ls[y]=ls[x]; ins(mid+1,r,rs[x],rs[y],v); }
}
int qry(int x,int y,int l,int r,int v){
	if (r==v) return sum[y]-sum[x]; int mid=(l+r)>>1;
	if (v<=mid) return qry(ls[x],ls[y],l,mid,v);
	else return sum[ls[y]]-sum[ls[x]]+qry(rs[x],rs[y],mid+1,r,v);
}
int main(){
	n=read(); m=read(); int cas=read(),k=read(),i;
	for (i=0; i<=n; i++) val[p[i]=i]=inf; cnt=n;
	for (i=1; i<=m; i++){
		int x=read(),y=read(); ex[i]=x; ey[i]=y;
		if (x==y){ pre[i]=i; continue; }
		if (getanc(x)==getanc(y)){
			int z=getmin(x,y); pre[i]=val[z];
			cut(ex[pre[i]],z); cut(ey[pre[i]],z);
		}
		val[++cnt]=i; p[cnt]=cnt;
		link(x,cnt); link(y,cnt);
	}
	for (i=1; i<=m; i++) ins(0,m,rt[i-1],rt[i],pre[i]);
	int ans=0;
	while (cas--){
		int u=read(),v=read();
		if (k){ u^=ans; v^=ans; }
		printf("%d\n",ans=n-qry(rt[u-1],rt[v],0,m,u-1));
	}
	return 0;
}

by lych

2016.2.27