bzoj3289（莫队+树状数组）

Description
Mato同学从各路神犇以各种方式（你们懂的）收集了许多资料，这些资料一共有n份，每份有一个大小和一个编号
。为了防止他人偷拷，这些资料都是加密过的，只能用Mato自己写的程序才能访问。Mato每天随机选一个区间[l,r
]，他今天就看编号在此区间内的这些资料。Mato有一个习惯，他总是从文件大小从小到大看资料。他先把要看的
文件按编号顺序依次拷贝出来，再用他写的排序程序给文件大小排序。排序程序可以在1单位时间内交换2个相邻的
文件（因为加密需要，不能随机访问）。Mato想要使文件交换次数最小，你能告诉他每天需要交换多少次吗？
Input
第一行一个正整数n，表示Mato的资料份数。
第二行由空格隔开的n个正整数，第i个表示编号为i的资料的大小。
第三行一个正整数q，表示Mato会看几天资料。
之后q行每行两个正整数l、r，表示Mato这天看[l,r]区间的文件。
n,q <= 50000
Output
q行，每行一个正整数，表示Mato这天需要交换的次数。

Sample Input
4

1 4 2 3

2

1 2

2 4
Sample Output
0

2

//样例解释：第一天，Mato不需要交换

第二天，Mato可以把2号交换2次移到最后。
HINT
Source
By taorunz

---

经典题目。
每次删除和加入直接在树状数组里进行统计即可

#include<bits/stdc++.h>
using namespace std;
int a[50100] , n , Q , m;
int tmp[50100];
int tr[50100] , belong[50100];
unsigned int ans,out[50100];
struct data
{
    int n , id;
}da[50100];
struct question
{
    int l , r , id;
}q[50100];
int read()
{
    int sum = 0;char c = getchar();bool flag = true;
    while( c < '0' || c > '9' ) {if(c == '-') flag = false;c = getchar();}
    while( c >= '0' && c <= '9' ) sum = sum * 10 + c - 48 , c = getchar();
    if(flag)  return sum;
     else return -sum;
}  
unsigned int find(int x)
{
    unsigned int sum = 0;
    for(;x;x -= x&(-x)){
        sum += tr[x];
        //if(c==5)printf("#%d %d\n",x,tr[x]);
    } 
    return sum;
}
void add(int x,int v)
{
    for(;x <= m;x += x&(-x))
        tr[x] += v;
    return;
}
bool mycmp(data a,data b)
{
    return a.n < b.n;
}
bool mycmp1(question a,question b)
{
    return belong[a.l] < belong[b.l] || 
          (belong[a.l] == belong[b.l] && a.r < b.r); 
}
void init()
{
    n = read();; 
    for(int i = 1;i <= n;++i)
        da[i].n = read() , da[i].id = i;
    sort(da + 1,da + n + 1,mycmp);
    for(int i = 1;i <= n;++i)
        if(da[i].n != da[i - 1].n) tmp[i] = ++m;
            else tmp[i] = m;
    for(int i = 1;i <= n;++i)
        a[da[i].id] = tmp[i];
    int k = sqrt(n);
    for(int i = 1,j = k,t = 1;i <= n;i += k,j += k,++t)
        if(j <= n)   for(int x = i;x <= j;++x)    belong[x] = t;
        else   for(int x = i;x <= n;++x) belong[x] = t;
    Q = read();
    for(int i = 1;i <= Q;++i)
        q[i].l = read() , q[i].r = read() , q[i].id = i;
    sort(q + 1,q + Q + 1,mycmp1);
    return;
}
void work()
{
    for(int i = q[1].r;i >= q[1].l;--i)
        ans += find(a[i] - 1) , add(a[i] , 1);
    out[q[1].id] = ans;
    int l = q[1].l , r = q[1].r;
    for(int x = 2;x <= Q;++x)
    {
        if(l != q[x].l)
        if(l < q[x].l)
            for(int i = l;i < q[x].l;++i)
                ans -= find(a[i] - 1) , add(a[i],-1);
        else
            for(int i = l - 1;i >= q[x].l;--i)
                ans += find(a[i] - 1) , add(a[i],1);
        l = q[x].l;
        if(r != q[x].r)
        if(r < q[x].r)
            for(int i = r + 1;i <= q[x].r;++i)
                ans += i - l - find(a[i]) , add(a[i],1);
        else
            for(int i = r;i > q[x].r;--i)
                ans -= i - l + 1 - find(a[i]) , add(a[i],-1);
        r = q[x].r;
        out[q[x].id] = ans;
    }  
    for(int i = 1;i <= Q;++i) 
        printf("%d\n",out[i]);
    return;
}
int main()
{
    init();
    work();
    return 0;
}