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
