[LCT 树状数组] HDU 5333 Undirected Graph

这个题跟 [LCT 主席树] BZOJ 3514 Codechef MARCH14 GERALD07加强版有异曲同工之妙
考虑询问按右端点排序
对于一条边 边的权值定义为较小的一端的编号
那么显然权值大的边更有利
对于每个右端点 我们加入所有新增的边 维护最大生成树
询问么 就是n-生成树中编号大于等于左端点 因为这样的边会把两个联通块连在一起

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<stack>
using namespace std;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
  return *p1++;
}
inline bool read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1; else if (c==EOF) return 0;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b; return 1;
}

const int N=100005;

namespace BIT{
  int c[N],maxn;
  inline void init(int n){ maxn=n; for (int i=1;i<=n;i++) c[i]=0; }
  inline void add(int x,int r){ for (int i=x;i<=maxn;i+=i&-i) c[i]+=r; }
  inline int sum(int x){ int ret=0; for (int i=x;i;i-=i&-i) ret+=c[i]; return ret; }
  inline int sum(int l,int r){ return sum(r)-sum(l-1); }
}

int us[N*2],vs[N*2],ws[N*3];
inline int Min(int a,int b){
  return ws[a]<ws[b]?a:b;
}

struct node{
  node *f,*ch[2];
  int idx,v; bool rev;
  bool isr() { return !f || (f->ch[0]!=this && f->ch[1]!=this); }
  bool dir() { return f->ch[1]==this; }
  void setc(node *c,int d) { ch[d]=c; if (c) c->f=this; }
  void push() {
    if (rev) { swap(ch[0],ch[1]); if(ch[0])ch[0]->rev^=1; if(ch[1])ch[1]->rev^=1; rev=0;}
  }
  void upd() {
    v=idx; if(ch[0])v=Min(v,ch[0]->v); if(ch[1])v=Min(v,ch[1]->v);
  }
}pool[N*3];
inline void rot(node *x){
  node *p=x->f; bool d=x->dir();
  if (!p->isr()) p->f->setc(x,p->dir()); else x->f=p->f;
  p->setc(x->ch[d^1],d); x->setc(p,d^1); p->upd();
}
stack<node*> sta;
inline void splay(node *x){
  node *q=x; for (;!q->isr();q=q->f) sta.push(q); sta.push(q);
  while (!sta.empty()) sta.top()->push(),sta.pop();
  while (!x->isr())
    if (x->f->isr()) rot(x);
    else if (x->dir()==x->f->dir()) rot(x->f),rot(x);
    else rot(x),rot(x);
  x->upd();
}
void expose(node *x){
  for (node *q=NULL;x;q=x,x=x->f) splay(x),x->ch[1]=q,x->upd();
}
void exert(node *x){ expose(x); splay(x); x->rev^=1; }
void link(node *x,node *y) { exert(y); y->f=x; }
void cut(node *x,node *y) { exert(x); expose(y); splay(y); y->ch[0]=x->f=NULL; y->upd(); }
node *frt(node *x){ expose(x),splay(x); for (;x->push(),x->ch[0];x=x->ch[0]); return x; }
int query(node *x,node *y){ exert(x); expose(y); splay(y); return y->v; }

int n,m;

struct abcd{
  int l,r,idx;
  bool operator < (const abcd &B) const{ return r<B.r; }
}eve[N*2];
int tot,ans[N];

inline void addedge(int i){
  if (frt(pool+m+us[i])==frt(pool+m+vs[i])){
    int p=query(pool+m+us[i],pool+m+vs[i]);
    if (ws[p]<ws[i]){
      cut(pool+p,pool+m+us[p]),cut(pool+p,pool+m+vs[p]); BIT::add(ws[p],-1);
      link(pool+i,pool+m+us[i]),link(pool+i,pool+m+vs[i]); BIT::add(ws[i],1);
    }
  }else
    link(pool+i,pool+m+us[i]),link(pool+i,pool+m+vs[i]),BIT::add(ws[i],1);
}

inline int query(int i){
  return n-BIT::sum(eve[i].l,eve[i].r);
}

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  while (read(n)&&read(m)&&read(tot)){
    for (int i=1;i<=m;i++) { read(eve[i].l),read(eve[i].r); if (eve[i].l>eve[i].r) swap(eve[i].l,eve[i].r); }
    sort(eve+1,eve+m+1);
    for (int i=1;i<=m;i++) us[i]=ws[i]=eve[i].l,vs[i]=eve[i].r;
    for (int i=m+1;i<=m+n;i++) ws[i]=1<<30;
    for (int i=1;i<=m+n;i++) pool[i].idx=pool[i].v=i,pool[i].rev=0,pool[i].f=pool[i].ch[0]=pool[i].ch[1]=NULL;
    for (int i=1;i<=tot;i++) read(eve[i].l),read(eve[i].r),eve[i].idx=i;
    sort(eve+1,eve+tot+1); int pnt1=1,pnt2=1;
    BIT::init(n);
    for (int i=1;i<=n;i++){
      while (pnt1<=m && vs[pnt1]==i)
    addedge(pnt1),pnt1++;
      while (pnt2<=tot && eve[pnt2].r==i)
    ans[eve[pnt2].idx]=query(pnt2),pnt2++;
    }
    for (int i=1;i<=tot;i++) printf("%d\n",ans[i]);
  }
  return 0;
}