bzoj 4025 二分图——线段树分治+LCT

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=4025

线段树分治，用 LCT 维护链的长度即可。不过很慢。

正常（更快）的方法应该是线段树分治+并查集（按秩合并，链长可以暴力爬）或者 LCT 维护删除时间最大生成树。就不写了。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define ls Ls[cr]
#define rs Rs[cr]
#define lc c[x][0]
#define rc c[x][1]
#define pb push_back
using namespace std;
int rdn()
{
  int ret=;bool fx=;char ch=getchar();
  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
  return fx?ret:-ret;
}
const int N=1e5+,M=2e5+;
int n,m,fa[N],c[N][],rev[N],siz[N],stk[N],tp;
int tot,Ls[M],Rs[M],top;
struct Node{
  int x,y;
  Node(int x=,int y=):x(x),y(y) {}
}sta[N];
vector<Node> vt[M];

void build(int l,int r,int cr)
{
  if(l==r)return; int mid=l+r>>;
  ls=++tot; build(l,mid,ls);
  rs=++tot; build(mid+,r,rs);
}
void ins(int l,int r,int cr,int L,int R,Node k)
{
  if(l>=L&&r<=R){vt[cr].pb(k);return;}
  int mid=l+r>>;
  if(L<=mid)ins(l,mid,ls,L,R,k);
  if(mid<R)ins(mid+,r,rs,L,R,k);
}

bool isrt(int x){return c[fa[x]][]!=x&&c[fa[x]][]!=x;}
void pshp(int x){siz[x]=siz[lc]+siz[rc]+;}
void Rev(int x){if(rev[x]){rev[x]=;rev[lc]^=;rev[rc]^=;swap(lc,rc);}}
void rotate(int x)
{
  int y=fa[x],z=fa[y],d=(x==c[y][]);
  if(!isrt(y))c[z][y==c[z][]]=x;
  fa[x]=z;
  fa[y]=x; fa[c[x][!d]]=y;
  c[y][d]=c[x][!d]; c[x][!d]=y;
  pshp(y); pshp(x);
}
void splay(int x)
{
  stk[tp=]=x;
  for(int k=x;!isrt(k);k=fa[k])stk[++tp]=fa[k];
  for(int i=tp;i;i--)Rev(stk[i]);
  int y,z;
  while(!isrt(x))
    {
      y=fa[x]; z=fa[y];
      if(!isrt(y))
    ( (x==c[y][])^(y==c[z][]) )?rotate(x):rotate(y);
      rotate(x);
    }
}
void access(int x)
{
  for(int t=;x;splay(x),rc=t,pshp(x),t=x,x=fa[x]);
}
void mkrt(int x)
{
  access(x); splay(x); rev[x]^=;
}
void split(int x,int y)
{
  mkrt(x); access(y); splay(y);
}
bool link(Node k)
{
  int x=k.x, y=k.y; split(x,y);
  int cr=y; while(c[cr][])cr=c[cr][];
  if(cr==x) return (siz[y]&);
  sta[++top]=k; fa[x]=y; return false;
}
void cut(Node k)
{
  int x=k.x, y=k.y; split(x,y);
  c[y][]=fa[x]=; pshp(y);
}
void solve(int l,int r,int cr)
{
  int sz=vt[cr].size(); bool flag=;
  for(int i=;i<sz;i++)
    {
      flag=link(vt[cr][i]); if(flag)break;
    }
  if(flag){for(int i=l;i<=r;i++)puts("No");return;}
  if(l==r){puts("Yes");return;}
  int mid=l+r>>,nw=top;
  solve(l,mid,ls); for(int& i=top;i>nw;i--)cut(sta[i]);
  solve(mid+,r,rs); for(int& i=top;i>nw;i--)cut(sta[i]);
}
int main()
{
  n=rdn();int T=rdn();m=rdn();
  tot=;build(,m,);
  for(int i=;i<=n;i++)siz[i]=;
  for(int i=,u,v,st,en;i<=T;i++)
    {
      u=rdn();v=rdn();st=rdn()+;en=rdn();
      ins(,m,,st,en,Node(u,v));
    }
  solve(,m,); return ;
}