莫队 入门（离线基础板子）

解决类型题：给定一个数列，进行区间的查询.

对于给定的数列，如果我们已经知道了一段区间要求的结果，设左右边界为L和R，然后让我们求qL到qR区间的结果ans，此时，我们可以让R和L移动，到qL和qR上，在移动的过程中，对ans进行修改即可，有点像尺取法对吧。我们要优化时间，就要尽量减少移动的次数，所以此时对整个数列进行分块，进行一个sort排序，如果左端点在一个块内，就以右端点排序，否则以左端点排序。至于时间复杂度，我也不是太理解，可能是因为不会分快吧。

例题：对于这些板子题，加减函数的ans计算还是比较灵活多变的

小Z的袜子(hose)bzoj 2038 

首先要知道，我们要求的是什么，随机对一段区间抽两双袜子，问是同一种颜色的概率为多少？还要化简分数，那我先确定一下分母是谁：long=(R-L+1)是区间长度，抽两双袜子，一共有多少种排列组合？long*(long-1)/2对吧。那我们可以用莫队算法来求一下分子：在这段区间中，某种颜色的袜子有x只我知道，所有颜色x的排列组合不就是分子嘛！！不过，x的变化跟分子并不是 y=x的关系，所以我有一个操作：先把颜色x的排列组合减去，等x 加减1之后，我再计算颜色x的排列组合，再把它加上去就okk了。

注意：区间已经进行了一个排序操作，由于最后还要求分母，我还要把顺序还原，所以就又有了一个乱七八糟的操作咯

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
using namespace std;
#define N 50000+20
struct node
{
    int l,r,num;
}q[N];
int block;
bool comp(node x,node y)
{
    return (x.l/block)==(y.l/block)?x.r<y.r:x.l<y.l;
}
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
int sum[N],cnt[N],a[N],n,m,id[N],ans;
int yao(int x)
{
    if(x%2)
        return (x-1)/2*x;
    return x/2*(x-1);
}
void add(int x)
{
    x=a[x];
    ans-=yao(cnt[x]);//减去原来的排列组合
    cnt[x]++;
    ans+=yao(cnt[x]);//加上数量变化后的排列组合
    return ;
}
void del(int x)
{
    x=a[x];
    ans-=yao(cnt[x]);
    cnt[x]--;
    ans+=yao(cnt[x]);
    return ;
}
int main()
{
    scanf("%d%d",&n,&m);
    memset(cnt,0,sizeof(cnt));
    block=(int)sqrt(n*1.0);
    for(int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&q[i].l,&q[i].r);
        q[i].num=i;
    }
    sort(q+1,q+1+m,comp);
    ans=0;
    int l=1,r=0,ql,qr;
    for(int i=1;i<=m;i++)
    {
        ql=q[i].l,qr=q[i].r;
        while(r<qr) add(++r);
        while(l<ql) del(l++);
        while(r>qr) del(r--);
        while(l>ql) add(--l);
        sum[q[i].num]=ans;
        id[q[i].num]=i;//乱七八糟的操作
    }
    for(int i=1;i<=m;i++)
    {
        int GCD=gcd(sum[i],yao(q[id[i]].r-q[id[i]].l+1));//乱七八糟的操作 
        printf("%d/%d\n",sum[i]/GCD,yao(q[id[i]].r-q[id[i]].l+1)/GCD);//求最大公因数
    }
    return 0;
}