[BZOJ][HEOI2012]采花 离线

Description

萧芸斓是Z国的公主,平时的一大爱好是采花。
今天天气晴朗,阳光明媚,公主清晨便去了皇宫中新建的花园采花。花园足够大,容纳了n朵花,花有c种颜色(用整数1-c表示),且花是排成一排的,以便于公主采花。公主每次采花后会统计采到的花的颜色数,颜色数越多她会越高兴!同时,她有一癖好,她不允许最后自己采到的花中,某一颜色的花只有一朵。为此,公主每采一朵花,要么此前已采到此颜色的花,要么有相当正确的直觉告诉她,她必能再次采到此颜色的花。由于时间关系,公主只能走过花园连续的一段进行采花,便让女仆福涵洁安排行程。福涵洁综合各种因素拟定了m个行程,然后一一向你询问公主能采到多少朵花(她知道你是编程高手,定能快速给出答案!),最后会选择令公主最高兴的行程(为了拿到更多奖金!)。

Input


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

Output

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

Sample Input

5 3 5
1 2 2 3 1
1 5
1 2
2 2
2 3
3 5

Sample Output

2
0 0 1 0
【样例说明】
询问[1, 5]:公主采颜色为1和2的花,由于颜色3的花只有一朵,公主不采;询问[1, 2]:颜色1和颜色2的花均只有一朵,公主不采;
询问[2, 2]:颜色2的花只有一朵,公主不采;
询问[2, 3]:由于颜色2的花有两朵,公主采颜色2的花;
询问[3, 5]:颜色1、2、3的花各一朵,公主不采。

HINT

【数据范围】

对于100%的数据,1 ≤ n ≤    10^6,c ≤ n,m ≤10^6。

Source

蹂躏,我数据结构学的真是毛线啊,这道题很明显是离线做的,看到题就像到了。
然后接着想。
然后继续想。
然后就没有然后了。
最后还是看了别人题解才做出来的(真不知道当时是怎么写得"生日礼物")。

先是对询问的区间按r[i]拍个序。
令上一次出现颜色为i的位置为last[i],上上次出现的位置为llast[i];
然后由1到n扫,每扫到一个点i,我们在llast[i]+1到last[i]上加一个一,然后回答询问r[j]为i的区间。
为什么这样就能成功呢?
我们这样观察,我们是由左向右扫的,而 题目要求的是至少两个颜色相同的点的个数,也就是说,对于任意一个颜色,他能算作被取的左端点是递增的,换句话说,因为现在和以后我们要回答的询问的右端点一直是往右走,所以一个颜色生效的最大右端点也是一直往右走的,所以我们每次维护一个上次的位置,上上次的位置,并在中间加一(不在1-上次的位置上加一是因为以前在1-上上次的位置上加过一了,所以每次只需维护一段距离最近的同种颜色的点就够了)

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<ctime>
using namespace std;
int n,c,m;
int col[2000000];
int lft[2000000],rgt[2000000],rank[2000000];
void kp(int l,int r){
	int i=l,j=r,mid=rgt[(l+r)/2];
	while (i<j){
		while (rgt[i]<mid) i++;
		while (mid<rgt[j]) j--;
		if (i<=j){
			swap(lft[i],lft[j]);
			swap(rgt[i],rgt[j]);
			swap(rank[i],rank[j]);
			i++;j--;
			}
		}
	if (l<j) kp(l,j);
	if (i<r) kp(i,r);
}
struct node{
	int l,r;
	int sum;
}	f[5000000];
void build(int k,int l,int r){
	f[k].l=l;f[k].r=r;
	if (l==r)	return;
	build(k*2,l,(l+r)/2);
	build(k*2+1,(l+r)/2+1,r);
}
void updata(int k){
	f[k].sum=f[k*2].sum+f[k*2+1].sum;
}
void modify(int k,int aim,int delta){
	if (f[k].l==aim&&f[k].r==aim){
		f[k].sum+=delta;
		return;
		}
	if (aim<=f[k*2].r) modify(k*2,aim,delta);
	if (f[k*2+1].l<=aim) modify(k*2+1,aim,delta);
	updata(k);
}
int query(int k,int l,int r){
	if (l<=f[k].l&&f[k].r<=r){
		return f[k].sum;
		}
	int cnt=0;
	if (l<=f[k*2].r) cnt+=query(k*2,l,r);
	if (f[k*2+1].l<=r) cnt+=query(k*2+1,l,r);
	return cnt;
}
int ans[2000000],last[2000000],llast[2000000];
int main(){
//	freopen("2743.in","r",stdin);
//	freopen("2743.out","w",stdout);
	scanf("%d%d%d",&n,&c,&m);
	build(1,0,n);
	for (int i=1;i<=n;i++) scanf("%d",&col[i]);
	for (int i=1;i<=m;i++){
		scanf("%d%d",&lft[i],&rgt[i]);
		rank[i]=i;
		}
	kp(1,m);
//	for (int i=1;i<=m;i++) printf("%d %d\n",lft[i],rgt[i]);
	
	int now=1;
	for (int i=1;i<=n;i++){
		modify(1,last[col[i]],1);
		modify(1,llast[col[i]],-1);
//		printf("<%d %d %d>\n",i,llast[col[i]]+1,last[col[i]]);
		llast[col[i]]=last[col[i]];
		last[col[i]]=i;
		while (rgt[now]==i){
			ans[rank[now]]=query(1,lft[now],rgt[now]);
			now++;
			}
		}
//	printf("<%d>\n",now);
	for (int i=1;i<=m;i++) printf("%d\n",ans[i]);
	return 0;
}


  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值