【莫队算法】小z的袜子（入门）

传说中，能解决一切可离线问题的莫队算法。一直没时间去认真看，看别人博客又看不懂，代码又感觉特别长。
最近小学弟去撸了下，发现看懂了，也顺带把我教会了=w=。

首先得要先懂得分块。
分块：
就是像桶排序那样，把一类东西丢进同一个桶里。

莫队里就是把问题的左区间分块了。

接下来讲莫队算法。
莫队的时间复杂度最差的情况下应该是$q*\sqrt{n}$
这种复杂度在大多数的$10^5$的询问，都是够用的

我们可以想象，对于一个区间的询问，我们可以拿两个指针去维护询问的左右区间。
假设现在有l指针和r指针，分别代表现在所在的左区间和右区间。

当我们询问到L和R时
如果发现l指针在L的左边，那么我们应该将l经过的位置的数值去除，然后l向L靠拢
如果发现l指针在L的右边，那么我们应该将l经过的位置的数值增加，然后l向L靠拢
R的话同理。

那这种操作的复杂度。最差情况就是。那人一次问你[1,n]一次问你[n/2,n/2]
这样操作法每次操作的复杂度都是O(n)，那么询问次数为q的时候则为O(n * q)

这样的复杂度是不能满足我们需求的

所以莫队将分块和以上那种操作结合到了一起，于是成了莫队算法

我们假设按照左指针分块，分成sqrt(n)块，右指针从小到大排序，那么左指针每次移动的范围为sqrt(n)，
右指针移动的范围为剩下的范围。
这里就不继续证明莫队的时间复杂度了。可以去看一下国家集训队莫队的论文

通常要修改的地方就是add和del函数

flag记录的是那个位置的计数器。

void add(int x) {
    ans += flag[a[x]];  
    flag[a[x]]++;
}

void del(int x) {
    flag[a[x]]--;
    ans -= flag[a[x]];
}

排序询问的时候如果同块，则按照r来排，如果不同块，按照块的顺序来排

int cmp(Question a, Question b) {
    if (pos[a.l] != pos[b.l]) {
        return pos[a.l] < pos[b.l];
    }
    return a.r < b.r;
}

再来说说这道题。
送上传送门：HYSBZ 2038

题意：
n只有序的袜子，每次询问一个区间，问随意拿一对袜子，能选到一样袜子颜色的概率，要求求最简分数。

思路：
我们对所有查询离线，进行莫队，每增加一只袜子，组合数就会增加原有袜子数。然后对区间也组合一下。再进行gcd一下就好了。

/*
@resources: HYSBZ 2038
@date: 2017-09-04
@author: QuanQqqqq
@algorithm: 莫队算法 
*/
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <math.h>

#define MAXN 50005
#define ll long long

using namespace std;

struct Question {
    int l, r, idx;  
} qsn[MAXN];

ll ans;
int pos[MAXN];
int flag[MAXN];
int a[MAXN];
ll lans[MAXN][2];

void add(int x) {
    ans += flag[a[x]];  
    flag[a[x]]++;
}

void del(int x) {
    flag[a[x]]--;
    ans -= flag[a[x]];
}

int cmp(Question a, Question b) {
    if (pos[a.l] != pos[b.l]) {
        return pos[a.l] < pos[b.l];
    }
    return a.r < b.r;
}

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

int main() {
    int n, q;
    while (~scanf("%d %d", &n, &q)) {
        memset(flag, 0, sizeof(flag));
        int k = sqrt(n);
        for (int i = 1; i <= n; i++) {
            scanf("%d", &a[i]);
            pos[i] = i / k;
        }
        for (int i = 1; i <= q; i++) {
            scanf("%d %d", &qsn[i].l, &qsn[i].r);
            qsn[i].idx = i;
        }
        sort(qsn + 1, qsn + q + 1, cmp);
        int l = 1, r = 0;
        ans = 0;
        for (int i = 1; i <= q; i++) {
            while (qsn[i].l < l) {
                add(--l);
            }
            while (qsn[i].l > l) {
                del(l++);
            }
            while (qsn[i].r < r) {
                del(r--);
            }
            while (qsn[i].r > r) {
                add(++r);
            }
            ll len = qsn[i].r - qsn[i].l + 1;
            len = len * (len - 1) >> 1;
            ll p = gcd(ans, len);
            lans[qsn[i].idx][0] = ans / p;
            lans[qsn[i].idx][1] = len / p;
        }
        for (int i = 1; i <= q; i++) {
            printf("%lld/%lld\n", lans[i][0], lans[i][1]);
        }
    }
}