51nod 1290 Counting Diff Pairs(莫队算法)

弱鸡表示看完题后，除了暴力计算移动L和R带来的影响，并没想到什么好法子。遂找大佬题解：http://blog.csdn.net/zsgg_acm/article/details/50262735
是通过树状数组处理的。
因为数据范围是1e9，所以要先离散化，不然没法处理。。离散化的时候载坑里了，调了老久。
举个例子：
假设现在的序列是:1,2,3，k=3，树状数组c全为0，要查询整个区间[1,3]内有多少对数字的差的绝对值<=k，此时L=1,R=0,L=l，不需要移动，R < r，右移，移到位置1，查询此时[1-3,1+3]区间和，res = 0，这时候加上1的贡献，c[1] = 1，然后R++，计算[2-3,2+3]区间和，res=1，这里有一对，然后再加上2的贡献，c[2]=1，然后R++，计算[3-3,3+3]的区间和，res=2，这里有两对，再加上3的贡献，c[3]=1，R==3，结束。这个询问共3对符合条件的数字。这组序列是连续，在题目中的数据，离散化之后也都是连续的了。还有1-3这些负数，离散化之后也都是正的了。
感觉我的叙述还是比较烂的，还是代码比较直观。

#include <bits/stdc++.h>using namespace std;typedef long long LL;const int MAXN = 5e5+10;pair<int,int> rec[MAXN];vector<int> vec;int num[MAXN],cur[MAXN];int Belong[MAXN];int block;LL res[MAXN];int n,k,q;struct Query{    int l,r,id;} Q[MAXN];bool cmp(const Query& a, const Query& b){    if(Belong[a.l] == Belong[b.l])        return a.r < b.r;    return Belong[a.l] < Belong[b.l];}int c[MAXN];void add(int i, int x){    while(i <= n)    {        c[i] += x;        i += (i&-i);    }}//计算区间和int query(int l, int r){    int s = 0;    while(r > 0)    {        s += c[r];        r -= (r&-r);    }    --l;    while(l > 0)    {        s -= c[l];        l -= (l&-l);    }    return s;}void solve(){    int L = 1, R = 0;    LL temp = 0;    for(int i = 1; i <= q; ++i)    {        while(L < Q[i].l)        {            add(cur[L],-1);            temp -= query(rec[L].first,rec[L].second);            ++L;        }        while(L > Q[i].l)        {            --L;            temp += query(rec[L].first,rec[L].second);            add(cur[L],1);        }        while(R < Q[i].r)        {            ++R;            temp += query(rec[R].first,rec[R].second);            add(cur[R],1);        }        while(R > Q[i].r)        {            add(cur[R],-1);            temp -= query(rec[R].first,rec[R].second);            --R;        }        res[Q[i].id] = temp;    }}//传来的flag主要用于判定是向左边的扩展，还是向右边的//向右边的扩展x，如果在数组中不存在的话，就会返回比他大的那个数字的迭代器，那是不符合条件的int tn;int getid(int x, int flag){    int id = lower_bound(vec.begin(),vec.end(),x)-vec.begin();    if(id == tn) return tn;    if(id == 0) return 1;    if(flag)        return id+1;    else    {        if(vec[id] == x) return id+1;        return id;    }}int main(){    scanf("%d %d %d",&n,&k,&q);    block = sqrt(n);    for(int i = 1; i <= n; ++i)    {        scanf("%d",&num[i]);        vec.push_back(num[i]);        Belong[i] = (i-1)/block+1;    }    sort(vec.begin(),vec.end());    vec.erase(unique(vec.begin(),vec.end()),vec.end());    tn = vec.size();    for(int i = 1; i <= n; ++i)    {        rec[i].first = getid(num[i]-k,1);        rec[i].second = getid(num[i]+k,0);        cur[i] = getid(num[i],1);    }    n = vec.size();    for(int i = 1; i <= q; ++i)    {        scanf("%d %d",&Q[i].l,&Q[i].r);        ++Q[i].l,++Q[i].r;        Q[i].id = i;    }    sort(Q+1,Q+1+q,cmp);    solve();    for(int i = 1; i <= q; ++i)        printf("%I64d\n",res[i]);    return 0;}