【bzoj4066】简单题

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