bzoj 4025: 二分图 (分治+图论)

我们注意到如果出现奇环就爆炸，所以我们维护并查集即可，但是如果要连u，v我们考虑u、v到根的距离分别为x，y如果（x+y）%2==0那么根直接合并，但是如果是奇数那么久新建一个点变成x和y的父亲，注意用深度维护平衡。

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
#define xx first
#define yy second
#define mp(a,b) make_pair(a,b)
using namespace std;
typedef long long LL;
typedef double db;
const int maxn=200005;
int ans[maxn];
struct pi{
	int u,v;
	int l,r;
}pp;
int fa[maxn];
int dep[maxn];
typedef pair<int ,int>p1;
int cnt;
void divide(int l,int r,vector<pi>g){
	int p=g.size();
	vector<p1>gg;
	gg.clear();
	int flag=0;
    int pre=cnt;
	for(int i=0;i<p;i++){
		if(g[i].l<=l&&g[i].r>=r){
			int x=g[i].u;
			int x1=0;
			int y=g[i].v;
			int y1=0;
			while(fa[x]!=x){
				x=fa[x];
				x1++;
			}
			while(fa[y]!=y){
				y=fa[y];
				y1++;
			}
			gg.push_back(mp(x,dep[x]));
			gg.push_back(mp(y,dep[y]));
			if(x==y){
				if((x1+y1)%2==0)
				flag=1;
			}
			else{
                if((x1+y1)%2==0){
                    if(dep[x]<dep[y]){
                        swap(x,y);
                    }
                    fa[y]=x;
                    dep[x]=max(dep[x],dep[y]+1);
                }
                else{
                    ++cnt;
                    fa[x]=cnt;
                    fa[y]=cnt;
                    fa[cnt]=cnt;
                    dep[cnt]=max(dep[x],dep[y]);
                }
			}
		}
	}
	if(flag){
		for(int i=l;i<=r;i++) ans[i]=1;
	}
	if(l==r){
		p=gg.size();
		for(int i=p-1;i>=0;i--){
			fa[gg[i].xx]=gg[i].xx;
			dep[gg[i].xx]=gg[i].yy;
		}
		gg.clear();
        cnt=pre;
		return ;
	}
	int mid=(l+r)/2;
	vector<pi>g1,g2;
	g1.clear();
	g2.clear();
	for(int i=0;i<p;i++){
		if(g[i].l<=l&&g[i].r>=r) continue;
		if(g[i].l<=mid) g1.push_back(g[i]);
		if(g[i].r>mid) g2.push_back(g[i]);
	}
	divide(l,mid,g1);
	divide(mid+1,r,g2);
	p=gg.size();
	for(int i=p-1;i>=0;i--){	
		fa[gg[i].xx]=gg[i].xx;
		dep[gg[i].xx]=gg[i].yy;
	}
	gg.clear();
    cnt=pre;
}
int main()
{
   int n,m,t;
   cin>>n>>m>>t;
   for(int i=1;i<=2*n;i++){
	   fa[i]=i;
	   dep[i]=1;
   }
   vector<pi>g;
   g.clear();
	for(int i=0;i<m;i++){
		scanf("%d%d%d%d",&pp.u,&pp.v,&pp.l,&pp.r);
		if(pp.l==pp.r) continue;
		pp.r--;
		g.push_back(pp);
	}
    cnt=n;
	divide(0,t-1,g);
   	for(int i=0;i<t;i++){
		if(ans[i]) printf("No\n");
		else printf("Yes\n");
	}
	return 0;
}