BZOJ 1176/2683 Mokia (三维偏序CDQ+树状数组)

题目大意：

洛谷传送门

三维偏序裸题。。

每次操作都看成一个三元组$<x,y,t>$，表示$x,y$坐标和操作时间$t $

询问操作拆成$4$个容斥

接下来就是$CDQ$了，外层按t排序，回溯时按$x$排序，用树状数组处理$y$这一维即可

 #include <vector>
 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #define N1 201000
 #define M1 2001000
 #define ll long long
 #define dd double
 #define inf 0x3f3f3f3f3f3f3f3fll
 using namespace std;

 int gint()
 {
     int ret=,fh=;char c=getchar();
     while(c<''||c>''){if(c=='-')fh=-;c=getchar();}
     while(c>=''&&c<=''){ret=ret*+c-'';c=getchar();}
     return ret*fh;
 }
 int n,m,ma;
 struct BIT{
 int s[M1];
 void update(int x,int w){ for(int i=x;i<=ma;i+=(i&(-i))) s[i]+=w; }
 int query(int x){ int ans=; for(int i=x;i;i-=(i&(-i))) ans+=s[i]; return ans; }
 void clr(int x){ for(int i=x;i<=ma;i+=(i&(-i))) s[i]=;}
 }s;
 struct node{int x,y,w,t,ans,p;}a[N1],tmp[N1];
 int cmp1(node s1,node s2){ if(s1.x!=s2.x) return s1.x<s2.x; return s1.y<s2.y; }
 int cmp2(node s1,node s2){ if(s1.x!=s2.x) return s1.x<s2.x; return s1.y<=s2.y; }
 int que[N1],tl;

 void CDQ(int L,int R)
 {
     if(R-L<=) return;
     int i,j,k,M=(L+R)>>;
     for(i=L,j=M,k=L;k<R;k++)
     {
         if(a[k].t<M) tmp[i++]=a[k];
         else tmp[j++]=a[k];
     }
     for(k=L;k<R;k++) a[k]=tmp[k];
     CDQ(L,M);
     for(i=L,j=M;i<M&&j<R;)
     {
         if(cmp2(a[i],a[j])){ if(!a[i].p) { s.update(a[i].y,a[i].w); que[++tl]=i; } i++; }
         else{ if(a[j].p) a[j].ans+=s.query(a[j].y); j++; }
     }
     while(i<M){ i++; }
     while(j<R){ a[j].ans+=s.query(a[j].y); j++; }
     while(tl){ s.clr(a[que[tl--]].y); }
     CDQ(M,R);
     for(i=L,j=M,k=;i<M&&j<R;)
     {
         if(cmp2(a[i],a[j])) tmp[++k]=a[i++];
         else tmp[++k]=a[j++];
     }
     while(i<M){ tmp[++k]=a[i++]; }
     while(j<R){ tmp[++k]=a[j++]; }
     for(k=L;k<R;k++) a[k]=tmp[k-L+];
 }
 int f[N1];

 int main()
 {
     int i,op,xx1,xx2,yy1,yy2,w; op=gint(); ma=gint();
     while(scanf("%d",&op)&&op!=)
     {
         if(op==) { n++; a[n].x=gint(); a[n].y=gint(); a[n].w=gint(); a[n].t=n;}
         if(op==)
         {
             xx1=gint(); yy1=gint(); xx2=gint(); yy2=gint(); m++;
             n++; a[n].x=xx1-; a[n].y=yy1-; a[n].t=n; a[n].w=m; a[n].p=;
             n++; a[n].x=xx2; a[n].y=yy1-; a[n].t=n; a[n].w=m; a[n].p=-;
             n++; a[n].x=xx1-; a[n].y=yy2; a[n].t=n; a[n].w=m; a[n].p=-;
             n++; a[n].x=xx2; a[n].y=yy2; a[n].t=n; a[n].w=m; a[n].p=;
         }
     }
     sort(a+,a+n+,cmp1);
     CDQ(,n+);
     for(i=;i<=n;i++)
     {
         if(!a[i].w) continue;
         f[a[i].w]+=a[i].p*a[i].ans;
     }
     for(i=;i<=m;i++) printf("%d\n",f[i]);
     return ;
 }