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;
}