BZOJ 3289: Mato的文件管理

3289: Mato的文件管理

Time Limit: 40 Sec  Memory Limit: 128 MB
Submit: 2368  Solved: 971
[Submit][Status][Discuss]

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]区间的文件。

Output

q行，每行一个正整数，表示Mato这天需要交换的次数。

Sample Input

4
1 4 2 3
2
1 2
2 4

Sample Output

0
2

HINT

Hint

n,q <= 50000

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

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

Source

By taorunz

[Submit][Status][Discuss]

小生的第二道莫队算法。易知一段区间内的最小交换次数就是这段区间内的逆序对数，这个可以通过用树状数组动态维护。树状数组是按权值组织的，当在区间右侧新加入一个元素时，新形成的逆序对数就是前面区间中严格大于该元素的元素个数，反之亦然，请客官自行推导吧(ﾉ*･ω･)ﾉ。

#include <bits/stdc++.h>

template <class T>
inline void read(T &num) {
    register int neg = false;
    register int bit = getchar();
    
    while (bit <= '0') {
        if (bit == '-')
            neg ^= neg;
        bit = getchar();
    }
    
    num = 0;
    
    while (bit >= '0') {
        num = num*10 
        + bit - '0';
        bit = getchar();
    }
    
    if (neg)num = -num;
}

const int siz = 50005;

int n, m, s, num[siz], tmp[siz], bit[siz], answer;

struct query {
    int l, r, id, ans;
}qry[siz];

inline bool cmp_lr(const query &A, const query &B) {
    if (A.l / s != B.l / s)
        return A.l < B.l;
    else
        return A.r < B.r;
}

inline bool cmp_id(const query &A, const query &B) {
    return A.id < B.id;
}

inline void add(int k, int v) {
    for (; k <= n; k += k&-k)bit[k] += v;
}

inline int ask(int k) {
    int ret = 0;
    for (; k >= 1; k -= k&-k)ret += bit[k];
    return ret;
}

signed main(void) {
    read(n);
    
    for (int i = 1; i <= n; ++i)
        read(num[i]);
        
    memcpy(tmp, num, sizeof(tmp));
    
    std::sort(tmp + 1, tmp + 1 + n);
    
    int tot = std::unique(tmp + 1, tmp + 1 + n) - tmp;
    
    for (int i = 1; i <= n; ++i)
        num[i] = std::lower_bound(tmp, tmp + tot, num[i]) - tmp;
    
    read(m);
    
    for (int i = 1; i <= m; ++i)
        read(qry[i].l), read(qry[i].r), qry[i].id = i;
        
    s = sqrt(n);
    
    std::sort(qry + 1, qry + 1 + m, cmp_lr);
    
    for (int i = 1, x = 1, y = 0; i <= m; ++i) {
        while (x < qry[i].l)answer -= ask(num[x] - 1), add(num[x], -1), ++x;
        while (x > qry[i].l)--x, answer += ask(num[x] - 1), add(num[x], +1);
        while (y < qry[i].r)++y, answer += ask(n) - ask(num[y]), add(num[y], +1);
        while (y > qry[i].r)answer -= ask(n) - ask(num[y]), add(num[y], -1), --y;
        qry[i].ans = answer;
    }
    
    std::sort(qry + 1, qry + 1 + m, cmp_id);
    
    for (int i = 1; i <= m; ++i)
        printf("%d
", qry[i].ans);
}

@Author: YouSiki