[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

4
1 4 2 3
2
1 2
2 4

Sample Output

0
2

HINT

n,q <= 50000

样例解释:第一天,Mato不需要交换

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

题解

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

#include<cmath>
#include<cstdio>
#include<algorithm>
#define buf 50001
using namespace std;
int F(){
    int x=0,c=getchar(),f=1;
    for(;c<48||c>57;c=getchar())
        if(!(c^45))
            f=-1;
    for(;c>47&&c<58;c=getchar())
        x=(x<<1)+(x<<3)+c-48;
    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=0;
    for(;x;x-=x&-x)
        res+=c[x];
    return res;
}
void solve(){
    sort(s+1,s+1+q);
    int l=1,r=0;
    for(int i=1;i<=q;i++){
        for(;l<s[i].l;l++)
            add(d[l],-1),now-=sum(d[l]-1);
        for(;r>s[i].r;r--)
            add(d[r],-1),now-=r-l-sum(d[r]);
        for(;l>s[i].l;)
            l--,add(d[l],1),now+=sum(d[l]-1);
        for(;r<s[i].r;)
            r++,add(d[r],1),now+=r-l+1-sum(d[r]);
        ans[s[i].id]=now;
    }
}
int main(){
    n=F();
    t=sqrt(n);
    for(int i=1;i<=n;i++)
        f[i]=d[i]=F();
    sort(f+1,f+1+n);
    for(int i=1;i<=n;i++)
        d[i]=lower_bound(f+1,f+1+n,d[i])-f;
    q=F();
    for(int i=1;i<=q;i++)
        s[i].id=i,
        s[i].l=F(),
        s[i].r=F(),
        s[i].blg=(s[i].l-1)/t+1;
    solve();
    for(int i=1;i<=q;i++)
        printf("%d\n",ans[i]);
    return 0;
}

 

posted @ 2016-11-29 17:23  keshuqi  阅读(168)  评论(0编辑  收藏  举报