【bzoj3514】GERALD07加强版 LCT+主席树

AC通道：http://www.lydsy.com/JudgeOnline/problem.php?id=3514

【题解】

来自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，我反正是用的函数式线段树

我就想不到啊，QAQ

果然是我弱。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
using namespace std;
#define FILE "read"
#define MAXN 400010
#define INF 1000000000
#define up(i,j,n) for(int i=j;i<=n;++i)
#define dn(i,j,n) for(int i=j;i>=n;--i)
namespace INIT{
	char buf[1<<15],*fs,*ft;
	inline char getc(){return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;}
	inline int read(){
		int x=0,f=1;  char ch=getc();
		while(!isdigit(ch))  {if(ch=='-')  f=-1;  ch=getc();}
		while(isdigit(ch))  {x=x*10+ch-'0';  ch=getc();}
		return x*f;
	}
}using namespace INIT;
int n,m,q,ans,tot,type,v[MAXN],st[MAXN];
namespace Link_Cut_Tree{
	int stack[MAXN],f[MAXN],vis[MAXN],minn[MAXN],son[MAXN][2];
	bool get(int x){return son[f[x]][1]==x;}
	bool isroot(int x){return son[f[x]][0]!=x&&son[f[x]][1]!=x;}
	void pushdown(int x){
		if(!vis[x])  return;
		swap(son[x][0],son[x][1]);
		vis[son[x][0]]^=1; vis[son[x][1]]^=1; vis[x]=0;
	}
	void updata(int x){
		minn[x]=x;
		minn[x]=v[minn[son[x][0]]]<v[minn[x]]?minn[son[x][0]]:minn[x];
		minn[x]=v[minn[son[x][1]]]<v[minn[x]]?minn[son[x][1]]:minn[x];
	}
	void rotate(int x){
		int y=f[x],z=f[y],which=get(x);
		if(!isroot(y))  son[z][son[z][1]==y]=x;
		son[y][which]=son[x][which^1]; f[son[y][which]]=y;
		son[x][which^1]=y;  f[y]=x; f[x]=z;
		updata(y);  updata(x);
	}
	void splay(int x){
		int top(0); stack[++top]=x;
		for(int i=x;!isroot(i);i=f[i]) stack[++top]=f[i];
		dn(i,top,1)  pushdown(stack[i]);
		for(int y=f[x];!isroot(x);rotate(x),y=f[x])
			if(!isroot(y))  rotate(get(x)==get(y)?y:x);
	}
	void access(int x){for(int temp(0);x;temp=x,x=f[x])splay(x),son[x][1]=temp,updata(x);}
	void reverse(int x){access(x);splay(x);vis[x]^=1;}
	void linkk(int x,int y){reverse(x);f[x]=y;}
	void cut(int x,int y){reverse(x);access(y);splay(y);f[x]=son[y][0]=0;}
	int find(int x){access(x);splay(x);while(son[x][0])x=son[x][0];return x;}
	int query(int x,int y){reverse(x);access(y);splay(y);return minn[y];}
}
namespace Chairman_Tree{
	int cnt,root[MAXN],sum[MAXN<<5],son[MAXN<<5][2];
	int newroot(int last){
		sum[++cnt]=sum[last]+1;
		son[cnt][0]=son[last][0]; 
		son[cnt][1]=son[last][1];
		return cnt;
	}
	void insert(int l,int r,int &root,int last,int x){
		root=newroot(last);
		if(l==r)  return;
		int mid=(l+r)>>1;
		if(x<=mid)  insert(l,mid,son[root][0],son[last][0],x);
		else insert(mid+1,r,son[root][1],son[last][1],x);
	}
	int ask(int l,int r,int x,int y,int v){
		if(r==v)  return sum[y]-sum[x];
		int mid=(l+r)>>1,temp=sum[son[y][0]]-sum[son[x][0]];
		if(mid>=v)  return ask(l,mid,son[x][0],son[y][0],v);
		else return temp+ask(mid+1,r,son[x][1],son[y][1],v);
	}
}
int main(){
	freopen(FILE".in","r",stdin);
	freopen(FILE".out","w",stdout);
	using namespace Link_Cut_Tree;
	n=read();  m=read();  q=read();  type=read();
	up(i,0,n)  minn[i]=i,v[i]=INF;  tot=n;
	up(i,1,m){
		int x=read(),y=read();
		if(x==y)  {st[i]=i; continue;}
		if(find(x)==find(y)){
			int t=query(x,y);  st[i]=v[t];
			cut(x,t);  cut(y,t);
		}
		minn[++tot]=tot;  v[tot]=i;
		linkk(x,tot);  linkk(y,tot);
	}
	using namespace Chairman_Tree;
	up(i,1,m)  insert(0,m,root[i],root[i-1],st[i]);
	up(i,1,q){
		int x=read(),y=read();
		if(type)  x^=ans,y^=ans;
		ans=n-ask(0,m,root[x-1],root[y],x-1);
		printf("%d\n",ans);
	}
	return 0;
}