P5787 二分图 /【模板】线段树分治（线段树分治、并查集）

关于什么是合理的实现

解析

本题把并查集写在了题面上
然而，我却一直沉浸在一个及其通用的判断二分图的方法中：

一个图是二分图的充要条件是它没有奇环

怎么维护这个玩意？带权并查集！
怎么套线段树分治？可持久化！
就这样，口胡完一个写起来几乎不可写的《正解》后，我却没有勇气把它实现…
因为可持久化的带权并查集再套个线段树…这玩意一听就没200行写不来啊…
而且最关键的是，本大聪明又动了动装着一片大海的脑瓜，发现这个东西似乎是3log了…
然后…就恬不知耻的点进了题解
然后…看完题解…就无地自容的点了出来…

但U1S1这确实令人有些耳目一新

首先回归二分图的本质，就是有连边的两点不能在同一集合
所以直接拿状态并查集就可以很轻松的维护了
这就已经把带权给去掉了
那怎么可持久化呢？
不必使用恶心的要死的主席树，可以直接维护一个栈存储并查集的合并操作，在递归结束时撤销
为了方便撤销，我们需要保留树的形态，因此不能路径压缩
但是本来主席树实现的可持久化也不行啊
这样的复杂度就是一个log的！
再加上线段树分治的log，总复杂度$nlog{n}^2$

以后在一些需要撤销并查集操作的地方都可以借鉴这个思路
虽然 （似乎是） 做不了可持久化并查集，但在绝大多数的时候都可以代替
不仅少个log，代码也好写的多

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define ull unsigned long long
const int N=2e5+100;
ll read() {
	ll x=0,f=1;char c=getchar();
	while(!isdigit(c)) {if(c=='-')f=-1;c=getchar();}
	while(isdigit(c)) {x=x*10+(c^48);c=getchar();}
	return x*f;
}

int n,m,k;
int fa[N],siz[N];
int find(int x){
	return x==fa[x]?x:find(fa[x]);
}
struct node{
	int x,y;
};
node zhan[N];
int top;
vector<node>v[N<<2];
bool vis[N];
#define mid ((l+r)>>1)
#define ls (k<<1)
#define rs (k<<1|1)
void add(int k,int l,int r,int x,int y,node o){
	if(x<=l&&r<=y){
		int xx=o.x,yy=o.y;
		//printf("xx=%d yy=%d\n",xx,yy);
		v[k].push_back((node){xx+n,yy});
		v[k].push_back((node){xx,yy+n});
		return;
	}
	if(x<=mid) add(ls,l,mid,x,y,o);
	if(y>mid) add(rs,mid+1,r,x,y,o);
	return;
}
int tp[N<<2];
#define oth(a) (a>n?a-n:a+n)
void give(int k,int l,int r){
	if(l==r){
		vis[l]=1;return;
	}
	give(ls,l,mid);
	give(rs,mid+1,r);
	return;
}
void work(int k,int l,int r){
	//printf("k=%d (%d %d)\n",k,l,r);
	tp[k]=top;int flag=0;
	for(int j=0,o=v[k].size();j<o;j++){
		int x=v[k][j].x,y=v[k][j].y;
		x=find(x);y=find(y);
		//printf("  merge:x=%d y=%d\n",x,y);
		if(x==oth(y)||y==oth(x)){
			//printf("  give!\n");
			give(k,l,r);flag=1;break;
		}
		if(siz[x]>siz[y]) swap(x,y);
		fa[x]=y;siz[y]+=siz[x];
		zhan[++top]=(node){x,y};
	}
	if(l<r&&!flag){
		work(ls,l,mid);work(rs,mid+1,r);
	}
	while(top!=tp[k]){
		int x=zhan[top].x,y=zhan[top].y;top--;
		fa[x]=x;siz[y]-=siz[x];
	}
	return;
}
int main(){
	n=read();m=read();k=read();
	for(int i=1;i<=2*n;i++){
		fa[i]=i;siz[i]=1;
	}
	for(int i=1;i<=m;i++){
		int x=read(),y=read(),l=read(),r=read();
		l++;
		add(1,1,k,l,r,(node){x,y});
	}
	work(1,1,k);
	for(int i=1;i<=k;i++){
		if(vis[i]) printf("No\n");
		else printf("Yes\n");
	}
	return 0;
}
/*
*/