hdu-1394-Minimum Inversion Number(树状数组)

本文介绍了一道关于寻找序列中倒置数最小值的问题。通过对序列进行特定操作并利用树状数组来高效计算每种可能排列中的倒置数,从而找到最小的倒置数。

                                         Minimum Inversion Number

                                                                                               Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
                                                                                               Total Submission(s): 17105    Accepted Submission(s): 10390


Problem Description
The inversion number of a given number sequence a1, a2, ..., an is the number of pairs (ai, aj) that satisfy i < j and ai > aj.

For a given sequence of numbers a1, a2, ..., an, if we move the first m >= 0 numbers to the end of the seqence, we will obtain another sequence. There are totally n such sequences as the following:

a1, a2, ..., an-1, an (where m = 0 - the initial seqence)
a2, a3, ..., an, a1 (where m = 1)
a3, a4, ..., an, a1, a2 (where m = 2)
...
an, a1, a2, ..., an-1 (where m = n-1)

You are asked to write a program to find the minimum inversion number out of the above sequences.
 

Input
The input consists of a number of test cases. Each case consists of two lines: the first line contains a positive integer n (n <= 5000); the next line contains a permutation of the n integers from 0 to n-1.
 

Output
For each case, output the minimum inversion number on a single line.
 

Sample Input
  
10 1 3 6 9 0 8 5 7 4 2
 

Sample Output
16


题意:给定n个数的整数序列,每次可将序列首个数放到最后一位,其余数相对位置不变,通过这样的操作可以得到n种序列,要求这n种序列中倒置数的最小个数,即满足下标i<j,a[i]>a[j]。
题目链接:Minimum Inversion Number
解题思路:
先看一个序列倒置数的求法,输入一个数,然后就立刻判断此时前面出现比他大的数的个数,因为序列的数满足0~n,所以可以用树状数组c[p]表示其区间[p-2^k-1,p]对应下标数字出现的个数,那么区间a[p]~a[n]就是我们要求的此时(满足下标i>j)前面出现的比他大的数的个数,注意,树状数组的下标都是从1开始,所以要把输入的a+1。
然后就讨论不同序列的情况,当某一时刻首位a[i]被移动到最后一位,那么原序列由a[i]组成的倒置数有a[i]-1个(此时a[i]加过1),也就是后面比他小的个数(因为是首位,前面不需要考虑)。同理,加到最后一个后,他前面有n-a[i]个数可以和a[i]形成倒置数。
代码:
//Minimum Inversion Number
#include<cstdio>
#include<cstring>
#define N 5010
using namespace std;

int c[N],n,a[N];

int lowbit(int t)  //计算树状数组t区间管理的a的个数
{
    return t&(-t); //t 二进制数后面k个0,返回值为2^k
}
void updata(int p,int f) //p在树状数组最先出现的地方是c[p]
{
    while(p <= n){ //将后面所有包含a[p]的树状数组更新
        c[p] += f;
        p += lowbit(p); 
    }
}
int getsum(int p) //求和 a[1] 到 a[p]
{
    int s = 0;
    while(p > 0) {
        s += c[p];   //c[p] 只管理lowbit(p) 个数
        p -= lowbit(p); //那么接下来要算a[1] 到 a[p] - lowbit(p) 的区间和
    }
    return s;
}
void init()
{
    for(int i = 0;i <= n;i++) c[i] = 0;
}

int main()
{
    int i,s,mn;
    while(~scanf("%d",&n)){
        init();
        s = 0;
        for(i = 1;i <= n;i++){
            scanf("%d",&a[i]);
            a[i]++;                        //树状数组从下标1开始
            s += getsum(n) - getsum(a[i]);  //a[i]前面比它大的数,即可生成的倒置数
           // printf("s = %d \n",s);
            updata(a[i],1);                //让所有包含a[i]的树状数组+1
           // print();
        }
        mn = s;
        for(i = 1;i <= n;i++){
            s += n - a[i] - (a[i] - 1);//将a[i]从首位换到最后一位后,倒置数先减少 a[i] - 1(注意a[i]+1处理过了)
            if(s < mn) mn = s;         //然后又增加了 n-a[i],比较取最小的
        }
        printf("%d\n",mn);
    }
    return 0;
}


评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值