Description
很久很久以前,在遥远的大陆上有一个美丽的国家。统治着这个美丽国家的国王是一个园艺爱好者,在他的皇家花园里种植着各种奇花异草。
有一天国王漫步在花园里,若有所思,他问一个园丁道: “最近我在思索一个问题,如果我们把花坛摆成六个六角形,那么……”
“那么本质上它是一个深度优先搜索,陛下”,园丁深深地向国王鞠了一躬。
“嗯……我听说有一种怪物叫九头蛇,它非常贪吃苹果树……”
“是的,显然这是一道经典的动态规划题,早在N元4002年我们就已经发现了其中的奥秘了,陛下”。
“该死的,你究竟是什么来头?”
“陛下息怒,干我们的这行经常莫名其妙地被问到和OI有关的题目,我也是为了预防万一啊!” 王者的尊严受到了伤害,这是不可容忍的。
看来一般的难题是难不倒这位园丁的,国王最后打算用车轮战来消耗他的实力: “年轻人,在我的花园里的每一棵树可以用一个整数坐标来表示,一会儿,我的骑士们会来轮番询问你某一个矩阵内有多少树,如果你不能立即答对,你就准备走人吧!”说完,国王气呼呼地先走了。
这下轮到园丁傻眼了,他没有准备过这样的问题。所幸的是,作为“全国园丁保护联盟”的会长——你,可以成为他的最后一根救命稻草。
Solution
其实只是一个二维偏序而已。。。
先用二维前缀和将每个询问转化一下,然后离线操作。
先按
x
x
排序,然后将依次插入树状数组,每次查询即可。
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;
}
//乱山残雪夜,孤烛异乡人。
// -- 崔涂《巴山道中除夜书怀》