CQUOJ 24920 覆盖次数 - 强化版（cdq分治）

http://acm.cqu.edu.cn:8888/oj/problem_show.php?pid=24920
这次是加的线段长度不是递增，有删除，然后问你每个线段覆盖的线段个数。
数据范围贼大，还有删除，cdq分治的套路，cdq分治就是两边分治，然后求解前面对后面的影响。
这题就是前面的加入并且没有删除的线段，对后面加入的线段的影响，就是按照一个端点排序，然后另一个端点树状数组查询就行了。还有个就是树状数组初始化的问题，如果每次都memset的话做logn次就TLE了，所以每次查询完之后再来一次，然后手动清除树状数组里的东西，复杂度$O(nlog^2n)$

---

代码：

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")

using namespace std;
#define   MAX           200005
#define   MAXN          6005
#define   maxnode       15
#define   sigma_size    30
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
const double pi    = acos(-1.0);
//const double inf   = 1e18;
const double eps   = 1e-8;
const LL    mod    = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
    char c;
    while((c=getchar())<'0' || c>'9');
    x=c-'0';
    while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

struct que{
    int op,id,x,y,f;
    bool operator < (const que&e)const{
        if(x==e.x) return op<e.op;
        return x>e.x;
    }
}q[MAX],a[MAX];
int ans[MAX];
int x[MAX],y[MAX];
int xx[MAX];
int tt;
int c[MAX];

void add(int x,int d){
    while(x<=tt){
        c[x]+=d;
        x+=lowbit(x);
    }
}

int sum(int x){
    int ret=0;
    while(x){
        ret+=c[x];
        x-=lowbit(x);
    }
    return ret;
}

void cdq(int l,int r){
    int mid=(l+r)/2;
    if(l>=r) return;
    cdq(l,mid);cdq(mid+1,r);
    int sz=0;
    for(int i=l;i<=mid;i++){
        a[sz++]=q[i];
    }
    for(int i=mid+1;i<=r;i++){
        if(q[i].op==0){
            a[sz++]=q[i];
            a[sz-1].op=2;
        }
    }
    sort(a,a+sz);
    for(int i=0;i<sz;i++){
        int pos=lower_bound(xx,xx+tt,a[i].y)-xx+1;
        if(a[i].op==2){
            ans[a[i].id]+=sum(pos);
        }
        else add(pos,a[i].f);
    }
    for(int i=0;i<sz;i++){
        int pos=lower_bound(xx,xx+tt,a[i].y)-xx+1;
        if(a[i].op==2){
            //ans[a[i].id]+=sum(pos);
        }
        else add(pos,-a[i].f);
    }
}

int main(){
    int t;
    cin>>t;
    while(t--){
        int Q;
        cin>>Q;
        int cnt=0;
        int tot=0;
        int num=0;
        mem(ans,0);
        for(int i=1;i<=Q;i++){
            int op;
            scanf("%d",&op);
            if(op==0){
                scanf("%d%d",&x[cnt],&y[cnt]);
                xx[num++]=x[cnt];
                xx[num++]=y[cnt];
                q[++tot]=(que){op,cnt+1,x[cnt],y[cnt],1};
                cnt++;
            }
            else{
                int o;
                scanf("%d",&o);
                q[++tot]=(que){op,cnt,x[o-1],y[o-1],-1};
            }
        }
        sort(xx,xx+num);
        tt=unique(xx,xx+num)-xx;
        cdq(1,tot);
        for(int i=1;i<=cnt;i++){
            printf("%d\n",ans[i]);
        }
    }
    return 0;
}