[BZOJ 2683] 简单题

最新推荐文章于 2022-04-29 23:18:56 发布

LPA20020220

最新推荐文章于 2022-04-29 23:18:56 发布

阅读量607

点赞数

分类专栏： CDQ 分治文章标签： CDQ 分治

本文链接：https://blog.csdn.net/lpa20020220/article/details/80363968

版权

分治同时被 2 个专栏收录

8 篇文章 0 订阅

订阅专栏

CDQ

7 篇文章 0 订阅

订阅专栏

BZOJ传送门(权限题)

题目描述

你有一个 $N*N$ 的棋盘，每个格子内有一个整数，初始时的时候全部为0，现在需要维护两种操作：
$1\ x\ y\ A$ ， $1<=x,y<=N$ ， $A$ 是正整数，将格子 $(x,y)$ 里的数字加上 $A$

$2\ x_1\ y_1\ x_2\ y_2$ ， $1\leq x_1\leq x_2\leq N, 1\leq y1\leq y_2 \leq N$ ，输出 $x_1\ y_1\ x_2\ y_2$ 这个矩形内的数字和

$3$ 终止程序

输入输出格式

输入格式

输入文件第一行一个正整数 $N$ 。接下来每行一个操作。

输出格式

对于每个2操作，输出一个对应的答案。

输入输出样例

输入样例

输出样例

3
5

数据范围

$1\leq N\leq 500000$ ,操作数不超过200000个，内存限制20M。对于100%的数据，操作1中的A不超过2000。

解题分析

看到 $500000$ 的数据范围和 $20M$ 就知道二维树状数组 / 线段树凉凉了… 那么我们只好用 $CDQ$ 分治做… 显然对于一个矩形，我们需要查询四次在树状数组上的前缀差分答案即可。

不过有个很迷的问题：博主差分时将所有点X坐标和Y坐标都加1（防止差分至0节点导致爆炸），但居然有两个点WA掉了，而不加1则AC。

代码如下：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cctype>
#include <cstdlib>
#include <cmath>
#define R register
#define IN inline
#define gc getchar()
#define W while
#define MX 500050
#define lowbit(i) (i & -i)
#define ll long long
template <class T>
IN void in (T &x)
{
    bool fu = false; x = 0; R char c = gc;
    W (!isdigit(c)) 
    {
        if(c == '-') fu = true;
        c = gc;
    }
    W (isdigit(c))
    x = (x << 1) + (x << 3) + c - 48, c = gc;
    if(fu) x = -x;
}
struct Que
{//博主将四个端点分为左右两部分， 存为两个询问， 最后相减即可。
    int x, y1, y2, del, typ, id, tim;
    short ord;
}eve[MX << 1], buf[MX << 1];
IN bool operator < (const Que &x, const Que &y)
{return x.x == y.x ? x.typ < y.typ : x.x < y.x;}
struct Res
{
    ll res[3];
}ans[MX];
int q, cnt, tot, arr, limit;
ll get1, get2, tree[MX];
IN void add(R int pos, const int &del)
{W (pos <= limit) tree[pos] += del, pos += lowbit(pos);}
IN void clear(R int pos)
{W (pos <= limit) {tree[pos] = 0, pos += lowbit(pos);}}
IN ll query(R int pos)
{
    R ll ret = 0;
    W (pos) ret += tree[pos], pos -= lowbit(pos);
    return ret;
}
void cdq(const int &lef, const int &rig)
{
    if(lef == rig) return;
    int mid = lef + rig >> 1;
    int lb = lef, rb = mid + 1;
    for (R int i = lef; i <= rig; ++i)//对时间排序
    {
        if(eve[i].tim > mid) buf[rb++] = eve[i];
        else buf[lb++] = eve[i];
    }
    for (R int i = lef; i <= rig; ++i) eve[i] = buf[i];
    lb = lef, rb = mid + 1;
    W (233)
    {
        if(eve[rb].typ ^ 1)
        {
            W (eve[lb].x <= eve[rb].x && lb <= mid)
            {
                if(eve[lb].typ & 1) add(eve[lb].y1, eve[lb].del);
                ++lb;
            }
            get1 = query(eve[rb].y1);
            get2 = query(eve[rb].y2);
            ans[eve[rb].id].res[eve[rb].ord] += get1 - get2;
        }
        ++rb;
        if(rb > rig) break;
    }
    for (R int i = lef; i < lb; ++i) if(eve[i].typ & 1) clear(eve[i].y1);
    cdq(lef, mid), cdq(mid + 1, rig);
}
int main(void)
{
    int a, b, c, d;
    in(limit);
    W (233)
    {
        in(eve[++cnt].typ);
        if(eve[cnt].typ == 3) {--cnt; break;}
        if(eve[cnt].typ ^ 1)
        {
            in(a), in(b), in(c), in(d);

            eve[cnt].ord = 1, eve[cnt].x = a - 1, eve[cnt].y2 = b - 1, eve[cnt].y1 = d, eve[cnt].id = ++tot, eve[cnt].ord = 1, eve[cnt].tim = ++arr;
            eve[++cnt].typ = 2, eve[cnt].x = c, eve[cnt].y2 = b - 1, eve[cnt].y1 = d, eve[cnt].id = tot, eve[cnt].ord = 2, eve[cnt].tim = ++arr;
        }
        else {in(eve[cnt].x), in(eve[cnt].y1), in(eve[cnt].del), eve[cnt].tim = ++arr;}
    }
    std::sort(eve + 1, eve + 1 + cnt);
    cdq(1, cnt);
    for (R int i = 1; i <= tot; ++i) printf("%lld\n", ans[i].res[2] - ans[i].res[1]);
}