算法学习->莫队算法（小z的袜子）

一、什么是莫队算法

莫队算法是一个对于区间、树或其他结构离线（在线）维护的算法，此算法基于一些基本算法，例如暴力维护，树状数组，分块，最小曼哈顿距离生成树，对其进行揉合从而产生的一个简单易懂且短小好写的算法。此算法在很多情况下可以很轻松的切掉一些复杂而且难写的数据结构问题。莫队算法其实只是一种将暴力枚举优化的算法，主要是多次统计任意一段区间内数的和，只能用于离线。当知道【l，r】区间的解时，可推出【l+1，r】、【l，r+1】、【l-1，r】、【l，r-1】的解。

二、莫队算法适用问题

莫队算法在解决离线区间询问几乎是无敌的。
几乎只要能离线，用分块的莫队算法都能取得一个令人满意的的解法（当有修改操作时，需增加一个数组t【】记录修改的时间）。
所以就有很多扩展（解决线段树等数据结构由于需要区间加和性而不能解决的问题），如区间众数，平均数什么的。

三、莫队算法模板

int s[maxn],vis[maxn],Ans[maxn],md[maxn];

struct Query{
    int lef,rig,id;
}Que[maxn];

bool cmp(struct Query a,struct Query b)
{
    if(md[a.lef]!=md[b.lef])return a.lef<b.lef;//当有时间数组时，先按l大小排序，
    return a.rig<b.rig;//在l相等的情况下，按r大小排序，在r相等的情况下，再按t数组排序 
}

void modui()
{
    int i,lef=0,rig=q+1,ans=0;
    for(i=1;i<=m;i++)
    {   //不断跳动 
        while(rig<Que[i].rig)   //[l,r]->[l,r+1] 
        {
            vis[s[rig]]--;
            if(vis[s[rig]]==0)
                ans--;
            rig++;
        }
        while(rig>Que[i].rig)   //[l,r]->[l,r-1] 
        {
            rig--;
            if(vis[s[rig]]==0)
                ans++;
            vis[s[rig]]++;
        }
        while(lef<Que[i].lef)   //[l,r]->[l+1,r] 
        {
            lef++;
            if(vis[s[lef]]==0)
                ans++;
            vis[s[lef]]++;
        }
        while(lef>Que[i].lef)   //[l,r]->[l-1,r] 
        {
            vis[s[lef]]--;
            if(vis[s[lef]]==0)
                ans--;
            lef--;
        }
        Ans[Que[i].id]=ans;
    }
    return ;
}

四、典型莫队算法题

小z的袜子
AC代码：

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdio>
typedef long long ll;
using namespace std;

ll gcd(ll a,ll b){
    return b==0?a:gcd(b,a%b);
}

ll n,m,col[50005],md[50005],ans1[50005],ans2[50005],block,cnt[50005];
struct Query{
    ll l,r,id;
}query[50005];

bool cmp(struct Query x,struct Query y){    //m个查询先按照第几个块排序，再按照 R排序。(不然会tle)
    return x.l/block==y.l/block?x.r<y.r:x.l<y.l;
}

void modui(){
    ll R=0,L=1,res=0;
    memset(cnt,0,sizeof(cnt));
    for(int i=1;i<=m;i++){
        while(R<query[i].r){
            R++;
            res+=(cnt[col[R]]+1)*(cnt[col[R]]+1)-cnt[col[R]]*cnt[col[R]];
            cnt[col[R]]++;
        }
        while(L<query[i].l){
            res-=cnt[col[L]]*cnt[col[L]]-(cnt[col[L]]-1)*(cnt[col[L]]-1);
            cnt[col[L]]--;
            L++;
        }
        while(R>query[i].r){
            res-=cnt[col[R]]*cnt[col[R]]-(cnt[col[R]]-1)*(cnt[col[R]]-1);
            cnt[col[R]]--;
            R--;
        }
        while(L>query[i].l){
            L--;
            res+=(cnt[col[L]]+1)*(cnt[col[L]]+1)-cnt[col[L]]*cnt[col[L]];
            cnt[col[L]]++;
        }
        ans1[query[i].id]=res-(R-L+1);
        ans2[query[i].id]=1LL*(R-L+1)*(R-L);
    }
}

int main(){
    scanf("%lld%lld",&n,&m);
    block=sqrt(n);
    for(ll i=1;i<=n;i++){
        scanf("%lld",&col[i]);
        md[i]=((i-1)/block)+1;
    }
    for(ll i=1;i<=m;i++){
        scanf("%lld%lld",&query[i].l,&query[i].r);
        query[i].id=i;
    }
    sort(query+1,query+m+1,cmp);
    modui();
    for(int i=1;i<=m;i++){  //最简化 
        int G=gcd(ans1[i],ans2[i]);
        ans1[i]/=G;
        ans2[i]/=G;
        if(!ans1[i]) ans2[i]=1;
    }
    for(ll i=1;i<=m;i++){
        if(ans1[i]==0){
            printf("0/%lld\n",ans2[i]);
        }else{
            printf("%lld/%lld\n",ans1[i],ans2[i]);
        }
    }
    return 0;
}