BZOJ 3289：Mato的文件管理（莫队算法+树状数组）

http://www.lydsy.com/JudgeOnline/problem.php?id=3289

题意：……

思路：求交换次数即求逆序对数。确定了这个之后，先离散化数组。然后在后面插入元素的话，就是在区间里面找比它大的元素数量，在前面插入元素的话，就是在区间里面找比它小的元素数量。删除操作类似。因为排序是从小到大排序，所以要找比它大的数量就是区间长度减去小于等于该元素的数量，所以是（R - L + 1 - sum(a[i])），要找比它小的数量就是（sum(a[i] - 1)）。然后用莫队算法处理区间询问就好了。

 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #include <cmath>
 using namespace std;
 #define N 50010
 struct node {
     int l, r, id, ans;
 } p[N];
 int bit[N], a[N], b[N], num[N], ans, kuai, n;

 bool cmp(const node &a, const node &b) {
     if(a.l / kuai == b.l / kuai) return a.r < b.r;
     return a.l / kuai < b.l / kuai;
 }
 bool cmpid(const node &a, const node &b) { return a.id < b.id; }
 int lowbit(int x) { return x & (-x); }
 void add(int x, int w) { for( ; x <= n; x += lowbit(x)) bit[x] += w; }
 int sum(int x) { int ans = ; for( ; x ; x -= lowbit(x)) ans += bit[x]; return ans; }

 int main() {
     while(~scanf("%d", &n)) {
         for(int i = ; i <= n; i++) scanf("%d", &a[i]), b[i] = a[i];
         int q; scanf("%d", &q);
         for(int i = ; i <= q; i++) scanf("%d%d", &p[i].l, &p[i].r), p[i].id = i;

         memset(bit, , sizeof(bit));
         ans = ; kuai = sqrt(n);
         sort(b + , b +  + n);
         int cnt = unique(b + , b +  + n) - b - ;
         for(int i = ; i <= n; i++) a[i] = lower_bound(b + , b +  + cnt, a[i]) - b;

         sort(p + , p +  + q, cmp);
         for(int L = , R = , i = ; i <= q; i++) {
             int l = p[i].l, r = p[i].r;
             for( ; L < l; L++) { ans -= sum(a[L] - ); add(a[L], -); }
             for( ; L > l; L--) { ans += sum(a[L-] - ); add(a[L-], ); }
             for( ; R < r; R++) { ans += R - L +  - sum(a[R+]); add(a[R+], ); }
             for( ; R > r; R--) { ans -= R - L +  - sum(a[R]); add(a[R], -); }
             p[i].ans = ans;
         }
         sort(p + , p +  + q, cmpid);
         for(int i = ; i <= q; i++) printf("%d\n", p[i].ans);
     }
     return ;
 }