Minimum Inversion Number(线段树求逆序数小总结)

首先,求逆序数对的思路:
1.得到整个数列后,从前往后扫,统计比a[i]小的,在a[i]后面的有多少个
这样做的话,应该是只有n2的暴力作法,没想到更好的方法
2.统计a[i]前面的,且比它大的数
这样做的话,就可以利用输入的时效性,每输入一个数,就把这个数的num[i]值加1,
然后统计比这个数大的数的num和,
因为这里的和一定是在这个数列中比a[i]大,且在它前面出现的数之和,
然后把把这个和加到总逆序数sum里。
这样做的话直接的暴力作法依然是n2,但是,
我们可以在, 统计 比这个数大的数的num和这一步进行优化,利用线段树求区间域值的复杂度是logn,
所以总体复杂度就降到了nlogn。

再来看这道题,求得初始数列的逆序数后,再求其他排列的逆序数有一个规律,就是
sum = sum + (n - 1 - a[i]) - a[i];
这个自行验证吧,相信很容易得出

最后,拓展一下,如果要求正序数怎么办?很简单,无非是大小调一下
再问,如果要求满足i<j<k,且a[i]>a[j]>a[k]的数对总数怎么办?

可以从中间的这个数入手,统计a[i]>a[j]的对数m,以及a[j]>a[k]的对数n,m*n就是。。。
要求a[i]>a[j]的个数还是一样的,那么a[j]>a[k]的个数呢?
两种思路:
1.得到a[i]>a[j]的对数后,将数列倒过来后再求a[j]<a[k]的对数
2.更简单的做法是,找到规律发现,n = 整个数列中比a[j]小的数 — 在a[j]前面已经出现的比a[j]小的数的个数
即(假设数列是从1开始的) n = (a[j] -1) - (j - 1 - m )
 
 
 
 
代码
#include <stdio.h>
#include <algorithm>

using namespace std;

int a[5005];
struct Node{
    int l,r,num;
}tree[50000];

void Build(int n,int x,int y){
    tree[n].l = x;
    tree[n].r = y;
    tree[n].num = 0;
    if(x == y){
        return;
    }
    int mid = (x + y) / 2;
    Build(2*n,x,mid);
    Build(2*n+1,mid+1,y);
}

void Modify(int n,int x){
    int l = tree[n].l;
    int r = tree[n].r;
    int mid = (l + r) / 2;
    if(x == l && x == r){
        tree[n].num = 1;
        return;
    }
    if(x <= mid)    Modify(2*n,x);
    else            Modify(2*n+1,x);
    tree[n].num = tree[2*n].num + tree[2*n+1].num;
}

int Query(int n,int x,int y){
    int l = tree[n].l;
    int r = tree[n].r;
    int mid = (l + r) / 2;
    int ans = 0;;
    if(x == l && y == r)
        return tree[n].num;
    if(x <= mid)   ans += Query(2*n,x,min(mid,y));
    if(y > mid)    ans += Query(2*n+1,max(mid+1,x),y);
    return ans;
}
int main(){
    int n,sum,ans;
    int i,j;

    while(scanf("%d",&n) != EOF){
        sum = 0;
        Build(1,0,n);
        for(i = 1;i <= n;i++){
            scanf("%d",&a[i]);
            Modify(1,a[i]);
            sum += Query(1,a[i]+1,n);
        }
        ans = sum;
        for(i = 1;i < n;i++){
            sum = sum + (n - 1 - a[i]) - a[i];
            if(sum < ans)
                ans = sum;
        }
        printf("%d\n",ans);
    }
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值