BZOJ---2038：小Z的袜子【莫队入门】

题目：

https://www.lydsy.com/JudgeOnline/problem.php?id=2038

分析：

对于这类区间统计问题，都会尽量利用上一次的统计来加速本次区间的统计，也就是说上一次和本次区间的重叠部分只统计一次

不难想到把区间按左端点排序，这样排序后，统计时左端点一直向右移动，但右端点可能会从头移到尾，从尾移到头，不够优秀

莫队算法：将1 - N分成 N/sqrt(M)个块，按块的位置来排序，复杂度证明参考：戳戳戳~~

代码：

#include <bits/stdc++.h>
 
using namespace std;
typedef long long LL;
#define pii pair<LL,LL>
const int MAXN = 5e4+54;
int n,m,block,a[MAXN],vis[MAXN];
LL num;
struct segment{
    int l,r,id; 
}s[MAXN];
pii ans[MAXN];
bool cmp(segment a,segment b){
    return a.r/block == b.r/block ? a.l < b.l : a.r < b.r;
}
inline void add(int x){
    num += vis[a[x]];
    vis[a[x]] ++;
}
inline void del(int x){
    vis[a[x]] --;
    num -= vis[a[x]];
}
int main(){
    cin >> n >> m;
    for(int i = 1;i <= n; ++i) cin >> a[i];
    for(int i = 1;i <= m; ++i) s[i].id = i,cin >> s[i].l >> s[i].r;
    block = sqrt(n*1.0);
    sort(s+1,s+m+1,cmp);
    int L = 1,R = 0;
    for(int i = 1;i <= m; ++i){
        while(L < s[i].l) del(L++);
        while(L > s[i].l) add(--L);
        while(R < s[i].r) add(++R);
        while(R > s[i].r) del(R--);
        LL sum = 1LL * (R-L+1) * (R-L) / 2;
        LL g = __gcd(num,sum);
        ans[s[i].id] = make_pair(num/g,sum/g);
    }
    for(int i = 1;i <= m; ++i) cout<<ans[i].first<<"/"<<ans[i].second<<'\n';
    return 0;
}