[JZOJ4769]graph

题目大意

有一个图，你要支持动态加边、删边。每个操作结束之后都要判断这个图是否是二分图。

点数$n\le3\times10^5$，操作数$m\le3\times10^5$

---

题目分析

怎么判断一个图是否是二分图？看有没有奇环就行了！
首先分析得到一个性质：如果一个（连通）图存在一个奇环，那么对于这个连通图的任意一棵生成树，都能找到一个奇环，它由一条树路径再加上一条连接路径两端的路径组成。这个反证一下，假设存在一个奇环不经过任何树路径，我们将其看成由三条非树边路径连接的三个点，然后分类讨论奇偶性得出矛盾即可。证明过程中我们还可以发现一个推论：就是这个一定存在一种新的奇环不经过原本奇环的某条边（甚至是某条路径）。
有了这个性质，我们就可以考虑这样一种做法：
离线得出每一条边的删除时间。然后按加入顺序加边，使用LCT维护图关于边删除时间的最大生成树。
对于加边操作，如果加入后图中形成了一个环，就删掉环上删除时间最小的边。破环之后，我们要看看这个环是否是奇环，如果是，就把删除的边丢进一个集合里面。
对于删边操作，如果是树边就删树边，如果在集合上就删掉集合中这个元素。
判断二分图呢？直接看看集合是否为空就好了。空就是二分图。
为什么这样是对的呢？如果我们删掉的边而且没有丢进集合的边后来和别的边构成了奇环，根据我们上面的理论，这个奇环一定能够调整为另外一个不经过这条边的奇环，且构成是生成树树边加上一个非树路径，当然，这个非树路径也一定可以调整成一条非树边。
注意这里的信息都在边上，我们要对边建点做LCT。
时间复杂度$\mathrm O((n+m)\log (n+m))$。
相当于$6\times10^5$个点做LCT，于是我T成狗，还不如$\mathrm O(m\log^2m)$的分治做法快。。。
看来我写的是LC“T”。
想过掉还是写分治吧qwq

---

代码实现

#include <algorithm>
#include <iostream>
#include <cstdio>
#include <cctype>
#include <stack>

using namespace std;

inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while (!isdigit(ch)) f=ch=='-'?-1:f,ch=getchar();
    while (isdigit(ch)) x=x*10+ch-'0',ch=getchar();
    return x*f;
}

const int N=300050;
const int M=300050;
const int V=N+M;

bool inset[M],ontree[M];
int n,m,setsize,etot;
int st[M],en[M],eid[M];
int opt[M][2];

struct link_cut_tree
{
    int fa[V],id[V],size[V],mi[V],par[V];
    int son[V][2];
    stack<int> st;
    bool mark[V];

    inline bool side(int x){return son[fa[x]][1]==x;}

    inline void update(int x)
    {
        size[x]=size[son[x][0]]+size[son[x][1]]+1;
        mi[x]=en[mi[son[x][0]]]<en[mi[son[x][1]]]?mi[son[x][0]]:mi[son[x][1]];
        mi[x]=en[mi[x]]<en[id[x]]?mi[x]:id[x];
    }

    inline void R(int x){swap(son[x][0],son[x][1]),mark[x]^=1;}

    inline void clear(int x)
    {
        if (mark[x])
        {
            if (son[x][0]) R(son[x][0]);
            if (son[x][1]) R(son[x][1]);
            mark[x]=0;
        }
    }

    inline void rotate(int x)
    {
        int y=fa[x];bool s=side(x);
        if (fa[y]) son[fa[y]][side(y)]=x;
        if (son[x][s^1]) fa[son[x][s^1]]=y;
        son[y][s]=son[x][s^1],son[x][s^1]=y;
        fa[x]=fa[y],fa[y]=x;
        if (par[y]) par[x]=par[y],par[y]=0;
        update(y),update(x);
    }

    inline void pushdown(int x,int y)
    {
        for (;x!=y;st.push(x),x=fa[x]);
        for (;!st.empty();clear(st.top()),st.pop());
    }

    inline void splay(int x,int y)
    {
        for (pushdown(x,y);fa[x]!=y;rotate(x))
            if (fa[fa[x]]!=y)
                if (side(x)==side(fa[x])) rotate(fa[x]);
                else rotate(x);
    }

    inline int access(int x)
    {
        int nxt=0;
        for (;x;update(nxt=x),x=par[x])
        {
            splay(x,0);
            if (son[x][1]) par[son[x][1]]=x,fa[son[x][1]]=0;
            if (nxt) par[nxt]=0,fa[nxt]=x;
            son[x][1]=nxt;
        }
        return nxt;
    }

    inline void makeroot(int x){R(access(x));}

    inline void link(int x,int y)
    {
        makeroot(x),access(x),splay(x,0);
        par[x]=y,access(x);
    }

    inline void cut(int x,int y)
    {
        makeroot(x),access(y),splay(y,0);
        fa[x]=son[y][0]=par[y]=0,update(y);
    }

    inline void query(int x,int y,int &len,int &minid)
    {
        makeroot(x),access(y),splay(y,0);
        len=size[y],minid=mi[y];
    }

    inline bool is_connect(int x,int y)
    {
        makeroot(x),access(y),splay(y,0);
        int tmp=par[y];
        par[y]=-1,splay(x,0);
        bool ret=par[x]==-1;
        return (ret?par[x]=tmp:par[y]=tmp),ret;
    }
}lct;

void process()
{
    for (int i=1;i<=n+m;++i) lct.size[i]=i;
    for (int i=1;i<=m;++i)
    {
        if (!opt[i][1])
        {
            int id=eid[opt[i][0]],x=opt[id][0],y=opt[id][1];
            if (inset[id]) inset[id]=0,--setsize;
            else if (ontree[id]) ontree[id]=0,lct.cut(x,id+n),lct.cut(y,id+n);
        }
        else
        {
            int x=opt[i][0],y=opt[i][1];
            if (lct.is_connect(x,y))
            {
                int len,id;
                lct.query(x,y,len,id);
                if (en[id]<en[i])
                {
                    ontree[id]=0,lct.cut(opt[id][0],id+n),lct.cut(opt[id][1],id+n);
                    ontree[i]=1,lct.id[i+n]=i,lct.link(x,i+n),lct.link(y,i+n);
                }else id=i;
                if (!((len-1>>1)&1)) inset[id]=1,++setsize;
            }else ontree[i]=1,lct.link(x,i+n),lct.link(y,i+n);
        }
        printf(setsize?"NO\n":"YES\n");
    }
}

int main()
{
    freopen("graph.in","r",stdin),freopen("graph.out","w",stdout);
    n=read(),m=read(),en[0]=m+2;
    for (int i=1,tp,x,y;i<=m;++i)
    {
        tp=read(),opt[i][0]=x=read()+1;
        if (!tp) en[i]=i;
        else
        {
            opt[i][1]=y=read()+1;
            st[i]=i,en[i]=m+1,eid[++etot]=i;
        }
    }
    process();
    fclose(stdin),fclose(stdout);
    return 0;
}