bzoj 1878: [SDOI2009]HH的项链 &&luogu P1972

1878: [SDOI2009]HH的项链

Time Limit: 4 Sec   Memory Limit: 64 MB
Submit: 4871   Solved: 2410
[ Submit][ Status][ Discuss]

Description

HH有一串由各种漂亮的贝壳组成的项链。HH相信不同的贝壳会带来好运，所以每次散步 完后，他都会随意取出一

段贝壳，思考它们所表达的含义。HH不断地收集新的贝壳，因此他的项链变得越来越长。有一天，他突然提出了一

个问题：某一段贝壳中，包含了多少种不同的贝壳？这个问题很难回答。。。因为项链实在是太长了。于是，他只

好求助睿智的你，来解决这个问题。

Input

第一行：一个整数N，表示项链的长度。 

第二行：N个整数，表示依次表示项链中贝壳的编号（编号为0到1000000之间的整数）。 

第三行：一个整数M，表示HH询问的个数。 

接下来M行：每行两个整数，L和R（1 ≤ L ≤ R ≤ N），表示询问的区间。

N ≤ 50000，M ≤ 200000。

Output

M行，每行一个整数，依次表示询问对应的答案。

Sample Input

6
1 2 3 4 3 5
3
1 2
3 5
2 6

Sample Output

2
2
4

HINT

Source

Day2

法一：（nlogn）

实测884ms

pre[i]表示i左边离i最近的j

使得a[i]==a[]j

（其实pre数组的作用就是让我们可以边扫边添加，并且使添加不重复）

将询问的右端点挂在序列上

扫一遍序列

每次将[pre[i]+1,i]这段区间+1

然后查询当前右端点下

询问左端点位置的值

#include<cstdio>
#include<cstring>
inline int read()
{
	int ans=0,f=1;char t=getchar();
	while(t<'0'||t>'9')	f=(t=='-'?-1:1),t=getchar();
	while(t>='0'&&t<='9')	ans=ans*10+t-'0',t=getchar();
	return ans*f;
}
inline int lowbit(int x)
{
	return x&(-x);
}
int n;
const int N=50007,M=2e5+7;
int c[N];
inline void modify(int x,int num)
{
	for( ;x<=n;x+=lowbit(x))	
	c[x]+=num;	
}
inline int query(int x)
{
	int ans=0;
	for( ;x;x-=lowbit(x))
	ans+=c[x];
	return ans;
}
int pos[1000007],pre[N];
struct node
{
	int l,next;
}e[M];
int first[N];int ans[M];int cnt;
inline void insert(int l,int r)
{
	e[++cnt]=(node){l,first[r]};first[r]=cnt;
}
int main()
{
	n=read();int u,v;
	for(int i=1;i<=n;i++)	
	{
		u=read();
		pre[i]=pos[u];
		pos[u]=i;
	}
	int m=read();
	for(int i=1;i<=m;i++)	u=read(),v=read(),insert(u,v);
	
	for(int i=1;i<=n;i++)
	{
		modify(pre[i]+1,1);	modify(i+1,-1);
		for(int k=first[i];k;k=e[k].next)	ans[k]=query(e[k].l);
	}
	for(int i=1;i<=m;i++)	printf("%d\n",ans[i]);
	return 0;
}

法二：直接莫队（n*sqrt(n) ）

实测1632ms

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
int read()
{
    int ans=0,f=1;char t=getchar();
    while(t<'0'||t>'9')   f=(t=='-'?-1:1),t=getchar();
    while(t>='0'&&t<='9') ans=ans*10+t-'0',t=getchar();
    return ans*f;
}
const int N=2e5+7,M=1e6+7;
int a[N],pos[N];
struct node
{
    int l,r,order;
}e[N];
inline bool cmp(node a,node b)
{
    if(pos[a.l]==pos[b.l])  return a.r<b.r;
    return pos[a.l]<pos[b.l];
}
int m;
int bo[M],final[N];
inline void work()
{
    int ans=0;
    for(int l=1,r=0,i=1;i<=m;i++)
    {
        while(r<e[i].r)  
        {
            r++;
            if(!bo[a[r]])   ans++;
            bo[a[r]]++;
        }
        while(r>e[i].r)
        {
            bo[a[r]]--;
            if(!bo[a[r]])   ans--;
            r--;
        }
        while(l<e[i].l)
        {
            bo[a[l]]--;
            if(!bo[a[l]])   ans--;
            l++;
        }
        while(l>e[i].l)  
        {
            l--;
            if(!bo[a[l]])   ans++;
            bo[a[l]]++;
        }
        final[e[i].order]=ans;
    }
}
int main()
{   
    int n=read();
    for(int i=1;i<=n;i++)    a[i]=read();
    m=read();
    for(int i=1;i<=m;i++)    e[i].l=read(),e[i].r=read(),e[i].order=i;
    int block=sqrt(n);
    for(int i=1;i<=n;i++)    pos[i]=(i-1)/block+1;
    std::sort(e+1,e+1+m,cmp);
    work();
    for(int i=1;i<=m;i++)    printf("%d\n",final[i]);
    return 0;
}