JZOJ4769. graph

最新推荐文章于 2019-03-20 21:25:51 发布

Akakii

最新推荐文章于 2019-03-20 21:25:51 发布

阅读量588

点赞数 1

分类专栏： CDQ分治

本文链接：https://blog.csdn.net/Akak__ii/article/details/52503306

版权

题解同时被 3 个专栏收录

161 篇文章 1 订阅

订阅专栏

-----------分治-----------

15 篇文章 0 订阅

订阅专栏

----------图论----------

14 篇文章 0 订阅

订阅专栏

题目大意

给定 $n$ 个点和 $m$ 个操作，每个操作是增加一条边 $(u,v)$ 或者删除之前的一条边。
询问每次操作后的图是否是二分图。

Data Constraint
$n,m\leq300000$

题解

这题是一个比较经典的动态二分图问题。做法挺多，这里讲一种分治的做法(还可以LCT？)。
首先，得到了加边与删边操作就等于我们知道了一条边存在的时间段。
考虑按时间分治。假设当前我们分治的时间段是 $[L,R]$ ，设在这个时段出现过的边集合为 $E$ 。如果存在 $(u,v)\in E$ 存在的时间恰好为 $[L,R]$ 我们就把这条边加入并查集(用并查集维护两点间的距离的奇偶性)，如果此时出现了奇环，那就表明 $[L,R]$ 的图都不是二分图。然后就把边分成两组递归下去求解。

时间复杂度： $O(mlog^2m)$

SRC

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
using namespace std ;

#define N 300000 + 10
#define VEDGE vector < Edge >
struct Edge {
    int u , v ;
    int st , ed ;
} ;

VEDGE E ;

int S[N] ;
int fa[N] , Len[N] , Rank[N] ;
int Ans[N] ;
int n , m , top , Cnt ;

int Get( int x ) { return fa[x] == x ? x : Get(fa[x]) ; }

inline int Find( int x ) {
    int ret = 0 ;
    while ( fa[x] != x ) {
        ret ^= Len[x] ;
        x = fa[x] ;
    }
    return ret ;
}

inline void ReSet( int Goal ) {
    while ( top > Goal ) {
        if ( S[top] < 0 ) {
            S[top] = -S[top] ;
            Rank[S[top]] -- ;
        } else {
            Len[S[top]] = 0 ;
            fa[S[top]] = S[top] ;
        }
        S[top] = 0 ;
        top -- ;
    }
}

void DIV( int l , int r , VEDGE G ) {
    if ( l > r ) return ;
    VEDGE LE , RE ;
    LE.clear() , RE.clear() ;
    int mid = (l + r) / 2 , Origin = top , Num = G.size() ;
    for (int i = 0 ; i < Num ; i ++ ) {
        Edge e = G[i] ;
        if ( e.st == l && e.ed == r ) {
            int fx = Get(e.u) ;
            int fy = Get(e.v) ;
            if ( fx == fy ) {
                if ( (Find(e.u) ^ Find(e.v) ^ 1) ) {
                    ReSet( Origin ) ;
                    return ;
                }
                continue ;
            }
            if ( Rank[fx] > Rank[fy] ) swap( fx , fy ) ;
            Len[fx] = (Find(e.u) ^ Find(e.v) ^ 1) ;
            fa[fx] = fy ;
            S[++top] = fx ;
            if ( Rank[fx] == Rank[fy] ) Rank[fy] ++ , S[++top] = -fy ;
            continue ;
        }
        if ( e.ed <= mid ) LE.push_back(e) ;
        else if ( e.st > mid ) RE.push_back(e) ;
        else {
            int oed = e.ed ;
            e.ed = mid ;
            LE.push_back(e) ;
            e.ed = oed ;
            e.st = mid + 1 ;
            RE.push_back(e) ;
        }
    }
    if ( l == r ) { Ans[l] = 1 ; return ; }
    DIV( l , mid , LE ) ;
    DIV( mid + 1 , r , RE ) ;
    ReSet( Origin ) ;
}

int main() {
    freopen( "graph.in" , "r" , stdin ) ;
    freopen( "graph.out" , "w" , stdout ) ;
    scanf( "%d%d" , &n , &m ) ;
    for (int i = 1 ; i <= n ; i ++ ) fa[i] = i ;
    for (int i = 1 ; i <= m ; i ++ ) {
        int op ;
        scanf( "%d" , &op ) ;
        if ( op ) {
            int a , b ;
            scanf( "%d%d" , &a , &b ) ;
            Edge e ;
            e.u = a + 1 ;
            e.v = b + 1 ;
            e.st = i ;
            e.ed = m ;
            E.push_back(e) ;
        } else {
            int a ;
            scanf( "%d" , &a ) ;
            E[a].ed = i - 1 ;
        }
    }
    DIV( 1 , m , E ) ;
    for (int i = 1 ; i <= m ; i ++ ) {
        if ( Ans[i] ) printf( "YES\n" ) ;
        else printf( "NO\n" ) ;
    }
    return 0 ;
}