[bzoj3289]Mato的文件管理

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

Sample Output

HINT

n,q <= 50000

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

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

题解

莫队算法，每次处理询问时维护区间的动态逆序对

#include<cmath>

#include<cstdio>

#include<algorithm>

#define buf 50001

using namespace std;

int F(){

    int x=,c=getchar(),f=;

    for(;c<||c>;c=getchar())

        if(!(c^))

            f=-;

    for(;c>&&c<;c=getchar())

        x=(x<<)+(x<<)+c-;

    return x*f;

}

struct ins{

    int l,r,blg,id;

    bool operator<(const ins h)const{

        if(!(blg^h.blg))

            return r<h.r;

        return blg<h.blg;

    }

}s[buf];

int n,q,t,f[buf],d[buf];

unsigned int ans[buf],c[buf],now;

inline void add(int x,int d){

    for(;x<=n;x+=x&-x)

        c[x]+=d;

}

inline unsigned int sum(int x){

    unsigned int res=;

    for(;x;x-=x&-x)

        res+=c[x];

    return res;

}

void solve(){

    sort(s+,s++q);

    int l=,r=;

    for(int i=;i<=q;i++){

        for(;l<s[i].l;l++)

            add(d[l],-),now-=sum(d[l]-);

        for(;r>s[i].r;r--)

            add(d[r],-),now-=r-l-sum(d[r]);

        for(;l>s[i].l;)

            l--,add(d[l],),now+=sum(d[l]-);

        for(;r<s[i].r;)

            r++,add(d[r],),now+=r-l+-sum(d[r]);

        ans[s[i].id]=now;

    }

}

int main(){

    n=F();

    t=sqrt(n);

    for(int i=;i<=n;i++)

        f[i]=d[i]=F();

    sort(f+,f++n);

    for(int i=;i<=n;i++)

        d[i]=lower_bound(f+,f++n,d[i])-f;

    q=F();

    for(int i=;i<=q;i++)

        s[i].id=i,

        s[i].l=F(),

        s[i].r=F(),

        s[i].blg=(s[i].l-)/t+;

    solve();

    for(int i=;i<=q;i++)

        printf("%d\n",ans[i]);

    return ;

}