bzoj4025 二分图（lct）

如何维护这个图是不是二分图呢，我们就来维护一下图中奇环的数量。我们把每条边按扫描线那样分成两个事件来处理。把所有事件按时间排序。我们维护一棵树，每条边的权值是这条边的删除时间，我们维护一棵类似最大生成树的树。这样加进一条边以后只需要与当前树上的边进行判断即可。遇到加边事件：
1)如果加的是一条树边，就直接加树边。
2)如果是一条非树边，替换掉环中最小的一条边，然后:
1.构成奇环，奇环个数cnt+1，最小的那条边标记1，删除这条边时奇环个数-1
2.构成偶环，嘛也不干。
遇到删边事件：
1)删掉一条树边，根据我们维护的最大生成树，一定没有非树边能够连接起这两个连通块了，所以我们只能断开它们。
2)删掉一条非树边，看标记，如果带着奇环，就奇环个数-1
如果奇环个数为0则是二分图，否则不是。
tips：对于边权的处理：类似魔法森林，化边权为点权。
要维护路径最小边和路径长度。

#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 300010
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,T,fa[N],c[N][2],v[N],mn[N],num=0,cnt=0,tid[N],q[N],sz[N];
bool rev[N];int mark[N];//0--树边，1--奇环的非树边，2--偶环的非树边
struct Event{
    int x,y,t,op,id;
    friend bool operator<(Event a,Event b){return a.t<b.t;}
}E[400010];
inline bool isroot(int x){
    return x!=c[fa[x]][0]&&x!=c[fa[x]][1];
}
inline void update(int x){
    int l=c[x][0],r=c[x][1];
    sz[x]=sz[l]+sz[r]+1;
    if(v[x]<=v[mn[l]]&&v[x]<=v[mn[r]]) mn[x]=x;
    else if(v[mn[l]]<=v[mn[r]]) mn[x]=mn[l];
    else mn[x]=mn[r];
}
inline void pushdown(int p){
    if(!rev[p]) return;rev[p]=0;int l=c[p][0],r=c[p][1];
    rev[l]^=1;rev[r]^=1;swap(c[p][0],c[p][1]);
}
inline void rotate(int x){
    int y=fa[x],z=fa[y],l=x==c[y][1],r=l^1;
    if(!isroot(y)) c[z][y==c[z][1]]=x;
    fa[c[x][r]]=y;fa[y]=x;fa[x]=z;
    c[y][l]=c[x][r];c[x][r]=y;update(y);update(x);
}
inline void splay(int x){
    int top=0;q[++top]=x;
    for(int xx=x;!isroot(xx);xx=fa[xx]) q[++top]=fa[xx];
    while(top) pushdown(q[top--]);
    while(!isroot(x)){
        int y=fa[x],z=fa[y];
        if(!isroot(y)){
            if(x==c[y][1]^y==c[z][1]) rotate(x);
            else rotate(y);
        }rotate(x);
    }
}
inline void access(int x){
    int y=0;
    while(x){splay(x);c[x][1]=y;update(x);y=x;x=fa[x];}
}
inline void makeroot(int x){
    access(x);splay(x);rev[x]^=1;
}
inline void link(int x,int y){
    makeroot(x);fa[x]=y;
}
inline void cut(int x,int y){
    makeroot(x);access(y);splay(y);c[y][0]=fa[x]=0;update(y);
}
inline int find(int x){
    access(x);splay(x);while(c[x][0]) x=c[x][0];return x;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();T=read();
    for(int i=0;i<=n;++i) v[i]=inf;
    for(int i=1;i<=m;++i){
        int x=read(),y=read(),t1=read(),t2=read();v[i+n]=t2;
        E[++num].t=t1;E[num].op=1;E[num].x=x;E[num].y=y;E[num].id=i;
        E[++num].t=t2;E[num].op=0;E[num].x=x;E[num].y=y;E[num].id=i;
    }sort(E+1,E+num+1);int now=1;
    for(int tim=1;tim<=T;++tim){
        while(now<=num&&E[now].t<tim){
            tid[E[now].id]=now;
            if(E[now].op){
                int x=E[now].x,y=E[now].y,c=E[now].id+n;
                int xx=find(x),yy=find(y);
                if(xx!=yy){link(x,c),link(c,y);++now;continue;}
                makeroot(x);access(y);splay(y);int z=mn[y];
                if((sz[y]>>1)+1&1){
                    if(v[z]<v[c]){cut(E[tid[z-n]].x,z);cut(E[tid[z-n]].y,z);link(x,c);link(c,y);mark[z-n]=1;}
                    else mark[c-n]=1;++cnt;
                }else{
                    if(v[z]<v[c]){cut(E[tid[z-n]].x,z);cut(E[tid[z-n]].y,z);link(x,c);link(c,y);mark[z-n]=2;}
                    else mark[c-n]=2;
                }
            }else{
                int c=E[now].id;
                if(!mark[c]) cut(E[now].x,c+n),cut(c+n,E[now].y);
                else if(mark[c]==1) --cnt;
            }++now;
        }puts(cnt?"No":"Yes");
    }return 0;
}