bzoj1176 [Balkan2007]Mokia（cdq分治+树状数组）

S据说一直都是0.。。好像是描述有误。。可以直接忽视了。
把每一个询问矩形容斥成4个，使得都是前缀矩形，然后cdq分治就好啦。注意树状数组不能每次O(N)大规模清空，这样复杂度其实是不对的。必须是O(（mid-l+1）*logn)的清空。原来我以前写的都是假的orz

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x7fffffff
#define N 2000010
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int S,n=0,m,tot=0,ans[10010],c[N];
struct oper{
    int op,x,y,id,val;
    friend bool operator<(oper a,oper b){
        return a.x==b.x?a.op<b.op:a.x<b.x;
    }
}q[200010],tmp[200010];
inline void add(int x,int val){
    for(;x<=m;x+=x&-x) c[x]+=val;
}inline int ask(int x){
    int res=0;for(;x;x-=x&-x) res+=c[x];return res;
}
inline void cdq(int l,int r){
    if(l==r) return;
    int mid=l+r>>1;cdq(l,mid);cdq(mid+1,r);
    int p1=l,p2=mid+1,owo=0;
    while(p1<=mid&&p2<=r){
        if(q[p1]<q[p2]){
            if(q[p1].op==1) add(q[p1].y,q[p1].val);
            tmp[++owo]=q[p1++];
        }else{
            if(q[p2].op==2) ans[q[p2].id]+=ask(q[p2].y)*q[p2].val;
            tmp[++owo]=q[p2++];
        }
    }int ed=p1-1;
    while(p1<=mid) tmp[++owo]=q[p1++];
    while(p2<=r){
        if(q[p2].op==2) ans[q[p2].id]+=ask(q[p2].y)*q[p2].val;
        tmp[++owo]=q[p2++];
    }for(int i=l;i<=ed;++i) if(q[i].op==1) add(q[i].y,-q[i].val);
    for(int i=l;i<=r;++i) q[i]=tmp[i-l+1];
}
int main(){
//  freopen("a.in","r",stdin);
    S=read();m=read();
    while(1){
        int op=read();if(op==3) break;int x=read(),y=read();
        if(op==1) q[++n].op=1,q[n].x=x,q[n].y=y,q[n].val=read();
        else{
            int x2=read(),y2=read();++tot;
            q[++n].op=2;q[n].x=x2;q[n].y=y2;q[n].id=tot;q[n].val=1;
            q[++n].op=2;q[n].x=x-1;q[n].y=y2;q[n].id=tot;q[n].val=-1;
            q[++n].op=2;q[n].x=x2;q[n].y=y-1;q[n].id=tot;q[n].val=-1;
            q[++n].op=2;q[n].x=x-1;q[n].y=y-1;q[n].id=tot;q[n].val=1;
        }
    }cdq(1,n);
    for(int i=1;i<=tot;++i) printf("%d\n",ans[i]);
    return 0;
}