【BZOJ1935】【SHOI2007】园丁的烦恼（树状数组，CDQ分治？）

Description

很久很久以前，在遥远的大陆上有一个美丽的国家。统治着这个美丽国家的国王是一个园艺爱好者，在他的皇家花园里种植着各种奇花异草。

有一天国王漫步在花园里，若有所思，他问一个园丁道： “最近我在思索一个问题，如果我们把花坛摆成六个六角形，那么……”

“那么本质上它是一个深度优先搜索，陛下”，园丁深深地向国王鞠了一躬。

“嗯……我听说有一种怪物叫九头蛇，它非常贪吃苹果树……”

“是的，显然这是一道经典的动态规划题，早在N元4002年我们就已经发现了其中的奥秘了，陛下”。

“该死的，你究竟是什么来头？”

“陛下息怒，干我们的这行经常莫名其妙地被问到和OI有关的题目，我也是为了预防万一啊！” 王者的尊严受到了伤害，这是不可容忍的。

看来一般的难题是难不倒这位园丁的，国王最后打算用车轮战来消耗他的实力： “年轻人，在我的花园里的每一棵树可以用一个整数坐标来表示，一会儿，我的骑士们会来轮番询问你某一个矩阵内有多少树，如果你不能立即答对，你就准备走人吧！”说完，国王气呼呼地先走了。

这下轮到园丁傻眼了，他没有准备过这样的问题。所幸的是，作为“全国园丁保护联盟”的会长——你，可以成为他的最后一根救命稻草。

Solution

其实只是一个二维偏序而已。。。

先用二维前缀和将每个询问转化一下，然后离线操作。
先按$x$排序，然后将$y$依次插入树状数组，每次查询即可。

Source

/************************************************
 * Au: Hany01
 * Date: Mar 13th, 2018
 * Prob: [BZOJ1935][SHOI2007] 园丁的烦恼
 * Email: hany01@foxmail.com
************************************************/

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;
typedef pair<int, int> PII;
#define rep(i, j) for (register int i = 0, i##_end_ = (j); i < i##_end_; ++ i)
#define For(i, j, k) for (register int i = (j), i##_end_ = (k); i <= i##_end_; ++ i)
#define Fordown(i, j, k) for (register int i = (j), i##_end_ = (k); i >= i##_end_; -- i)
#define Set(a, b) memset(a, b, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define fir first
#define sec second
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define ALL(a) (a).begin(), (a).end()
#define SZ(a) ((int)(a).size())
#define INF (0x3f3f3f3f)
#define INF1 (2139062143)
#define Mod (1000000007)
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define y1 wozenmezhemecaia

template <typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, 1 : 0; }
template <typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, 1 : 0; }

inline int read()
{
    register int _, __; register char c_;
    for (_ = 0, __ = 1, c_ = getchar(); c_ < '0' || c_ > '9'; c_ = getchar()) if (c_ == '-') __ = -1;
    for ( ; c_ >= '0' && c_ <= '9'; c_ = getchar()) _ = (_ << 1) + (_ << 3) + (c_ ^ 48);
    return _ * __;
}

inline void File()
{
#ifdef hany01
    freopen("bzoj1935.in", "r", stdin);
    freopen("bzoj1935.out", "w", stdout);
#endif
}

struct Data {
    int x, y, id;
    Data(int x = 0, int y = 0, int id = 0): x(x), y(y), id(id) {}
};
bool operator < (const Data& A, const Data& B) {
    return A.x < B.x ? 1 : (A.x == B.x && A.y < B.y);
}
struct Question { int x1, y1, x2, y2; };

const int maxn = 500005;

int n, m, cnt, totx, toty, lsx[maxn << 1], lsy[maxn << 1], Ans[maxn << 2];
Question Q[maxn];
Data a[maxn], q[maxn << 2];

inline void Init()
{
    n = read(), m = read();
    For(i, 1, n) a[i].x = read() + 1, a[i].y = read() + 1;
    sort(a + 1, a + 1 + n);

    For(i, 1, m)
        Q[i].x1 = read() + 1, Q[i].y1 = read() + 1, Q[i].x2 = read() + 1, Q[i].y2 = read() + 1,
        q[++ cnt] = Data(Q[i].x1 - 1, Q[i].y1 - 1, i), q[++ cnt] = Data(Q[i].x2, Q[i].y2, i),
        q[++ cnt] = Data(Q[i].x1 - 1, Q[i].y2, -i), q[++ cnt] = Data(Q[i].x2, Q[i].y1 - 1, -i),
        lsx[++ totx] = Q[i].x1 - 1, lsx[++ totx] = Q[i].x2, lsy[++ toty] = Q[i].y1 - 1, lsy[++ toty] = Q[i].y2;
    sort(q + 1, q + 1 + cnt);

    sort(lsx + 1, lsx + 1 + totx), sort(lsy + 1, lsy + 1 + toty);
    totx = unique(lsx + 1, lsx + 1 + totx) - lsx - 1, toty = unique(lsy + 1, lsy + 1 + toty) - lsy - 1;
    For(i, 1, cnt) q[i].x = lower_bound(lsx + 1, lsx + 1 + totx, q[i].x) - lsx, q[i].y = lower_bound(lsy + 1, lsy + 1 + toty, q[i].y) - lsy;
}

struct BIT {
    int c[maxn << 1];
#define lb(x) ((x) & -(x))
    inline void update(int x) { for ( ; x <= n; x += lb(x)) ++ c[x]; }
    inline int query(int x) { register int Ans = 0; for ( ; x; x -= lb(x)) Ans += c[x]; return Ans; }
}bit;

inline void Solve()
{
    int cur1 = 1, cur2 = 1;
    For(i, 1, totx) {
        while (cur1 <= n && a[cur1].x <= lsx[i]) bit.update(lower_bound(lsy + 1, lsy + 1 + toty, a[cur1 ++].y) - lsy);
        while (cur2 <= cnt && q[cur2].x <= i) Ans[cur2] = bit.query(q[cur2].y), ++ cur2;
    }
}

inline void Answer()
{
    static int Prt[maxn];
    For(i, 1, cnt) Prt[abs(q[i].id)] += (q[i].id > 0 ? 1 : -1) * Ans[i];
    For(i, 1, m) printf("%d\n", Prt[i]);
}

int main()
{

    File();

    Init();

    Solve();

    Answer();

    return 0;

}
//乱山残雪夜，孤烛异乡人。
//    -- 崔涂《巴山道中除夜书怀》