应该算是入门题,只不过我很久没做这种矩形题了,所以没有想到一个关键点导致我模型没转化过来。
还是像之前一样把询问和赋值多添加一个维,表示时间,然后添加很显然就不说了,每次直接bit加入就好,问题是查询。
一个子矩阵x1 y1 x2 y2 的答案等于ans[x2][y2]+ans[x1][y1]-ans[x1][y2]-ans[x2][y1]
注意一下边界就好了,然后直接cdq。
%GEOTCBRL
还有那个S是没有意义的,题目错了。
#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=2e6+5;
typedef int ll;
int n;
struct node
{
int id,x,y,v,t;
node(){}
node(int id,int x,int y,int v,int t):id(id),x(x),y(y),v(v),t(t){}
bool operator < (const node a)const
{
if (x==a.x&&y==a.y)return t<a.t;
if (x==a.x)return y<a.y;
return x<a.x;
}
}q[N],tmp[N];
inline int lowbit(int x)
{
return x&(-x);
}
ll val[N];
inline void add(int x,int y)
{
while (x<=n)
{
val[x]+=y;
x+=lowbit(x);
}
}
inline ll ask(int x)
{
ll ret=0;
while (x>0)
{
ret+=val[x];
x-=lowbit(x);
}
return ret;
}
int tot,ques;
ll ans[N];
inline void ins(int x,int y,int v,int t)
{
++tot;
q[tot]=node(tot,x,y,v,t);
}
inline void cdq(int l,int r)
{
if (l==r)return;
int mid=(l+r)>>1;
fo(i,l,r)
{
if (q[i].id<=mid&&!q[i].t)add(q[i].y,q[i].v);
if (q[i].id>mid&&q[i].t)ans[q[i].t]+=ask(q[i].y)*q[i].v;
}
fo(i,l,r)if (q[i].id<=mid&&!q[i].t)add(q[i].y,-q[i].v);
int l1=l,l2=mid+1;
fo(i,l,r)
if (q[i].id<=mid)tmp[l1++]=q[i];
else tmp[l2++]=q[i];
fo(i,l,r)q[i]=tmp[i];
cdq(l,mid);
cdq(mid+1,r);
}
int main()
{
scanf("%d",&n);
int opt;
while (scanf("%d",&opt)&&opt!=3)
{
int x1,y1,x2,y2;
scanf("%d%d%d",&x1,&y1,&x2);
if (opt==1)ins(x1,y1,x2,0);
else
{
++ques;
scanf("%d",&y2);
ins(x1-1,y1-1,1,ques);
ins(x1-1,y2,-1,ques);
ins(x2,y1-1,-1,ques);
ins(x2,y2,1,ques);
}
}
sort(q+1,q+tot+1);
cdq(1,tot);
fo(i,1,ques)printf("%d\n",ans[i]);
return 0;
}