1. 程式人生 > >BZOJ3809: Gty的二逼妹子序列 BZOJ3236: [Ahoi2013]作業

BZOJ3809: Gty的二逼妹子序列 BZOJ3236: [Ahoi2013]作業

BZOJ3809: Gty的二逼妹子序列

又是一道許可權題。。。

本蒟蒻沒錢氪金。。。

附上洛谷題面:

洛谷P4867 Gty的二逼妹子序列

題目描述

Autumn和Bakser又在研究Gty的妹子序列了!但他們遇到了一個難題。

對於一段妹子們,他們想讓你幫忙求出這之內美麗度\in [a,b][a,b]的妹子的美麗度的種類數。

為了方便,我們規定妹子們的美麗度全都在[1,n][1,n]中。
給定一個長度為$n(1 \le n \le 100000)$的正整數序列$s(1 \le si \le n)$,對於$m(1 \le m \le 1000000)$次詢問l,r,a,b,每次輸出$s_l \cdots s_r$中,權值$\in [a,b]$的權值的種類數。

輸入輸出格式

輸入格式:

 

第一行包括兩個整數$n,m(1 \le n \le 100000,1 \le m \le 1000000)$,表示數列ss中的元素數和詢問數。

第二行包括nn個整數$s1…sn(1 \le si \le n)$

接下來mm行,每行包括44個整數$l,r,a,b(1 \le l \le r \le n,1 \le a \le b \le n)$,意義見題目描述。

保證涉及的所有數在C++的int內。保證輸入合法。

 

輸出格式:

 

對每個詢問,單獨輸出一行,表示$s_l \cdots s_r$中權值$\in [a,b]$的權值的種類數。

 

輸入輸出樣例

輸入樣例#1:  複製
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
輸出樣例#1:  複製
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。

建議使用輸入/輸出優化。


題解Here!

怎麼跟這題好像:

BZOJ3236: [Ahoi2013]作業

一看,只是把第一個詢問刪掉了。。。

所以還是莫隊+樹狀陣列。

正解當然是莫隊+分塊,但是我不想寫啊。。。

然後這個莫隊過百萬系列是個什麼操作???

附程式碼:

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#define MAXN 100010
#define MAXM 1000010
using namespace std;
int n,m,q,block;
int val[MAXN],num[MAXN],ans[MAXM];
struct Question{
	int l,r,a,b,id;
	friend bool operator <(const Question &p,const Question &q){
		return (p.r/block==q.r/block?(((p.r/block)&1)?p.l>q.l:p.l<q.l):p.r<q.r);
	}
}que[MAXM];
inline int read(){
	int date=0,w=1;char c=0;
	while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}
	while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();}
	return date*w;
}
namespace BIT{
	int bit[MAXN];
	inline int lowbit(int x){return x&(-x);}
	inline void add(int x,int v){for(;x<=n;x+=lowbit(x))bit[x]+=v;}
	inline int sum(int x){int s=0;for(;x;x-=lowbit(x))s+=bit[x];return s;}
}
inline void add(int x){
	if(!num[x])BIT::add(x,1);
	num[x]++;
}
inline void del(int x){
	num[x]--;
	if(!num[x])BIT::add(x,-1);
}
void work(){
	int left=1,right=0;
	for(int i=1;i<=m;i++){
		while(left<que[i].l)del(val[left++]);
		while(left>que[i].l)add(val[--left]);
		while(right<que[i].r)add(val[++right]);
		while(right>que[i].r)del(val[right--]);
		ans[que[i].id]=BIT::sum(que[i].b)-BIT::sum(que[i].a-1);
	}
	for(int i=1;i<=m;i++)printf("%d\n",ans[i]);
}
void init(){
	n=read();m=read();
	for(int i=1;i<=n;i++)val[i]=read();
	block=sqrt(n);
	for(int i=1;i<=m;i++){
		que[i].l=read();que[i].r=read();que[i].a=read();que[i].b=read();
		que[i].id=i;
	}
	sort(que+1,que+m+1);
}
int main(){
	init();
	work();
    return 0;
}