[AHOI2013]作业

这个题一看就首选莫队，然而询问的东东·有些坑·，如果裸的莫队要跑140多s（极限数据，我试过，嘿嘿），所以当然要优化了，常规的优化方法是套一个树状数组来维护ans，

不过我要说的不是树状数组而是分块（分块是个好东东，秒杀树状数组，线段树。。。），对数值再分一次块。只要搞一个类似延迟标记的东东，f[],g[],f存每个块中数值个数，

g存每个块中数值种类个数。散点暴力，整块直接加上f，g，然后就可以A掉了。

代码如下：

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
const int N=100000+10;
const int M=1000000+10;
const int qrt=1000;
int n,m,bl,qr;
int block[N],beg[qrt+10],bend[qrt+10],val[N],f[qrt+10],g[qrt+10],cnt[N];
struct DATE{
	int l,r,mn,mx,id;
}date[M];
struct P{
	int o,t;
}ans[M];
inline int read(){
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}
inline bool cmp(DATE a,DATE b){
	if(block[a.l]==block[b.l]) return a.r<b.r;
	return block[a.l]<block[b.l];
}
inline void change(int pos,int va){
	int col=val[pos];int bll=block[col];
	cnt[col]+=va;f[bll]+=va;
	if((cnt[col]==1&&va>0)||(!cnt[col])) g[bll]+=va;
}
inline P query(int l,int r){
	int x=0,y=0;
	int le=block[l],ri=block[r];
	if(le==ri){
		for(int i=l;i<=r;i++){
			if(cnt[i]){
				x+=cnt[i];
				y+=1;
			}
		}
		P temp=(P){x,y};
		return temp;
	}
	for(int i=le+1;i<=ri-1;i++){
		x+=f[i];
		y+=g[i];
	}
	for(int i=l;i<=bend[le];i++){
		if(cnt[i]){
				x+=cnt[i];
				y+=1;
			}
	}
	for(int i=beg[ri];i<=r;i++){
		if(cnt[i]){
				x+=cnt[i];
				y+=1;
			}
	}
	P temp=(P){x,y};
	return temp;
}
inline void init(){
	for(int i=1;i<=bl;i++){
		beg[i]=qr*(i-1)+1;
		bend[i]=min(n,i*qr);
	}
}
int main(){
	freopen("ahoi2013_homework.in","r",stdin);
    	freopen("ahoi2013_homework.out","w",stdout);
	n=read();m=read();
	qr=sqrt(n+0.5);
	for(int i=1;i<=n;i++){
		val[i]=read();
		block[i]=(i-1)/qr+1;
	}
	bl=block[n];init();
	for(int i=1;i<=m;i++){
		date[i].l=read();date[i].r=read();
		date[i].mn=read();date[i].mx=read();
		date[i].id=i;
	}
	sort(date+1,date+m+1,cmp);
	int l=1,r=0;
	for(int i=1;i<=m;i++){
		for(;r<date[i].r;r++) change(r+1,1);
		for(;r>date[i].r;r--) change(r,-1);
		for(;l<date[i].l;l++) change(l,-1);
		for(;l>date[i].l;l--) change(l-1,1);
		int bh=date[i].id;
		ans[bh]=query(date[i].mn,date[i].mx);
	}
	for(int i=1;i<=m;i++)  printf("%d %d\n",ans[i].o,ans[i].t);
	return 0;
}