BZOJ4025二分图题解（线段树分治+并查集）

题目：BZOJ4025.
题目大意：给定一张$n$个点$m$条边的无向图，以及每条边的出现时间段$[l_i,r_i]$，问每个时间点的图是否是二分图.
$1\le n\le 10^5,1\le m\le 2\ast 10^5$，结束时间$\le 10^5$.

LCT做法.

考虑二分图的性质，想到一张图为二分图的条件为没有奇环，往这个方向往下想.

考虑只有加边的情况，发现这个问题可以很容易用带权并查集来解决，具体见NOIP2010关押罪犯.

然后重新看题，发现这个问题中边是在某个时间段出现的，这提示我们使用线段树分治来解决.

具体就是对时间建立线段树，在每个节点上开一个vector记录当前区间会出现的边，把每一条边覆盖上去后直接dfs一遍这棵树，每次到一个节点就加入贡献，退出就撤销贡献.

发现这个问题现在只有加边和撤销了，加边已经可以做了，但是撤销怎么办？

想到可以让并查集用启发式合并而不用路径压缩，就可以支持撤销了.

时间复杂度$O(n{\log }^{2}n)$.

代码如下：

#include<bits/stdc++.h>
  using namespace std;

#define Abigail inline void
typedef long long LL;

const int N=100000,M=200000;

struct side{
  int x,y,s,t;
}e[M+9];
int ans[N+9],n,m,t;
vector<int>tr[N*4+9];

void Add_side(int L,int R,int id,int l,int r,int k=1){
  if (L>R) return;
  if (L==l&&R==r) {tr[k].push_back(id);return;}
  int mid=l+r>>1;
  if (R<=mid) Add_side(L,R,id,l,mid,k<<1);
  else if (L>mid) Add_side(L,R,id,mid+1,r,k<<1|1);
    else Add_side(L,mid,id,l,mid,k<<1),Add_side(mid+1,R,id,mid+1,r,k<<1|1);
}

int fa[N+9],d[N+9],siz[N+9];
int st[M+9],ts,flag[M+9],cnt,sx[M+9],sy[M+9];

void Insert(int id){
  int x=e[id].x,y=e[id].y,fx=x,fy=y,cx=d[x],cy=d[y];
  while (fx^fa[fx]) fx=fa[fx],cx^=d[fx];
  while (fy^fa[fy]) fy=fa[fy],cy^=d[fy];
  if (fx==fy){
    if (cx==cy) flag[id]=3,++cnt;
    else flag[id]=2;
  }else{
  	flag[id]=1;
  	if (siz[fx]<siz[fy]) swap(fx,fy),swap(cx,cy);
    siz[fx]+=siz[fy];d[fy]^=cx^cy^1;fa[fy]=fx;
  }
  st[++ts]=id;sx[ts]=fx;sy[ts]=fy;
}

void Erase(){
  int id=st[ts],x=sx[ts],y=sy[ts],tmp=flag[id];--ts;
  flag[id]=0;
  if (tmp==3) --cnt;
  if (tmp^1) return;
  fa[y]=y;d[y]=0;siz[x]-=siz[y];
}

void dfs(int l,int r,int k){
  int len=tr[k].size();
  for (int i=0;i<len;++i)
    Insert(tr[k][i]);
  int mid=l+r>>1;
  if (l==r) ans[l]=cnt;
  else dfs(l,mid,k<<1),dfs(mid+1,r,k<<1|1);
  for (int i=0;i<len;++i)
    Erase();
}

Abigail into(){
  scanf("%d%d%d",&n,&m,&t);
  for (int i=1;i<=m;++i){
  	scanf("%d%d%d%d",&e[i].x,&e[i].y,&e[i].s,&e[i].t);
  	Add_side(e[i].s+1,e[i].t,i,1,t);
  }
}

Abigail work(){
  for (int i=1;i<=n;++i) fa[i]=i,siz[i]=1;
  dfs(1,t,1);
}

Abigail outo(){
  for (int i=1;i<=t;++i)
    puts(ans[i]==0?"Yes":"No");
}

int main(){
  into();
  work();
  outo();
  return 0;
}