解题报告：D - 小Z的袜子(hose) HYSBZ - 2038 莫队算法

题目链接

学习博客来源：点击打开链接

这题为莫队亲自出的莫队算法的练习题，必做的题，写下来记录一下莫队算法以便自己以后进行查看。

莫队算法可以用来解决静态离线的区间查询问题。

莫队算法的关键在于 已知[L,R]区间的情况下，能否快速得到其相邻4个区间的答案。

莫队算法的高效在于对询问的离线分块排序函数，以块为第一关键词，y为第二关键词，块的长度直接决定莫队算法的效率，为了保持最快一般采用sqrt()，如果空间有较严格的要求可以考虑n^(2/3) 。

莫队算法也可以用来解决动态离线的区间查询问题，复杂度大约为n^(5/3)，详细过程见大佬的博客。

由于会大量用到更新操作，所以对更新操作的计算优化会大大提升代码效率。

代码：

#include<bits/stdc++.h>

using namespace std;

int A[50005];
int num[50005];
long long ans[50005][2];
int pos[50005];
long long up,down;
int len , n , q ;
struct event{
    int l,r,pos;
}Q[50005];
bool cmp(event a,event b){
    if(pos[a.l]!=pos[b.l])return pos[a.l]<pos[b.l];
    return a.r<b.r;
}
long long C(int x){
    if(x<=1)return 0;
    return 1LL*x*(x-1)/2;
}

inline void det(int x){
    up = up - C(num[x]) + C(--num[x]);
}
inline void add(int x){
    up = up - C(num[x]) + C(++num[x]);
}


int main()
{
    while(scanf("%d%d",&n,&q)==2){
        int k = sqrt(n+0.5);
        for(int i=1;i<=n;i++){
            scanf("%d",&A[i]);
            pos[i] = (i-1)/k;
        }for(int i=1;i<=q;i++){
            scanf("%d%d",&Q[i].l,&Q[i].r);
            Q[i].pos = i;
        }memset(num,0,sizeof(num));
        sort(Q+1,Q+1+q,cmp);
        int l = 1 ,r = 2;
        num[A[l]]++;num[A[r]]++;
        len = 2 , up = (A[l]==A[r]) , down = 1;
        for(int i=1;i<=q;i++){
            while(r<Q[i].r)add(A[++r]);
            while(l>Q[i].l)add(A[--l]);
            while(r>Q[i].r)det(A[r--]);
            while(l<Q[i].l)det(A[l++]);
            int k = Q[i].pos ;
            down = C(r-l+1);
            long long gcd = __gcd(up,down);
            ans[k][0] = up/gcd;
            ans[k][1] = down/gcd;
        }for(int i=1;i<=q;i++){
            printf("%lld/%lld\n",ans[i][0],ans[i][1]);
        }
    }return 0;
}