这个题一看就首选莫队,然而询问的东东·有些坑·,如果裸的莫队要跑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;
}