【bzoj3809】Gty的二逼妹子序列

*题目描述：

Autumn和Bakser又在研究Gty的妹子序列了！但他们遇到了一个难题。

对于一段妹子们，他们想让你帮忙求出这之内美丽度∈[a,b]的妹子的美丽度的种类数。

为了方便，我们规定妹子们的美丽度全都在[1,n]中。

给定一个长度为n(1<=n<=100000)的正整数序列s(1<=si<=n)，对于m(1<=m<=1000000)次询问“l,r,a,b”，每次输出sl…sr中，权值∈[a,b]的权值的种类数。

*输入：

第一行包括两个整数n,m(1<=n<=100000,1<=m<=1000000),表示数列s中的元素数和询问数。

第二行包括n个整数s1…sn(1<=si<=n)。

接下来m行,每行包括4个整数l,r,a,b(1<=l<=r<=n,1<=a<=b<=n),意义见题目描述。

保证涉及的所有数在C++的int内。

保证输入合法。

*输出：

对每个询问，单独输出一行，表示sl…sr中权值∈[a,b]的权值的种类数。

样例输入：

10 10

4 4 5 1 4 1 5 1 2 1

5 9 1 2

3 4 7 9

4 4 2 5

2 3 4 7

5 10 4 4

3 9 1 1

1 4 5 9

8 9 3 3

2 2 1 6

8 9 1 4

*样例输出：

2

0

0

2

1

1

1

0

1

2

*提示：

样例的部分解释：

5 9 1 2

子序列为4 1 5 1 2

在[1,2]里的权值有1,1,2，有2种，因此答案为2。

3 4 7 9

子序列为5 1

在[7,9]里的权值有5，有1种，因此答案为1。

4 4 2 5

子序列为1

没有权值在[2,5]中的，因此答案为0。

2 3 4 7

子序列为4 5

权值在[4,7]中的有4,5，因此答案为2。

建议使用输入/输出优化。

*题解：

莫队+分块。将询问离线莫队分块以后，再按a和b的权值分块来统计区间的个数。时间复杂度为O(n*sqrt(n))。

*代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>

#ifdef WIN32
#define LL "%I64d"
#else
#define LL "%lld"
#endif

#ifdef CT
#define debug(...) printf(__VA_ARGS__)
#define setfile()
#else
#define debug(...)
#define filename ""
#define setfile() freopen(filename".in", "r", stdin); freopen(filename".out", "w", stdout);
#endif

#define R register
#define getc() (S == T && (T = (S = B) + fread(B, 1, 1 << 15, stdin), S == T) ? EOF : *S++)
#define dmax(_a, _b) ((_a) > (_b) ? (_a) : (_b))
#define dmin(_a, _b) ((_a) < (_b) ? (_a) : (_b))
#define cmax(_a, _b) (_a < (_b) ? _a = (_b) : 0)
#define cmin(_a, _b) (_a > (_b) ? _a = (_b) : 0)
char B[1 << 15], *S = B, *T = B;
inline int FastIn()
{
R char ch; R int cnt = 0; R bool minus = 0;
while (ch = getc(), (ch < '0' || ch > '9') && ch != '-') ;
ch == '-' ? minus = 1 : cnt = ch - '0';
while (ch = getc(), ch >= '0' && ch <= '9') cnt = cnt * 10 + ch - '0';
return minus ? -cnt : cnt;
}
#define maxn 100010
#define maxm 1000010
int a[maxn], block, n, id[maxn], cnt[maxn], bcnt[maxn], ans[maxm];
struct Query
{
int l, r, a, b, id;
}q[maxm];
inline bool operator < (const Query &i, const Query &j)
{
return id[i.l] != id[j.l] ? i.l < j.l : (id[i.l] & 1 ? i.r > j.r : i.r < j.r);
}
inline int query(R int aa, R int bb)
{
R int tmp = 0;
if (id[aa] == id[bb])
{
for (R int i = aa; i <= bb; ++i)
if (cnt[i]) ++tmp;
}
else
{
for (R int i = aa; i < (id[aa] + 1) * block; ++i)
if (cnt[i]) ++tmp;
for (R int i = id[bb] * block; i <= bb; ++i)
if (cnt[i]) ++tmp;
for (R int i = id[aa] + 1; i < id[bb]; ++i)
tmp += bcnt[i];
}
return tmp;
}
inline void add(R int x)
{
++cnt[x];
if (cnt[x] == 1) ++bcnt[id[x]];
}
inline void del(R int x)
{
--cnt[x];
if (!cnt[x]) --bcnt[id[x]];
}
int main()
{
//  setfile();
n = FastIn(); R int m = FastIn(); block = sqrt(n);
for (R int i = 1; i <= n; ++i) a[i] = FastIn(), id[i] = i / block;
for (R int i = 1; i <= m; ++i)
{
q[i] = (Query) {FastIn(), FastIn(), FastIn(), FastIn(), i};
}
std::sort(q + 1, q + m + 1);
R int l = 1, r = 0;
for (R int i = 1; i <= m; ++i)
{
while (l < q[i].l) del(a[l]), ++l;
while (r > q[i].r) del(a[r]), --r;
while (l > q[i].l) --l, add(a[l]);
while (r < q[i].r) ++r, add(a[r]);
ans[q[i].id] = query(q[i].a, q[i].b);
}
for (R int i = 1; i <= m; ++i) printf("%d\n", ans[i] );
return 0;
}