題目傳送門

題目大意：

　　二維平面上有q次操作，每次操作可以是添加一個點，也可以是添加一個矩形，問每次操作後，有多少 點-矩形 這樣的pair，pair的條件是點被矩形覆蓋（邊緣覆蓋也算）。

思路：

　　cdq分治，由於加點和加矩形都既是修改操作又是查詢操作，而且兩種方式完全不一樣，所以用兩部分cdq來寫。

　　先將矩形拆成四個點，並且向左下角擴展一個單元，左下角和右上角的點的權值賦為1，左上角和右下角賦為-1。

　　對於加矩形的操作，遇到加的點則修改樹狀數組，遇到矩形的點查詢小於這個矩形的值，並且乘以這個矩形點的權值。

　　對於加點的操作，由於往左下角擴展了，所以應該按x從右往左處理，碰到一個矩形點，更新樹狀數組，碰到加的點，則ans+=sum（m）-sum（y-1），m是y的上界，因為sum(m)是剛好抵消的情況（等於0），而sum（y-1）則代表了有幾個矩形的右下角在這個點的下方。

　　（還是看代碼比較好寫，註意樹狀數組的上界，sum（m）不要作死的用sum（m+1）代替，因為這個自閉了）。

#include<bits/stdc++.h> 
using namespace std;
const int maxn=100010;
int q,tot,m,brr[maxn<<2],maxd;
long long ans[maxn],res[maxn<<2];
struct node{
    int id,x,y,type,val;
}a[maxn<<2],b[maxn<<2];
inline void Hash(){
    sort(brr
 
+1,brr+1+m);
    m=unique(brr+1,brr+1+m)-brr-1;
    for(int i=1;i<=tot;i++)
    {
        a[i].x=lower_bound(brr+1,brr+1+m,a[i].x)-brr;
        a[i].y=lower_bound(brr+1,brr+1+m,a[i].y)-brr;
        maxd=max(maxd,a[i].y);
    }
}
inline void add(int x,int v){
    while(x<=m){res[x]+=v, x+= x & (-x);}
}
inline 
 
long long  sum(int x){
    long long tu=0;
    while(x>0){tu+=res[x], x -= x & (-x);}
    return tu;
}
inline bool cmpx(const node &a,const node &b){
    if(a.x!=b.x)return a.x<b.x;
    return a.id<b.id;
}
inline bool cmpx2(const node &a,const node &b){
    if(a.x!=b.x)return a.x>b.x;
    return a.id<b.id;
}

inline void cdq1(int l,int r){//加點 
    if(l==r)return ;
    
    int mid=l+((r-l)>>1);
//    printf("l:%d r:%d  mid:%d\n",l,r,mid);
    cdq1(l,mid),cdq1(mid+1,r);
    int cc=0;
    for(int i=l;i<=mid;i++)
    {
        b[++cc]=a[i];
        b[cc].id=0;
    }
    for(int i=mid+1;i<=r;i++)
    {
        b[++cc]=a[i];
    }
    sort(b+1,b+1+cc,cmpx2);
    for(int i=1;i<=cc;i++)
    {
        if(b[i].id==0){
            if(b[i].type==1)continue;
            add(b[i].y,b[i].val);
        //    printf("y:%d  val:%d\n",b[i].y,b[i].val);
        }else{
            if(b[i].type==2)continue;
            ans[b[i].id]+=sum(m)-sum(b[i].y-1);//**sum(m+1) 
        }
    }
    for(int i=1;i<=cc;i++)
    {
        if(b[i].id==0&&b[i].type==2)add(b[i].y,-b[i].val);
    }
}

inline void cdq2(int l,int r){//加矩形 
    if(l==r)return ;
    int mid=l+((r-l)>>1);
    cdq2(l,mid),cdq2(mid+1,r);
    int cc=0;
    for(int i=l;i<=mid;i++)
    {
        b[++cc]=a[i];
        b[cc].id=0;
    }
    for(int i=mid+1;i<=r;i++)
    {
        b[++cc]=a[i];
    }
    sort(b+1,b+1+cc,cmpx);
    for(int i=1;i<=cc;i++)
    {
        if(b[i].id==0){
            if(b[i].type==2)continue;
            add(b[i].y,1);
        }else{
            if(b[i].type==1)continue;
            ans[b[i].id]+=sum(b[i].y)*b[i].val;
        }
    }
    for(int i=1;i<=cc;i++)
    {
        if(b[i].id==0&&b[i].type==1)add(b[i].y,-1);
    }
}

int main(){
    scanf("%d",&q);
    tot=0;
    for(int i=1;i<=q;i++)
    {
        int type,x1,y1,x2,y2;
        scanf("%d",&type);
        a[++tot].type=type;
        a[tot].id=i;
        if(a[tot].type==1)
        {
            scanf("%d%d",&x1,&y1);
            a[tot].x=x1,a[tot].y=y1;
            brr[++m]=x1,brr[++m]=y1;
        }
        else
        {
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            x1--,y1--;//tuo zhan ju xin
            brr[++m]=x1,brr[++m]=x2,brr[++m]=y1,brr[++m]=y2;
            a[tot].x=x1,a[tot].y=y1,a[tot].val=1;
            a[++tot].x=x1,a[tot].y=y2,a[tot].val=-1,a[tot].type=type,a[tot].id=i;
            a[++tot].x=x2,a[tot].y=y2,a[tot].val=1,a[tot].type=type,a[tot].id=i;
            a[++tot].x=x2,a[tot].y=y1,a[tot].val=-1,a[tot].type=type,a[tot].id=i;
        }
    }
    Hash();

    cdq2(1,tot); 
    cdq1(1,tot);
    
    for(int i=1;i<=q;i++)
    {
        ans[i]=ans[i]+ans[i-1]; 
    printf("%lld\n",ans[i]);
    }

}

gym101964G Matrix Queries seerc2018k題 cdq分治