bzoj3289 mato的文件管理

3289: Mato的文件管理

Time Limit: 40 Sec  Memory Limit: 128 MB
Submit: 3007  Solved: 1272
[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次移到最后。

 

————————————————————————————————————————————————————————————————————————————

思路：

　　交换次数也就是逆序对的个数。

　　明显莫队并以树状数组或线段树来维护逆序对的个数。莫队每增加或减少一个值时，计算增加或减少逆序对的个数，OK!

　　道理很简单，但是我用了多种线段树的写法，都TLE，树状数组一次过了！

　　哪位大神看一下我的线段树，那里导致的超时！感谢！

————————————————————————————————————————————————————————————————————————————

 

#include<bits/stdc++.h>
  
using namespace std;
const int maxn=5e4+7;
typedef long long ll;
int n,m;
int zl[maxn],hash[maxn];
struct que
{
    int l,r,id,belong;
}q[maxn];
ll da[maxn],ans;
int l,r,siz,xx;
int sz[maxn];
int dat[maxn<<2],lll[maxn<<2],rr[maxn<<2];
int js=0;
bool cmp(que a,que b)
{
    return a.belong<b.belong || (a.belong==b.belong && a.r<b.r);
}
void readint(int &x)
{
    char c=getchar();
    int f=1;
    for(;c>'9'||c<'0';c=getchar())if(c=='-')f=-f;
    x=0;
    for(;c<='9' && c>='0';c=getchar())x=(x<<1)+(x<<3)+c-'0';
    x=x*f;
}
inline int lowbit(int x)
{
    return x&(-x);
}
int fin(int ps)
{
    int tp=0;
    for(;ps>0;ps-=lowbit(ps))tp+=sz[ps];
    return tp;
}
void update(int ps,int x)
{
    for(;ps<=siz;ps+=lowbit(ps))sz[ps]+=x;
}
void add0(int ps)
{
    ps=zl[ps];
    ans+=fin(ps-1);
    update(ps,1);
}
void del0(int ps)
{
    ps=zl[ps];
    ans-=fin(ps-1);
    update(ps,-1);
}
void add1(int ps)
{
    ps=zl[ps];
    ans+=xx-fin(ps);
    update(ps,1);
}
void del1(int ps)
{
    ps=zl[ps];
    ans-=xx-fin(ps);
    update(ps,-1);
}
int main()
{
    readint(n);
    for(int i=1;i<=n;i++)
    {
        readint(zl[i]);hash[i]=zl[i];
    }
    sort(hash+1,hash+1+n);
    siz=unique(hash+1,hash+1+n)-hash-1;
    for(int i=1;i<=n;i++)
    {
        zl[i]=lower_bound(hash+1,hash+1+siz,zl[i])-hash;
    }
    readint(m);
    int bks=sqrt(n+0.5);
    for(int i=0;i<m;i++)
    {
        readint(q[i].l);readint(q[i].r);
        q[i].id=i;q[i].belong=(q[i].l-1)/bks;
    }
    sort(q,q+m,cmp);
    memset(sz,0,sizeof(sz));
    l=2;r=1;ans=0;
    for(int i=0;i<m;i++)
    {
        while(l<q[i].l)
        {
            del0(l);l++;xx=r-l+1;
        }
        while(l>q[i].l)
        {
            l--;add0(l);xx=r-l+1;
        }
        while(r<q[i].r)
        {
            r++;add1(r); xx=r-l+1;
        } 
        while(r>q[i].r)
        {
            del1(r);r--;xx=r-l+1;
        }
        da[q[i].id]=ans;
    }
    for(int i=0;i<m;i++)printf("%d\n",da[i]);
    return 0;
}

下面是线段树没过的，请大神指点！

#include<bits/stdc++.h>
  
using namespace std;
const int maxn=5e4+7;
typedef long long ll;
int n,m;
int zl[maxn],hash[maxn];
struct que
{
    int l,r,id,belong;
}q[maxn];
ll da[maxn],ans;
int l,r,siz;
int dat[maxn<<2],lll[maxn<<2],rr[maxn<<2];
int js=0;
bool cmp(que a,que b)
{
    return a.belong<b.belong || (a.belong==b.belong && a.r<b.r);
}
void readint(int &x)
{
    char c=getchar();
    int f=1;
    for(;c>'9'||c<'0';c=getchar())if(c=='-')f=-f;
    x=0;
    for(;c<='9' && c>='0';c=getchar())x=(x<<1)+(x<<3)+c-'0';
    x=x*f;
}
void build(int cur,int l,int r)
{
    dat[cur]=0;
    lll[cur]=l;rr[cur]=r;
    if(l==r)return ;
    int mid=(l+r)/2;
    build(cur<<1,l,mid);
    build((cur<<1)|1,mid+1,r);  
}
int fin(int cur,int l,int r)
{
    if(r<l)return 0;
    if(l<=lll[cur] && rr[cur]<=r)return dat[cur];
    int mid=(lll[cur]+rr[cur])>>1;
    int tp=0;
    if(l<=mid)tp+=fin(cur<<1,l,r);
    if(r>mid)tp+=fin((cur<<1)|1,l,r);
    return tp;
}
void update(int cur,int ps,int x)
{
    dat[cur]+=x;
    if(lll[cur]==rr[cur])return ;
    int mid=(lll[cur]+rr[cur])>>1;
    if(ps<=mid)update(cur<<1,ps,x);
    else update((cur<<1)|1,ps,x); 
}
void add0(int ps)
{
    ps=zl[ps];
    ans+=fin(1,1,ps-1);
    update(1,ps,1);
}
void del0(int ps)
{
    ps=zl[ps];
    ans-=fin(1,1,ps-1);
    update(1,ps,-1);
}
void add1(int ps)
{
    ps=zl[ps];
    ans+=fin(1,ps+1,siz);
    update(1,ps,1);
}
void del1(int ps)
{
    ps=zl[ps];
    ans-=fin(1,ps+1,siz);
    update(1,ps,-1);
}
int main()
{
    readint(n);
    for(int i=1;i<=n;i++)
    {
        readint(zl[i]);hash[i]=zl[i];
    }
    sort(hash+1,hash+1+n);
    siz=unique(hash+1,hash+1+n)-hash-1;
    for(int i=1;i<=n;i++)
    {
        zl[i]=lower_bound(hash+1,hash+1+siz,zl[i])-hash;
    }
    readint(m);
    int bks=sqrt(n+0.5);
    for(int i=0;i<m;i++)
    {
        readint(q[i].l);readint(q[i].r);
        q[i].id=i;q[i].belong=(q[i].l-1)/bks;
    }
    sort(q,q+m,cmp);
    build(1,1,siz);
    l=2;r=1;ans=0;
    for(int i=0;i<m;i++)
    {
        while(l<q[i].l)
        {
            del0(l);l++;
        }
        while(l>q[i].l)
        {
            l--;add0(l);
        }
        while(r<q[i].r)
        {
            r++;add1(r); 
        } 
        while(r>q[i].r)
        {
            del1(r);r--;
        }
        da[q[i].id]=ans;
    }
    for(int i=0;i<m;i++)printf("%d\n",da[i]);
    return 0;
}