BZOJ 3514 Codechef MARCH14 GERALD07加强版 LCT+主席树

题意:链接

方法: LCT

解析:

被这个名字长到恶心的东西恶心了一个晚上。

不过我终于A掉他辣哈哈哈哈哈哈哈哈哈哈哈

第一遍读完题我就不知道说啥好了，这tm什么玩意。

这咋求啊？

然后还顺便百度了下连通块=-=

题太神不知道咋做啊。

然后有个特别奇♂怪的ntr算法。

就是每一次加边的时候，如果自环那就ntr自己，对答案不会有影响。肯定嘛这是，相当于没连边。

非自环呢？

首先判断这条边加上去之后对答案有没有什么影响。

如果加上去不产生环的话，那显然是使连通块-1。

如果产生环的话，那显然是对答案没有影响。

如果产生环，那这条边的ntr值就是环上最早出现的边的值。

代表它弹掉了这条边，当然我们需要弹掉这条边。

这证明了什么！

在这条被弹掉的边之后出现，并且在弹掉他的边之前出现的边显然是不能与弹掉他的边构成环的！

也就是说如果区间左端点大于这个被弹掉的边的话，弹掉他的边的出现一定会使答案减1。

如果区间里包含被弹掉的边的与弹掉他的边话，那么弹掉他的边显然是对答案没有影响的！他的出现是不会使连通块-1的。

脑补不出来就画图就好辣！并没有辣么难！

然后问题转化成什么了？

先用LCT预处理出所有ntr值。

然后在ntr数组上[l,r]的区间询问ntr值小于l的数的个数！

主席树裸上！

然后我就T了。T_T

总结下T的经验吧。

我T在两个地方。

刚开始空间开的是够的，然后find_root太暴力被卡。

后来莫名其妙缩了下空间？find_root改写了下又被卡。

然后我就一直在研究奇怪的东西。

后来发现对拍为甚么太大会错？

…..主席树的空间啊混蛋！nm让我莫名其妙缩成了nm/2。

好题！

代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 4001000
#define INF 0x3f3f3f3f
using namespace std;
int n,m,k,type,size;
int rev[N],ch[N][2],fa[N],root[N];
int mi[N],val[N],rt[N],sum[N],lson[N],rson[N];
int ntr[N];
struct node
{
    int x,y;
}edge[N];
inline int read()
{
    int x=0;
    char ch=getchar();
    while(ch<'0'||ch>'9'){ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return x;
}
inline void pushup(int x)
{
    if(!x)return;
    mi[x]=x;
    if(ch[x][0]!=0&&val[mi[ch[x][0]]]<val[mi[x]])mi[x]=mi[ch[x][0]];
    if(ch[x][1]!=0&&val[mi[ch[x][1]]]<val[mi[x]])mi[x]=mi[ch[x][1]];
}
inline void reverse(int x)
{
    if(!x)return;
    swap(ch[x][0],ch[x][1]);
    rev[x]^=1;
}
inline void pushdown(int x)
{
    if(!x)return;
    if(rev[x])
    {
        reverse(ch[x][0]);
        reverse(ch[x][1]);
        rev[x]^=1;
    }
}
inline void down(int x)
{
    if(!x)return;
    if(!rt[x])down(fa[x]);
    pushdown(x);
}
inline void rotate(int x)
{
    int y=fa[x],kind=ch[y][1]==x;
    ch[y][kind]=ch[x][!kind];
    fa[ch[y][kind]]=y;
    ch[x][!kind]=y;
    fa[x]=fa[y];
    fa[y]=x;
    if(rt[y])rt[y]=0,rt[x]=1;
    else ch[fa[x]][ch[fa[x]][1]==y]=x;
    pushup(y);
}
inline void splay(int x)
{
    down(x);
    while(!rt[x])
    {
        int y=fa[x],z=fa[y];
        if(rt[y])rotate(x);
        else if((ch[y][1]==x)==(ch[z][1]==y))rotate(y),rotate(x);
        else rotate(x),rotate(x);
    }
    pushup(x);
}
inline void access(int x)
{
    int y=0;
    while(x)
    {
        splay(x);
        rt[ch[x][1]]=1,rt[y]=0;
        ch[x][1]=y;
        pushup(x);
        y=x,x=fa[x];
    }
}
inline void movetoroot(int x)
{
    access(x);
    splay(x);
    reverse(x);
}
inline int find_root(int x)
{
    access(x),splay(x);
    while(ch[x][0])x=ch[x][0];
    return x;
}
inline void link(int x,int y)
{
    movetoroot(x);
    fa[x]=y;
}
inline void cut(int x,int y)
{
    movetoroot(x);
    access(y);
    splay(y);
    ch[y][0]=0;
    fa[x]=0;
    rt[x]=1;
    pushup(y); 
}
inline void init()
{
    for(int i=1;i<=n+m;i++)rt[i]=1,val[i]=INF;
}
inline void insert(int l,int r,int x,int &y,int val)
{
    y=++size;
    sum[y]=sum[x]+1;
    if(l==r)return;
    lson[y]=lson[x],rson[y]=rson[x];
    int mid=(l+r)>>1;
    if(val<=mid)insert(l,mid,lson[x],lson[y],val);
    else insert(mid+1,r,rson[x],rson[y],val);
}
inline int query(int l,int r,int x,int y,int val)
{
    if(r==val)return sum[y]-sum[x];
    int mid=(l+r)>>1;
    if(val<=mid)return query(l,mid,lson[x],lson[y],val);
    else return sum[lson[y]]-sum[lson[x]]+query(mid+1,r,rson[x],rson[y],val);
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&k,&type);
    init();
    for(int i=1;i<=m;i++)
    {
        int x,y;
        edge[i].x=read();edge[i].y=read();
        x=edge[i].x,y=edge[i].y;
        if(x==y){ntr[i]=i;}
        else if(find_root(x)==find_root(y))
        {
            movetoroot(x),access(y),splay(y);
            int tmp=mi[y];
            ntr[i]=tmp-n;
            cut(edge[tmp-n].x,tmp),cut(tmp,edge[tmp-n].y);
            val[n+i]=i,mi[n+i]=n+i;
            link(edge[i].x,i+n),link(i+n,edge[i].y);
        }else
        {
            val[n+i]=i,mi[n+i]=n+i;
            link(edge[i].x,i+n),link(i+n,edge[i].y);
        }
    }
    for(int i=1;i<=m;i++)
    {
        insert(0,m,root[i-1],root[i],ntr[i]);
    }
    int ans=0;
    for(int i=1;i<=k;i++)
    {
        int l,r;
        l=read();r=read();
        if(type)l^=ans,r^=ans;
        ans=n-query(0,m,root[l-1],root[r],l-1);
        printf("%d\n",ans);
    }
}