Description
你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作:
命令 参数限制 内容
1 x y A 1<=x,y<=N,A是正整数 将格子x,y里的数字加上A
2 x1 y1 x2 y2 1<=x1<= x2<=N
1<=y1<= y2<=N
输出x1 y1 x2 y2这个矩形内的数字和
3 无 终止程序
Input
输入文件第一行一个正整数N。
接下来每行一个操作。每条命令除第一个数字之外,
均要异或上一次输出的答案last_ans,初始时last_ans=0。
Output
对于每个2操作,输出一个对应的答案。
Sample Input
4
1 2 3 3
2 1 1 3 3
1 1 1 1
2 1 1 0 7
3
Sample Output
3
5
HINT
数据规模和约定
1<=N<=500000,操作数不超过200000个,内存限制20M,保证答案在int范围内并且解码之后数据仍合法。
样例解释见OJ2683
新加数据一组,但未重测—-2015.05.24
题解
kd-tree
加入多少次后暴力重构一下。(还是跑的很慢)
代码
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define inf 1000000
#define ll long long
using namespace std;
int Y1,n,opt,x1,x2,y2,cnt,tot,D,num,rt;
ll lastans;
ll read()
{
ll 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;
}
struct P
{
int d[2],mx[2],mn[2],v,l,r;
ll sum;
}p[200005],t[200005];
bool operator<(P a,P b)
{
return a.d[D]<b.d[D];
}
void update(int k)
{
int l=t[k].l,r=t[k].r;
for(int i=0;i<2;i++)
{
t[k].mn[i]=t[k].mx[i]=t[k].d[i];
if(l)t[k].mn[i]=min(t[k].mn[i],t[l].mn[i]);
if(l)t[k].mx[i]=max(t[k].mx[i],t[l].mx[i]);
if(r)t[k].mn[i]=min(t[k].mn[i],t[r].mn[i]);
if(r)t[k].mx[i]=max(t[k].mx[i],t[r].mx[i]);
}
t[k].sum=t[k].v+t[l].sum+t[r].sum;
}
int build(int l,int r,int now)
{
if (l>r) return 0;
D=now;
int mid=(l+r)>>1;
nth_element(p+l,p+mid,p+r+1);
t[mid]=p[mid];
t[mid].l=build(l,mid-1,now^1);
t[mid].r=build(mid+1,r,now^1);
update(mid);
return mid;
}
void insert(int &k,P a,int now)
{
if (!k)
{
k=++tot;
t[k].d[0]=t[k].mn[0]=t[k].mx[0]=a.d[0];
t[k].d[1]=t[k].mn[1]=t[k].mx[1]=a.d[1];
}
if (a.d[0]==t[k].d[0]&&a.d[1]==t[k].d[1])
{
t[k].v+=a.v;
t[k].sum+=a.v;
return;
}
if (a.d[now]<t[k].d[now]) insert(t[k].l,a,now^1);
else insert(t[k].r,a,now^1);
update(k);
}
bool in(int x1,int y1,int x2,int y2,int X1,int Y1,int X2,int Y2)
{
return x1<=X1&&X2<=x2&&y1<=Y1&&Y2<=y2;
}
bool out(int x1,int y1,int x2,int y2,int X1,int Y1,int X2,int Y2)
{
return x1>X2||x2<X1||y1>Y2||y2<Y1;
}
ll query(int k,int x1,int y1,int x2,int y2)
{
if(!k)return 0;
ll tmp=0;
if(in(x1,y1,x2,y2,t[k].mn[0],t[k].mn[1],t[k].mx[0],t[k].mx[1]))return t[k].sum;
if(out(x1,y1,x2,y2,t[k].mn[0],t[k].mn[1],t[k].mx[0],t[k].mx[1]))return 0;
if(in(x1,y1,x2,y2,t[k].d[0],t[k].d[1],t[k].d[0],t[k].d[1]))tmp+=t[k].v;
tmp+=query(t[k].l,x1,y1,x2,y2)+query(t[k].r,x1,y1,x2,y2);
return tmp;
}
int main()
{
n=read();opt=read();
while (opt!=3)
{
if (opt==1)
{
p[++cnt].d[0]=read()^lastans;
p[cnt].d[1]=read()^lastans;
p[cnt].v=read()^lastans;
num++;insert(rt,p[cnt],0);
if (num==10000)
{
num=0;
for (int i=1;i<=tot;i++)
p[i]=t[i];
rt=build(1,tot,0);
}
}
else
{
x1=read()^lastans;Y1=read()^lastans;
x2=read()^lastans;y2=read()^lastans;
lastans=query(rt,x1,Y1,x2,y2);
printf("%lld\n",lastans);
}
opt=read();
}
return 0;
}