BZOJ 4025 二分图 LCT

维护关于删除时间的最大生成树
一个奇环的影响仅取决于环上删除时间最小的边
对于每个时刻
加入一条边时
若它为树边直接加
若形成环，则弹出删除时间最小的边，若同时为奇环，则将最小边加入集合，表示这条边存在时，图中总有奇环
（偶环不用加这条边，因为若之后加入的边与这条边形成了奇环，那么新边一定能与其他边构成奇环）
删去一条边时
若为树边直接删
若在集合中，在集合中删去
这样，在每一时刻，若集合大小为0则图为二分图

一些细节：删边时加一个是否为树边的判断要快一些
边的sz为1，权值为删除时间
点sz为0，无权值
时间是0~T-1

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<queue>
#define For(i,j,k) for(int i=(j);i<=(int)k;i++)
#define Forr(i,j,k) for(int i=(j);i>=(int)k;i--)
#define Set(a,b) memset(a,b,sizeof(a))
#define L(i) (T[i].s[0])
#define R(i) (T[i].s[1])
#define F(i) (T[i].fa)
#define sz(i) (T[i].sz)
#define mi(i) (T[i].mi)
#define V(i) (T[i].val)
#define Loc(i) (R(F(i))==i)
#define getchar getchar_unlocked
using namespace std;
const int N=100010,M=200010,INF=0x3f3f3f3f;
template <class T>inline void  read(T &x){
    x=0;char c=getchar();int f(0);
    while(c<'0'||c>'9')f|=(c=='-'),c=getchar();
    while(c>='0'&&c<='9')x=(x<<3)+(x<<1)+(c^48),c=getchar();
    x=f?-x:x;
}
int n,m,T;
template <class T>inline void chkmin(T &a, T b){a=a>b?b:a;}
template <class T>inline void chkmax(T &a, T b){a=a<b?b:a;}
struct time{
    int Begin[N],to[M],Next[M],e;
    inline void init(){Set(Begin,-1);}
    inline void add(int x,int y){
        to[++e]=y,Next[e]=Begin[x],Begin[x]=e;
    }
}t1,t2;
struct Edge{
    int u,v,st,ed;
    Edge(int u=0,int v=0,int st=0,int ed=0):u(u),v(v),st(st),ed(ed){}
}E[M];
struct node{
    int s[2],fa,sz,mi,val,rev;
    node(){val=INF;}
};
struct LCT{
    node T[N+M];
    int In[M],tot,On[M];
    LCT(){T[0].sz=0;}
    inline void pushup(int x){
        mi(x)=x;
        if(L(x)&&V(mi(L(x)))<V(mi(x)))mi(x)=mi(L(x));
        if(R(x)&&V(mi(R(x)))<V(mi(x)))mi(x)=mi(R(x));
        sz(x)=sz(L(x))+sz(R(x))+(x>n);
    }
    inline void pushdown(int x){
        if(T[x].rev){
            T[x].rev^=1,T[L(x)].rev^=1,T[R(x)].rev^=1;
            swap(L(x),R(x));
        }
    }
    inline bool isrt(int x){
        return R(F(x))!=x&&L(F(x))!=x;
    }
    void Pushdown(int x){
        if(!isrt(x))Pushdown(F(x));
        pushdown(x);
    }
    inline void Rotate(int x){
        int A=F(x),B=F(A),l=Loc(x),r=l^1,d=Loc(A);
        if(!isrt(A))T[B].s[d]=x;F(x)=B;
        F(A)=x,F(T[x].s[r])=A,T[A].s[l]=T[x].s[r],T[x].s[r]=A;
        pushup(A);pushup(x);
    }
    inline void splay(int x){
        Pushdown(x);
        while(!isrt(x)){
            if(!isrt(F(x)))Rotate(x);
            Rotate(x);
        }
        pushup(x);
    }
    inline void access(int x){
        for(int i=0;x;i=x,x=F(x))
            splay(x),R(x)=i,pushup(x);
    }
    inline int findrt(int x){
        access(x),splay(x);
        while(L(x))x=L(x);
        return x;
    }
    inline void reverse(int x){
        access(x),splay(x),T[x].rev^=1;
    }
    inline void split(int x,int y){
        reverse(x),access(y),splay(y);
    }
    inline void cut(int x,int y){
        split(x,y);
        F(L(y))=0,L(y)=0;pushup(y);
    }
    inline void link(int x,int y){
        reverse(x),F(x)=y;
    }
    inline void del(int x){
        if(In[x-n])In[x-n]=0,tot--;
        else if(On[x-n])cut(x,E[x-n].u),cut(x,E[x-n].v);
    }
    inline int query1(int x,int y){
        split(x,y);
        return T[y].sz;
    }
    inline int query2(int x,int y){
        split(x,y);
        return mi(y);
    }
    inline void ins(int x){
        int u=E[x-n].u,v=E[x-n].v,d;
        if(u==v)In[x-n]=1,++tot;
        else if(findrt(u)!=findrt(v))link(u,x),link(x,v),On[x-n]=1;
        else if (V(d=query2(u,v))<V(x)){
            if(query1(u,v)&1^1)In[d-n]=1,++tot;
            cut(d,E[d-n].u),cut(d,E[d-n].v),On[d-n]=0;
            link(x,u),link(x,v);On[x-n]=1;
        }else if(query1(u,v)&1^1)In[x-n]=1,++tot,On[x-n]=0;
    }
}t;
int main(){
    read(n),read(m),read(T);
    t1.init(),t2.init();
    For(i,1,m){
        int a,b,c,d;
        read(a),read(b),read(c),read(d);
        E[i]=Edge(a,b,c,d);
        t1.add(c,i);t2.add(d,i);
        t.T[i+n].val=d;t.T[i+n].sz=1;
    }
    For(i,0,T-1){
        for(int j=t1.Begin[i];j!=-1;j=t1.Next[j])
            t.ins(t1.to[j]+n);
        for(int j=t2.Begin[i];j!=-1;j=t2.Next[j])
            t.del(t2.to[j]+n);
        puts(!t.tot?"Yes":"No");
    }
    return 0;
}