题意:加入边,删除边,查询当前图是否为二分图
本来想练lct,然后发现了线段树分治的做法,感觉好厉害。
lct做法的核心就是维护删除时间的最大生成树
首先口胡一个分块做法,和hnoi2016第一题类似的偏序关系,一样做。
线段树分治
数据结构题中如果使用对时间cdq分治,要求每个操作独立,不能很好的处理撤销(删除)操作。
采取线段树区间标记的思想
对于一个操作,它的存在时间是\([l,r]\)
我们模仿线段树打标记的过程进行分治,\(cdq(l,r,S)\)表示当前处理时间\([l,r]\),操作集合为\(S\)
如果区间就是当前区间,那么进行操作
否则继续递归
对于本题,用启发式合并 不路径压缩的并查集实现加边和撤销
越卡常越慢是smg
#include#include #include #include #include using namespace std;typedef long long ll;const int N=2e5+5;inline int read(){ char c=getchar(); int x=0,f=1; while(c<'0' || c>'9') {if(c=='-')f=-1; c=getchar();} while(c>='0' && c<='9') {x=x*10+c-'0'; c=getchar();} return x*f;} int n, m, T, u, v, l, r;struct edge { int u, v, l, r; bool operator <(const edge &a) const {return l == a.l ? r < a.r : l < a.l;}};typedef vector meow;meow a; int top;namespace ufs { struct node {int fa, val, size;} t[N]; struct info {int x, y; node a, b;} st[N]; inline void init() {for(int i=1; i<=n; i++) t[i] = (node){i, 0, 1};} inline int find(int x) {while(t[x].fa != x) x = t[x].fa; return x;} inline int dis(int x) { int ans=0; while(t[x].fa != x) ans ^= t[x].val, x = t[x].fa; return ans; } inline void link(int x, int y) { int val = dis(x) ^ dis(y) ^ 1; x = find(x); y = find(y); st[++top] = (info) {x, y, t[x], t[y]}; if(t[x].size > t[y].size) swap(x, y); t[x].fa = y; t[x].val = val; t[y].size += t[x].size; } inline void recov(int bot) { while(top != bot) { info &now = st[top--]; t[now.x] = now.a; t[now.y] = now.b; } }} using namespace ufs; void cdq(int l, int r, meow &a) { int mid = (l+r)>>1, bot = top; meow b, c; for(int i=0; i<(int)a.size(); i++) { edge &now = a[i]; int x = now.u, y = now.v; if(now.l == l && now.r == r) { int p = find(x), q = find(y); if(p == q) { int val = dis(x) ^ dis(y); if(val == 0) { for(int i=l; i<=r; i++) puts("No"); recov(bot); return; } } else link(x, y); } else if(now.r <= mid) b.push_back(now); else if(mid < now.l) c.push_back(now); else b.push_back( (edge){now.u, now.v, now.l, mid} ), c.push_back( (edge){now.u, now.v, mid+1, now.r} ); } if(l == r) puts("Yes"); else cdq(l, mid, b), cdq(mid+1, r, c); recov(bot);} int main() { //freopen("in", "r", stdin); n=read(); m=read(); T=read(); for(int i=1; i<=m; i++) { u=read(), v=read(), l=read()+1, r=read(); if(l > r) continue; a.push_back((edge){u, v, l, r}); } init(); cdq(1, T, a);}