graph

题目大意

维护一个图,若干操作每次加入一条边或删除一条已有边,每次操作后都输出该图是否为二分图。

CDQ

预处理每条边的存在时间。
考虑一个时刻,此时存在的边是所有存在时间区间包含该时刻的。对于每个时刻做的话,就是把所有包含该时刻的时间区间对应边加入并查集。
现在考虑把所有时刻一起做。
solve(l,r)表示解决[l,r]所有询问,先把所有时间区间刚好为[l,r]的加入并查集,对于剩余边,要么分到左区间,要么分到右区间,要么拆成两个时间区间分别下放左右区间。
然后就solve(l,mid)和solve(mid+1,r)
由于退出时需要把这个区间中加入的边都删除,因此可以开一个栈存储加入并查集的边,退出时把该区间加入的边退栈并删除。如何使并查集能删除呢?我们就要按秩合并了。
并查集维护很简单,维护每个节点的颜色。因为是按秩合并,还有删除,所以其实就是每个节点维护一个是否和父亲同色,然后查询一个点的颜色也用一个log的时间向上。
假如在维护并查集时发现非法,黑白染色失败意味着无法形成二分图,那么[l,r]全部不是二分图,可以直接退出。

#include<cstdio>
#include<algorithm>
#include<vector>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
const int maxn=300000+10;
struct dong{
    int l,r,x,y;
};
vector<dong> dl[maxn*4];
dong zlt,xdl;
int rank[maxn],co[maxn],fa[maxn],sta[maxn][4];
int edge[maxn][4];
bool bz[maxn];
int i,j,k,l,t,n,m,tot,top;
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;
}
int getfa(int x){
    while (fa[x]) x=fa[x];
    return x;
}
int getcolor(int x){
    int t=0;
    while (fa[x]){
        t^=co[x];
        x=fa[x];
    }
    return t;
}
void split(int id){
    int x=sta[id][0],y=sta[id][1];
    fa[x]=0;
    co[x]=0;
    rank[x]=sta[id][2];
    rank[y]=sta[id][3];
}
void solve(int p,int l,int r){
    int i,j,x,y,a,b,mid=(l+r)/2,size=dl[p].size(),cnt=0;
    fo(i,0,size-1){
        zlt=dl[p][i];
        if (zlt.l==l&&zlt.r==r){
            x=getfa(zlt.x);y=getfa(zlt.y);
            a=getcolor(zlt.x);b=getcolor(zlt.y);
            if (x==y&&a==b){
                //fo(j,l,r) bz[j]=0;
                while (cnt--) split(top--);
                return;
            }
            else if (x!=y){
                cnt++;
                if (rank[x]>rank[y]) swap(x,y);
                fa[x]=y;
                sta[++top][0]=x;
                sta[top][1]=y;
                sta[top][2]=rank[x];
                sta[top][3]=rank[y];
                if (a==b) co[x]=1;
                if (rank[x]==rank[y]) rank[y]++;
            }
        }
        else if (zlt.r<=mid) dl[p*2].push_back(zlt);
        else if (zlt.l>mid) dl[p*2+1].push_back(zlt);
        else{
            xdl=zlt;
            zlt.r=mid;
            xdl.l=mid+1;
            dl[p*2].push_back(zlt);
            dl[p*2+1].push_back(xdl);
        }
    }
    if (l==r){
        bz[l]=1;
        while (cnt--) split(top--);
        return;
    }
    solve(p*2,l,mid);
    solve(p*2+1,mid+1,r);
    while (cnt--) split(top--);
}
int main(){
    freopen("graph.in","r",stdin);freopen("graph.out","w",stdout);
    n=read();m=read();
    fo(i,1,m){
        t=read();
        if (t){
            j=read();k=read();
            j++;k++;
            edge[++tot][0]=j;
            edge[tot][1]=k;
            edge[tot][2]=i;
            edge[tot][3]=m;
        }
        else{
            j=read();
            j++;
            edge[j][3]=i-1;
        }
    }
    fo(i,1,m){
        zlt.l=edge[i][2];zlt.r=edge[i][3];
        zlt.x=edge[i][0];zlt.y=edge[i][1];
        dl[1].push_back(zlt);
    }
    //fo(i,1,m) bz[i]=1;
    solve(1,1,m);
    fo(i,1,m) if (bz[i]) printf("YES\n");else printf("NO\n");
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值