分治——求排列的逆序数

问题描述

• 考虑1,2,…,n(n<=100000)的排列$i_1,i_2,...,i_n$，如果其中存在j，k，满足$j<k$且$i_j>i_k$，那么就称$(i_j,i_k)$是这个排列的一个逆序。
• 一个排列含有的逆序的个数成为这个排列的逆序数。例如排列263451含有8个逆序(2,1),(6,3),(6,4),(6,5),(6,1),(3,1),(4,1),(5,1)，因此该排列的逆序数就是8.
• 现给定1,2,…,n的一个排列，求它的逆序数。

题目解析

• 1、笨办法：$\mathrm{O}(n^2)$
• 2、分治$\mathrm{O}(nlogn)$:
  • 1）将数组分成两半，分别求出左半边的逆序数和右半边的逆序数
  • 2）再算有多少逆序是由左半边取一个数和右半边取一个数构成(要求$\mathrm{O}(n)$实现)。

注：2）的关键：左半边和右半边都是排好序的。比如，都是从大到小排序的。这样，左右半边只需要从头到尾各扫一遍，就可以找出由两边各取一个数构成的逆序个数。

• 如下图所示：

每当al[i]>ar[j]时则说明存在逆序数，逆序数的个数为ar[j]及其以后的数：len2-j

• C语言代码示例

#include<stdio.h>

int mergeSortAndCount(int a[], int left, int mid, int right);

int merge(int a[], int left, int right)
{
  int mid = left + (right-left)/2;
  int count = 0;
  if(left < right)
  {
     count += merge(a, left, mid);
     count += merge(a, mid+1, right);
     count += mergeSortAndCount(a, left, mid, right);
  }
  return count;
}
int mergeSortAndCount(int a[], int left, int mid, int right)
{
  int len1 = mid - left + 1;
  int len2 = right - mid;
  int al[len1], ar[len2], i, j, k;
  int count = 0; //逆序数
  for(i = 0; i < len1; i++)
    al[i] = a[left+i];
  for(i = 0; i < len2; i++)
    ar[i] = a[mid+1+i];
  i = 0;
  j = 0;
  for(k = left; k <= right; k++)
  {
     if(i >= len1)
     {
       a[k] = ar[j];
       j++;
     }else if(j >= len2)
     {
       a[k] = al[i];
       i++;
     }else if(al[i] > ar[j]){
       a[k] = al[i];
       i++;
       count += len2-j;
     }else{
       a[k] = ar[j];
       j++;
     }
  }
  return count;
}

int main()
{
  int a[100];
  int i = 0, n;
  printf("请输入数组的长度：");
  scanf("%d", &n);
  printf("请输入%d个整数：\n", n);
  for(i=0; i< n; i++)
     scanf("%d", &a[i]);
  printf("逆序数为：%d\n", merge(a, 0, n-1));
  return 0;
}

• 运行示例

小结

• 通过此次实验，更加了解归并排序算法的实现过程，同时也解决了个人的思维定式。在求取逆序数个数的时候，没想到“ar[j]后的数都满足逆序”这个点，以至于困扰了一会儿。。。。。。

• 注：文中问题及代码参考 MOOC——《程序设计与算法》（北京大学 郭炜）

• 算法系列目录：

  • 一、递归形式的问题
    • 汉诺塔问题（Hanoi)
    • N皇后问题
    • 逆波兰表达式
    • 四则运算表达式求值
    • 爬楼梯
    • 放苹果
    • 算24
  • 二、程序或算法的时间复杂度
  • 三、二分算法
    • 二分查找
    • 二分法求方程的根
    • 寻找指定和的整数对
  • 四、分治算法
    • 归并排序、快速排序
    • 输出前m大的数
    • 求排列的逆序数