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