解题：AHOI 2013 作业

题面

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

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

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

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