[HEOI2012]采花树状数组 BZOJ 2743

题目描述

萧薰儿是古国的公主，平时的一大爱好是采花。

今天天气晴朗，阳光明媚，公主清晨便去了皇宫中新建的花园采花。

花园足够大，容纳了n朵花，花有c种颜色（用整数1-c表示），且花是排成一排的，以便于公主采花。公主每次采花后会统计采到的花的颜色数，颜色数越多她会越高兴！同时，她有一癖好，她不允许最后自己采到的花中，某一颜色的花只有一朵。为此，公主每采一朵花，要么此前已采到此颜色的花，要么有相当正确的直觉告诉她，她必能再次采到此颜色的花。

由于时间关系，公主只能走过花园连续的一段进行采花，便让女仆福涵洁安排行程。福涵洁综合各种因素拟定了m个行程，然后一一向你询问公主能采到多少朵花（她知道你是编程高手，定能快速给出答案！），最后会选择令公主最高兴的行程（为了拿到更多奖金！）。

输入输出格式

输入格式：

第一行四个空格隔开的整数n、c以及m。接下来一行n个空格隔开的整数，每个数在[1, c]间，第i个数表示第i朵花的颜色。接下来m行每行两个空格隔开的整数l和r（l ≤ r），表示女仆安排的行程为公主经过第l到第r朵花进行采花。

输出格式：

共m行，每行一个整数，第i个数表示公主在女仆的第i个行程中能采到的花的颜色数。

输入输出样例

输入样例#1：
复制

输出样例#1： 复制

说明

对于100%的数据，1 ≤ n ≤ 2∗1062*10^62∗106，c ≤ n，m ≤2∗1062*10^62∗106。

本题有两个subtask

subtask1保证n,m,c≤3∗105n,m,c \leq 3*10^5n,m,c≤3∗105，占100分

subtask2保证n,m,c≤2∗106n,m,c \leq 2*10^6n,m,c≤2∗106，占100分

与HH的项链相似，同样是用树状数组；

用nxt表示下一个的位置，nnxt表示下一个的下一个位置；

#include<iostream>

#include<cstdio>

#include<algorithm>

#include<cstdlib>

#include<cstring>

#include<string>

#include<cmath>

#include<map>

#include<set>

#include<vector>

#include<queue>

#include<bitset>

#include<ctime>

#include<deque>

#include<stack>

#include<functional>

#include<sstream>

//#include<cctype>

//#pragma GCC optimize(2)

using namespace std;

#define maxn 2000005

#define inf 0x7fffffff

//#define INF 1e18

#define rdint(x) scanf("%d",&x)

#define rdllt(x) scanf("%lld",&x)

#define rdult(x) scanf("%lu",&x)

#define rdlf(x) scanf("%lf",&x)

#define rdstr(x) scanf("%s",x)

typedef long long  ll;

typedef unsigned long long ull;

typedef unsigned int U;

#define ms(x) memset((x),0,sizeof(x))

const long long int mod = 1e6 + 7;

#define Mod 1000000000

#define sq(x) (x)*(x)

#define eps 1e-4

typedef pair<int, int> pii;

#define pi acos(-1.0)

//const int N = 1005;

#define REP(i,n) for(int i=0;i<(n);i++)

typedef pair<int, int> pii;

inline int rd() {

	int x = 0;

	char c = getchar();

	bool f = false;

	while (!isdigit(c)) {

		if (c == '-') f = true;

		c = getchar();

	}

	while (isdigit(c)) {

		x = (x << 1) + (x << 3) + (c ^ 48);

		c = getchar();

	}

	return f ? -x : x;

}

ll gcd(ll a, ll b) {

	return b == 0 ? a : gcd(b, a%b);

}

int sqr(int x) { return x * x; }

/*ll ans;

ll exgcd(ll a, ll b, ll &x, ll &y) {

	if (!b) {

		x = 1; y = 0; return a;

	}

	ans = exgcd(b, a%b, x, y);

	ll t = x; x = y; y = t - a / b * y;

	return ans;

}

*/

int col[maxn];

int n, k, m;

int c[maxn];

int b[maxn];

struct node {

	int l, r, id;

}q[maxn];

bool cmp(node a, node b) {

	if(a.l!=b.l)return a.l < b.l;

	return a.r < b.r;

}

void add(int x, int y) {

	while (x <= n) {

		c[x] += y; x += x & -x;

	}

}

int query(int x) {

	int sum = 0;

	while (x > 0) {

		sum += c[x]; x -= x & -x;

	}

	return sum;

}

int fir[maxn];

int nxt[maxn], nnxt[maxn];

int ct[maxn];

int ans[maxn];

int main() {

//	ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);

	rdint(n); rdint(k); rdint(m);

	for (int i = 1; i <= n; i++)rdint(col[i]);

	for (int i = 1; i <= m; i++) {

		rdint(q[i].l); rdint(q[i].r); q[i].id = i;

	}

	for (int i = n; i >= 1; i--) {

		nxt[i] = fir[col[i]]; fir[col[i]] = i;

	}

	for (int i = 1; i <= n; i++)nnxt[i] = nxt[nxt[i]];

	for (int i = 1; i <= n; i++)if ((++ct[col[i]]) == 2)add(i, 1);

	sort(q + 1, q + 1 + m, cmp);

	int cur = 1;

	for (int i = 1; i <= m; i++) {

		for (; cur < q[i].l; cur++) {

			if (nxt[cur])add(nxt[cur], -1);

			if (nnxt[cur])add(nnxt[cur], 1);

		}

		ans[q[i].id] = query(q[i].r) - query(q[i].l - 1);

	}

	for (int i = 1; i <= m; i++)printf("%d\n", ans[i]);

	return 0;

}