[BZOJ4537][Hnoi2016][分块]最小公倍数

很神的做法……

第一眼以为是cdq分治加lct（排序a,分治b）敲完lct模板发现处理不了某条边a很大b很小的情况……

把边按a排序，分成B块，第i块的a的权值范围为[Li,Ri]，那么对于询问j，aj在[Li,Ri]的范围内，把这一块以及这一块之前ai<=aj且bi<=bj的边加入图中，用并查集维护、查询。

这样复杂度可以降到$O(n\sqrt n ~log n)$

分块的大小取$\sqrt{nlogn}$时间要快一倍……

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>
#define N 100010

using namespace std;

int n,m,q,block,sz;
char Ans[N];
int f[N],mxa[N],mxb[N],size[N];

struct stp{
  int x,y,a,b,g;
}E[N],Q[N],tmp[N];
struct retuse{
  int g,mxa,mxb,isz;
}r[N<<2];

inline bool cmpa(stp a,stp b){
  if(a.a==b.a) return a.b<b.b;
  return a.a<b.a;
}

inline bool cmpb(stp a,stp b){
  if(a.b==b.b) return a.a<b.a;
  return a.b<b.b;
}

inline void reaD(int &x){
  char c=getchar();x=0;
  for(;c>57||c<48;c=getchar());for(;c>=48&&c<=57;x=x*10+c-48,c=getchar());
}

int find(int x){
  if(x==f[x]) return x;
  return find(f[x]);
}

inline void merge(int x,int y,int a,int b){
  int i=find(x),j=find(y);
  r[++sz].g=i; r[sz].mxa=mxa[i]; r[sz].mxb=mxb[i]; r[sz].isz=size[i];
  r[++sz].g=j; r[sz].mxa=mxa[j]; r[sz].mxb=mxb[j]; r[sz].isz=size[j];
  if(i==j){
    mxa[i]=max(mxa[i],a);
    mxb[i]=max(mxb[i],b);
    return;
  }
  if(size[i]>size[j]) swap(i,j);
  f[i]=j; size[j]+=size[i];
  mxa[j]=max(max(mxa[j],mxa[i]),a);
  mxb[j]=max(max(mxb[j],mxb[i]),b);
}

inline int query(int x,int y,int a,int b){
  int i=find(x),j=find(y);
  if(i!=j) return 0;
  return mxa[i]==a&&mxb[i]==b;
}

inline void Return(){
  for(;sz;--sz){
    f[r[sz].g]=r[sz].g;
    mxa[r[sz].g]=r[sz].mxa;
    mxb[r[sz].g]=r[sz].mxb;
    size[r[sz].g]=r[sz].isz;
  }
}

int main(){
  reaD(n); reaD(m);
  for(int i=1;i<=m;i++)
    reaD(E[i].x),reaD(E[i].y),reaD(E[i].a),reaD(E[i].b);
  reaD(q);
  for(int i=1;i<=q;i++){
    reaD(Q[i].x),reaD(Q[i].y),reaD(Q[i].a),reaD(Q[i].b);
    Q[i].g=i;
  }
  sort(E+1,E+1+m,cmpa);
  block=(int)sqrt(m*20);
  for(int i=1;i<=m;i+=block){
    sort(E+1,E+i,cmpb); int tp=0;
    for(int j=1;j<=n;j++) f[j]=j,mxa[j]=mxb[j]=-1,size[j]=1;
    for(int j=1,k=1;j<=q;j++)
      if(Q[j].a>=E[i].a&&(i+block>m||Q[j].a<E[i+block].a))
    tmp[++tp]=Q[j];
    sort(tmp+1,tmp+1+tp,cmpb);
    for(int j=1,k=1;j<=tp;j++){
      while(k<i&&E[k].b<=tmp[j].b) merge(E[k].x,E[k].y,E[k].a,E[k].b),k++;
      sz=0;
      for(int l=i;l<i+block&&l<=m;l++)
    if(E[l].a<=tmp[j].a&&E[l].b<=tmp[j].b)
      merge(E[l].x,E[l].y,E[l].a,E[l].b);
      Ans[tmp[j].g]=query(tmp[j].x,tmp[j].y,tmp[j].a,tmp[j].b);
      Return();
    }
  }
  for(int i=1;i<=q;i++) puts(Ans[i]?"Yes":"No");
}