解题：AHOI 2013 作业

题面

emmm......我把莫队扔到了杂题里，因为感觉局限挺大的=。=

这题是莫队维护信息+分块查询答案，都是两者的基本操作，复杂度$O(m$ $sqrt(n)+n$ $sqrt(m))$

所以为啥要写这水题的题解来着

 #include<cmath>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 const int N=,Sq=;
 struct a
 {
     long long ans,num;
     int l,r,xx,yy,v,id;
 }mo[N];
 int b[N],blo[N],cnt[N],exi[Sq],bkt[Sq],pts[Sq][];
 int n,m,t1,t2,t3,t4,lp,rp,sqr,srt,xnt,maxx;
 bool cmp(a x,a y)
 {
     return x.v==y.v?x.r<y.r:x.v<y.v;
 }
 bool com(a x,a y)
 {
     return x.id<y.id;
 }
 void change(int val,int typ)
 {
     if(typ)
     {
         bkt[blo[val]]++;
         exi[blo[val]]+=(++cnt[val]==);
     }
     else
     {
         bkt[blo[val]]--;
         exi[blo[val]]-=(!(--cnt[val]));
     }
 }
 int query1(int x,int y)
 {
     int ret=;
     if(blo[x]!=blo[y])
     {
         for(int i=x;i<=pts[blo[x]][];i++) ret+=cnt[i];
         for(int i=pts[blo[y]][];i<=y;i++) ret+=cnt[i];
         for(int i=blo[x]+;i<=blo[y]-;i++) ret+=bkt[i];
     }
     else for(int i=x;i<=y;i++) ret+=cnt[i];
     return ret;
 }
 int query2(int x,int y)
 {
     int ret=;
     if(blo[x]!=blo[y])
     {
         for(int i=x;i<=pts[blo[x]][];i++) ret+=(cnt[i]>);
         for(int i=pts[blo[y]][];i<=y;i++) ret+=(cnt[i]>);
         for(int i=blo[x]+;i<=blo[y]-;i++) ret+=exi[i];
     }
     else for(int i=x;i<=y;i++) ret+=(cnt[i]>);
     return ret;
 }
 int main ()
 {
     scanf("%d%d",&n,&m),sqr=sqrt(n);
     for(int i=;i<=n;i++)
         scanf("%d",&b[i]),maxx=max(maxx,b[i]);
     for(int i=;i<=m;i++)
     {
         scanf("%d%d%d%d",&t1,&t2,&t3,&t4);
         mo[i].v=(t1-)/sqr+,mo[i].id=i,maxx=max(maxx,t4);
         mo[i].l=t1,mo[i].r=t2,mo[i].xx=t3,mo[i].yy=t4;
     }
     pts[xnt=][]=,srt=sqrt(maxx);
     for(int i=;i<=n;i++)
     {
         blo[i]=(i-)/srt+;
         if(i%srt==)
         {
             pts[xnt++][]=i;
             pts[xnt][]=i+;
         }
     }
     pts[xnt][]=maxx,lp=;
     sort(mo+,mo++m,cmp);
     for(int i=;i<=m;i++)
     {
         while(lp<mo[i].l) change(b[lp++],);
         while(lp>mo[i].l) change(b[--lp],);
         while(rp<mo[i].r) change(b[++rp],);
         while(rp>mo[i].r) change(b[rp--],);
         mo[i].ans=query1(mo[i].xx,mo[i].yy);
         mo[i].num=query2(mo[i].xx,mo[i].yy);
     }
     sort(mo+,mo++m,com);
     for(int i=;i<=m;i++)
         printf("%lld %lld\n",mo[i].ans,mo[i].num);
     return ;
 }