莫队算法学习小记 bzoj3781

所谓莫队，就是。。。好吧我说不出来具体释义。。
这个算法常用于序列的处理中，对于序列处理有很大的作用。
能够使用莫队算法的题目都具有一个特点：
对于每一个询问的区间【l,r】都能O(1)推导到l-1,r.l,r+1.l+1,r.l,r-1.
使用莫队时要先对所有的询问分块，每一个块的长度为sqrt（n），然后按块的顺序排序，询问的r作为第二关键字（具体看题目，一般来说是以r为第二关键字）。
是一个离线算法。
下面给出一道例题：bzoj3781（权限题）
Description

小B有一个序列，包含N个1~K之间的整数。他一共有M个询问，每个询问给定一个区间[L..R]，求Sigma(c(i)^2)的值,其中i的值从1到K，其中c(i)表示数字i在[L..R]中的重复次数。小B请你帮助他回答询问。
Input

第一行，三个整数N、M、K。
第二行，N个整数，表示小B的序列。
接下来的M行，每行两个整数L、R。
Output

M行，每行一个整数，其中第i行的整数表示第i个询问的答案。
Sample Input

6 4 3
1 3 2 1 1 3
1 4
2 6
3 5
5 6
Sample Output

6
9
5
2
HINT

对于全部的数据，1<=N、M、K<=50000
这就是到裸莫队，每一次推倒的时候可以知道(n+1)^2-n^2=2n+1 ,然后每一次ans+=2n+1就可以了，有时不一定要从头开始，可以先预处理l[1],r[1]，然后再推导。代码（鲜明易懂）

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
#define fo(i,a,b) for (int i=a;i<=b;i++)
const int maxn=1e6;
#define ll long long
struct node{
    int l,r,id,bl;
};
int n,m,k;
int a[maxn];
ll res[maxn],bz[maxn];
node b[maxn];
int l,r;
bool cmp(node i,node j)
{
    return (i.bl<j.bl)||((i.l==j.l)&& (i.r<j.r));
}
int main()
{
    scanf("%d%d%d",&n,&m,&k);
    int size=sqrt(n);
    fo(i,1,n) scanf("%d",&a[i]);
    fo(i,1,m) scanf("%d%d",&b[i].l,&b[i].r),b[i].id=i,b[i].bl=(b[i].l-1)/size+1;
    sort(b+1,b+1+m,cmp);
    int l=1;
    int r=0;
    ll ans=0;
    fo(i,1,m)
    {
        while (l>b[i].l)l--,bz[a[l]]++,ans+=2*bz[a[l]]-1;
        while (r<b[i].r)r++,bz[a[r]]++,ans+=2*bz[a[r]]-1;
        while (l<b[i].l)bz[a[l]]--,ans-=2*bz[a[l]]+1,l++;
        while (r>b[i].r)bz[a[r]]--,ans-=2*bz[a[r]]+1,r--;
        res[b[i].id]=ans;
    }
    fo(i,1,m)printf("%lld\n",res[i]);
    return 0;
}