BZOJ 2038 小Z的袜子（莫队算法）

莫队算法的话，看这个比较好懂：
http://blog.csdn.net/bossup/article/details/39236275
不过这里面的博客的代码貌似有点问题，所以代码的话，
http://blog.csdn.net/lwt36/article/details/50583757
这个是对的，看我的也行。
什么时候用莫队？
1、区间不会被修改
2、可以离线
3、对于[L,R]，可以直接退出[L,R+1],[L,R-1],[L+1,R],[L-1,R]
4、数据规模合适，莫队的复杂度nsqrt(n)

题意：

题意我就不说了，其实这题对于公式的化简也算是一个难点，不过，都把他
当做莫队的模板题，所以也就忽视了推导过程。
对于L,R的询问。设其中颜色为x,y,z....的袜子的个数为a,b,c。。。
那么答案即为(a*(a-1)/2+b*(b-1)/2+c*(c-1)/2....)/((R-L+1)*
(R-L)/2)化简得:(a^2+b^2+c^2+...x^2-(a+b+c+d+.....))/
((R-L+1)*(R-L))
即：(a^2+b^2+c^2+...x^2-(R-L+1))/((R-L+1)*(R-L))
这就是推导过程，恩，看了就懂了。

莫队有一个关键，就是块内的r是上升的，一定要记住，不然和暴力
就没区别了。
我已开始排序的时候写成这样：

struct Query
{
    int l,r,id;
    Query(int l,int r,int id):l(l),r(r),id(id){}
    bool operator < (const Query& rhs) const {
        return l==rhs.l?r.r:l.l;
    }
};

结果T到死。
代码：

//
//  Created by  CQU_CST_WuErli
//  Copyright (c) 2016 CQU_CST_WuErli. All rights reserved.
//
// #include
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#define CLR(x) memset(x,0,sizeof(x))
#define OFF(x) memset(x,-1,sizeof(x))
#define MEM(x,a) memset((x),(a),sizeof(x))
#define For_UVa if (kase!=1) cout << endl
#define BUG cout << "I am here" << endl
#define lookln(x) cout << #x << "=" << x << endl
#define SI(a) scanf("%d",&a)
#define SII(a,b) scanf("%d%d",&a,&b)
#define SIII(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define rep(flag,start,end) for(int flag=start;flag<=end;flag++)
#define Rep(flag,start,end) for(int flag=start;flag>=end;flag--)
#define Lson l,mid,rt<<1
#define Rson mid+1,r,rt<<1|1
#define Root 1,n,1
#define BigInteger bign
const int MAX_L=2005;// For BigInteger
const int INF_INT=0x3f3f3f3f;
const long long INF_LL=0x7fffffff;
const int MOD=1e9+7;
const double eps=1e-9;
const double pi=acos(-1);
typedef long long  ll;
using namespace std;

inline int read() {
    char c = getchar();
    while(!isdigit(c)) c = getchar();
    int x = 0;
    while(isdigit(c)) {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x;
}

const int N=5e4+10;
const int M=250;

int n,m;
int c[N];
int num[N];
ll ans,up[N],dw[N];

struct Query
{
    int l,r,id;
    Query(int l,int r,int id):l(l),r(r),id(id){}
    bool operator < (const Query& rhs) const {
        return rvector block[M];

void update(int x,int d) {
    ans-=num[c[x]]*num[c[x]];
    num[c[x]]+=d;
    ans+=num[c[x]]*num[c[x]];
}

int main() {
#ifdef LOCAL
    freopen("C:\\Users\\john\\Desktop\\in.txt","r",stdin);
    // freopen("C:\\Users\\john\\Desktop\\out.txt","w",stdout);
#endif
    while(SII(n,m)==2) {
        rep(i,0,250) block[i].clear();
        rep(i,0,n-1) c[i]=read();
        rep(i,1,m) {
            int l,r;
            l=read();r=read();
            l--,r--;
            block[l/M].push_back(Query(l,r,i));
        }
        ans=0LL;
        CLR(num);
        rep(i,0,n/M+2) sort(block[i].begin(),block[i].end());
        int pL,pR;
        rep(i,0,n/M+2) if (block[i].size()) {
            pL=i*M;
            pR=i*M-1;
            rep(j,0,block[i].size()-1) {
                Query &q=block[i][j];
                if (pR1,q.r) update(k,1);
                else Rep(k,pR,q.r+1) update(k,-1);
                pR=q.r;
                if (pL1) update(k,-1);
                else Rep(k,pL-1,q.l) update(k,1);
                pL=q.l;
                up[q.id]=ans-q.r+q.l-1;
                dw[q.id]=1LL*(q.r-q.l)*(q.r-q.l+1);
                // cout << up[q.id] << ' ' << dw[q.id] << endl;
                ll tmp=__gcd(up[q.id],dw[q.id]);
                up[q.id]/=tmp;
                dw[q.id]/=tmp;
            }
            rep(j,pL,pR) update(j,-1);
        }
        rep(i,1,m) printf("%lld/%lld\n",up[i],dw[i]);
    }
    return 0;
}