关于什么是合理的实现
解析
本题把并查集写在了题面上
然而,我却一直沉浸在一个及其通用的判断二分图的方法中:
一个图是二分图的充要条件是它没有奇环
怎么维护这个玩意?带权并查集!
怎么套线段树分治?可持久化!
就这样,口胡完一个写起来几乎不可写的《正解》后,我却没有勇气把它实现…
因为可持久化的带权并查集再套个线段树…这玩意一听就没200行写不来啊…
而且最关键的是,本大聪明又动了动装着一片大海的脑瓜,发现这个东西似乎是3log了…
然后…就恬不知耻的点进了题解
然后…看完题解…就无地自容的点了出来…
但U1S1这确实令人有些耳目一新
首先回归二分图的本质,就是有连边的两点不能在同一集合
所以直接拿状态并查集就可以很轻松的维护了
这就已经把带权给去掉了
那怎么可持久化呢?
不必使用恶心的要死的主席树,可以直接维护一个栈存储并查集的合并操作,在递归结束时撤销
为了方便撤销,我们需要保留树的形态,因此不能路径压缩
但是本来主席树实现的可持久化也不行啊
这样的复杂度就是一个log的!
再加上线段树分治的log,总复杂度nlogn2nlogn^2nlogn2
以后在一些需要撤销并查集操作的地方都可以借鉴这个思路
虽然 (似乎是) 做不了可持久化并查集,但在绝大多数的时候都可以代替
不仅少个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;
}
/*
*/