莫队算法

【前言】

莫队算法（Mo’s algorithm）是由莫涛队长发明的，一种处理区间问题的离线算法，由于其代码简便，常数巨小，适合各种打暴力。

另：关于带修改的莫队，可以看这里

【做法】

前提：如果已知区间[L,R]的答案，可以很快（log或常数级别）地得到区间[L±1,R±1]的答案。
对于区间询问[L,R]，我们可以先读进来所有询问，按照某种顺序依次处理（线性移动L，R端点），就可以得到所有询问的答案。

那么，处理各个询问的顺序就十分重要了。
我们可以利用分块的思想，把1~n分成√n块，对于L在同一块的询问，按R排序。
下面证明复杂度：
对于L，最坏是每次走√n，即O(Q√n)
对于R，有√n块，每个块最坏走n次，即O(n√n)
所以总复杂度是O(n√n)

【例题】

BZOJ 2038

#include<cstdio>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const int maxn=50005;
int n,q,a[maxn],h[maxn],hsh[maxn];
LL now;
LL c2(int x) {return (LL)x*(x-1)/2;}
LL gcd(LL x,LL y){
    if (y==0) return x;
    return gcd(y,x%y);
}
struct data{
    int l,r,id;
    bool operator<(const data&b)const{
        if (h[l]==h[b.l]) return r<b.r;
        return l<b.l;
    }
}que[maxn];
struct la{
    LL x,y;
}ans[maxn];
inline int red(){
    int tot=0,f=1;char ch=getchar();
    while (ch<'0'||'9'<ch) {if (ch=='-') f=-f;ch=getchar();}
    while ('0'<=ch&&ch<='9') tot=tot*10+ch-48,ch=getchar();
    return f*tot;
}
void blocker(){
    int k=sqrt(n);
    for (int i=1;i<=n;i++) h[i]=(i-1)/k+1;
}
void update(int x,int d){
    now-=c2(hsh[a[x]]);
    hsh[a[x]]+=d;
    now+=c2(hsh[a[x]]);
}
int main(){
    n=red(),q=red();
    for (int i=1;i<=n;i++) a[i]=red();
    blocker();
    for (int i=1;i<=q;i++) que[i].l=red(),que[i].r=red(),que[i].id=i;
    sort(que+1,que+1+q);
    update(1,1);
    for (int i=1,L=1,R=1;i<=q;i++){
        while (L<que[i].l) update(L++,-1);
        while (L>que[i].l) update(--L,1);
        while (R<que[i].r) update(++R,1);
        while (R>que[i].r) update(R--,-1);
        ans[que[i].id]=(la){now,c2(que[i].r-que[i].l+1)};
    }
    for (int i=1;i<=q;i++){
        LL t=gcd(ans[i].x,ans[i].y);
        ans[i].x/=t;ans[i].y/=t;
        if (ans[i].x==0) printf("0/1\n");else printf("%lld/%lld\n",ans[i].x,ans[i].y);
    }
    return 0;
}