一道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
",ans);
ans=0,cnt=0;
}
}
return 0;
} 1176: [Balkan2007]Mokia
吐槽一句為什麼都是權限題呀。。