一道CDQ分治的模板題，然而我De了一上午Bug......

按時間分成左右兩半，按x坐標排序然後把y坐標丟到樹狀數組裏，掃一遍遇到左邊的就add,遇到右邊的query

幾個弱智出了bug的點，

一是先分了左右兩半再排序，保證的是這次的左右是上次沒有計算過的貢獻，

for(int i=l;i<=r;i++) qs[i].k=(i>mid);
    sort(qs+l,qs+r+1,cmp2);

然後時間的先後是因為一開始就是按時間排好序的已經保證了。

二是矩陣的一個經典的套路就是拆成四部分差分查詢。

三是很智障的，樹狀數組要註意y==0的情況

四是比較重要的，比較的時候條件要加充分，之前寫陌上花開也出現了這個問題，一維還是兩維相同就不管了，會爆炸。。。一定要考慮清楚所有會影響的維，都要加入排序條件中

bool cmp2(const node&a,const node&b) {
    return a.x<b.x||(a.x==b.x&&(a.y<b.y||(a.y==b.y&&bo[a.id]<bo[b.id])));
}

還有，數組大小要開夠，在下RE好像大半都是數組問題。。

//Twenty
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include
 
<cstring>
#include<queue>
#include<vector>
#define mid ((l+r)>>1)
using namespace std;
const int maxn=200000+299;
const int maxw=2000000+29;
int tot,s,w,opt,x,y,x2,y2,a,sg[maxw],bo[maxn],ans,cnt;
struct node{
    int x,y,v,id,sum,k;
    node(){sum=0;};
    node(int x,int y,int v,int id):x(x),y(y),v(v),id(id){sum=0
 
;}
}qs[maxn];
bool cmp1(const node&a,const node&b) {
    return a.id<b.id;
}
bool cmp2(const node&a,const node&b) {
    return a.x<b.x||(a.x==b.x&&(a.y<b.y||(a.y==b.y&&bo[a.id]<bo[b.id])));//!!!!!!!
}
void add(int x,int v) {
    for(int i=x;i<=w;i+=(i&(-i))) 
       sg[i]+=v; 
}
int qry(int x) {
    int res=0;
    for(int i=x;i>0;i-=(i&(-i))) 
        res+=sg[i];
    return res;
}
void solve(int l,int r) {
    if(l==r) return ;
    solve(l,mid); solve(mid+1,r);
    for(int i=l;i<=r;i++) qs[i].k=(i>mid);
    sort(qs+l,qs+r+1,cmp2);
    for(int i=l;i<=r;i++) {
        if(bo[qs[i].id]&&qs[i].k) {
            if(qs[i].y) qs[i].sum+=qry(qs[i].y);
        }
        else if(!bo[qs[i].id]&&!qs[i].k) {
            if(qs[i].y ) add(qs[i].y,qs[i].v);
        }
    }
    for(int i=l;i<=r;i++) 
        if(!bo[qs[i].id]&&!qs[i].k) {
            if(qs[i].y) add(qs[i].y,-qs[i].v); 
        }
}
int main()
{
    scanf("%d%d",&s,&w);
    for(;;) {
        scanf("%d",&opt); 
        if(opt==3) break;
        if(opt==1) {
            scanf("%d%d%d",&x,&y,&a);
            qs[++tot]=node(x,y,a,tot);
        }
        else {
            scanf("%d%d%d%d",&x,&y,&x2,&y2);
            qs[++tot]=node(x-1,y-1,1,tot);
            qs[++tot]=node(x-1,y2,-1,tot-1);
            qs[++tot]=node(x2,y-1,-1,tot-2);
            qs[++tot]=node(x2,y2,1,tot-3);
            bo[tot-3]=1;//是詢問 
        }
        
    }
    solve(1,tot);
    sort(qs+1,qs+tot+1,cmp1);
    for(int i=1;i<=tot;i++) 
    if(bo[qs[i].id]){
        cnt++;
        ans+=qs[i].sum*qs[i].v;
        if(cnt==4) {
            printf("%d\n",ans);
            ans=0,cnt=0;
          }  
    }
    return 0;
}

吐槽一句為什麽兩道模板都是權限題呀。。

BZOJ 1176: [Balkan2007]Mokia