求逆序数 树状数组

在一个排列中,如果一对数的前后位置与大小顺序相反,即前面的数大于后面的数,那么它们就称为一个逆序。一个排列中逆序的总数就称为这个排列的逆序数。


思路1:暴力查询,双重for循环

for(i=1;i<=n;i++) 
  for(j=i+1;j<=n;j++) 
      if(a[i]>a[j]) 
      { 
           k++; 
      }


如果数据吧比较多,更新和查找速度变得非常慢。


思路2:排序 
比如:2 4 3 1 
先查找最小的数1,在它前面有3个数据,全部都比他大,把这个1变成0 
再查找除去1最小的2,在它前面有0个数据比他大,然后把这个2变成0 
同理,查找3,他前面只有一个数据比他大(2已经变成0) 
所以,2 4 3 1这个数列的逆序数为3+0+1=4

如果数据吧比较多,更新和查找会比较繁琐。


思路3:结构体排序 + 树状数组
比如:2 4 3 1 
继续优化思路2: 
2 4 3 1(a) 
1 1 1 1(b) 
现在只需要查找当前数列中最小的数,查找他之前有多少个1,然后把当前位置的1变成0

2 4 3 1 
1 1 1 1 
1 2 1 4    //这个是把1放到树状数组里后的树状数组里的数
优化对1查找之后的累加 
利用树状数组,数组b初始化为1,求b的树状数组的前n-1项和,也就是sum(n-1)或者sum(n)-1,然后把n的位置减去1,也就是add(n,-1);

每次查找当前数列最小值,数据多的时候,查找就不那么方便,所以,我们利用结构体进行排序进行查找过程的优化

利用 t[].val 记录数据 
利用 t[].id   记录位置 
2 4 3 1 (t.val) 
1 2 3 4 (t.id) 
进行排序之后变成 
1 2 3 4 (t.val) 
4 1 3 2 (t.id ) 
1 1 1 1 ( b )//需要把b数组里的数放到树状数组 

然后我们就可以从a.val的第一个开始查找,依次向后移动,我们利用a.id记录位置,去判断和求和; 

比如,我们去判断a.val中的1前面有几个数,就可以用树状数组:sum(a.id)-1或者sum(a.id-1);把1变成0可以用add(a.id,-1);




题目链接: https://www.51nod.com/onlineJudge/submitList.html#!userId=13358&problemId=1019

这是A过的代码




#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cstring>
#define MAX 1000001
using namespace std;
long long n,k;
long long c[MAX],b[MAX];
struct node
{
    long long a,p;
} t[MAX];
bool cmp(node x,node y)
{
    return x.a<y.a;
}
int Lowbit(int n)
{
    return n&(n^(n-1));
}
int sum(int n)
{
    int ans=0;
    while(n>0)
    {
        ans+=c[n];
        n-=Lowbit(n);
    }
    return ans;
}
void add(int pos,int num)
{
    while(pos<=n)
    {
        c[pos]+=num;
        pos+=Lowbit(pos);
    }
}
int main()
{
    while(cin>>n)
    {
        for(int i=1; i<=n; i++)
        {
            cin>>t[i].a;
            t[i].p=i;
            b[i]=1;
        }
        for(int i=1; i<=n; i++)
        {
            c[i]=0;
            for(int j=i-Lowbit(i)+1; j<=i; j++)
                c[i]+=b[j];
            //cout<<c[i]<<" ";
        }
        //cout<<endl;
        sort(t+1,t+n+1,cmp);
        k=0;
        for(int i=1; i<=n; i++)
        {
            k+=sum(t[i].p)-1;
            add(t[i].p,-1);
        }
        cout<<k<<endl;
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值