论文题,简要说一下做法吧。
用
solve(l,r)
表示对于每一个
query
操作
i
,将
那么对于
solve(l,r)
,依旧是分治处理
solve(l,mid)
和
solve(mid+1,r)
,然后我们考虑区间
[l,mid]
中的点对区间
[mid+1,r]
中的矩形的影响,问题转化为求一些点在一些矩形范围内的权值和,我们将第一维排序,第二维用树状数组维护前缀和,然后把一个query拆为四个操作加加减减求解就行了。
最后再膜CDQ
#include<iostream>
#include<cstdio>
#include<algorithm>
#define N 200005
#define lowbit(i) (i&(-i))
using namespace std;
int s,w,cnt,total;
struct node {int x,y,id,dfn,val,opt;} q[N],nq[N];
int tree[N*10],ans[N];
inline int read()
{
int a=0,f=1; char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
return a*f;
}
inline bool operator<(node a,node b)
{
if (a.x==b.x&&a.y==b.y) return a.opt<b.opt;
return a.x==b.x?a.y<b.y:a.x<b.x;
}
inline void add(int x,int val)
{
for (int i=x;i<=w;i+=lowbit(i)) tree[i]+=val;
}
inline int query(int x)
{
int tmp=0;
for (int i=x;i;i-=lowbit(i)) tmp+=tree[i];
return tmp;
}
void solve(int l,int r)
{
if (l==r) return;
int mid=l+r>>1,p1=l,p2=mid+1;
for (int i=l;i<=r;i++)
{
if (q[i].id<=mid&&!q[i].opt) add(q[i].y,q[i].val);
if (q[i].id>mid&&q[i].opt) ans[q[i].dfn]+=q[i].val*query(q[i].y);
}
for (int i=l;i<=r;i++)
if (q[i].id<=mid&&!q[i].opt) add(q[i].y,-q[i].val);
p1=l; p2=mid+1;
for (int i=l;i<=r;i++)
if (q[i].id<=mid) nq[p1++]=q[i]; else nq[p2++]=q[i];
for (int i=l;i<=r;i++) q[i]=nq[i];
solve(l,mid); solve(mid+1,r);
}
int main()
{
s=read(); w=read();
for (;;)
{
int opt=read();
if (opt==3) break;
if (opt==1) q[++cnt].x=read(),q[cnt].y=read(),q[cnt].val=read(),q[cnt].id=cnt;
else
{
total++;
int x1=read(),y1=read(),x2=read(),y2=read();
q[++cnt].x=x1-1; q[cnt].y=y1-1; q[cnt].id=cnt; q[cnt].dfn=total; q[cnt].val=1; q[cnt].opt=1;
q[++cnt].x=x1-1; q[cnt].y=y2; q[cnt].id=cnt; q[cnt].dfn=total; q[cnt].val=-1; q[cnt].opt=1;
q[++cnt].x=x2; q[cnt].y=y1-1; q[cnt].id=cnt; q[cnt].dfn=total; q[cnt].val=-1; q[cnt].opt=1;
q[++cnt].x=x2; q[cnt].y=y2; q[cnt].id=cnt; q[cnt].dfn=total; q[cnt].val=1; q[cnt].opt=1;
}
}
sort(q+1,q+cnt+1);
solve(1,cnt);
for (int i=1;i<=total;i++) printf("%d\n",ans[i]);
return 0;
}