求数组A[1:n]的逆序对数, 如果i<j时 A[i]>A[j]则(i,j)则称为一个逆序对。
数组<2,3,8,6,1>的逆序对数就是5.
求逆序对数最简单的算法就是把数组的元素两两比较下,判断是否逆序对,复杂度为O(n^2).
有没有复杂度更低的算法呢?是否可以从排序算法的角度去思考逆序数?
排序的过程就是要把逆序的数的次序进行纠正的过程,复杂程度很有可能与逆序对的数量是有关的。
1. 逆序对的数量与数组中重复元素的数量是无关系的,因此,不稳定的排序方法如快排等不在考虑范围内
2. 稳定排序算法中,先看冒泡排序,每次进行冒泡交换数组元素时,相邻元素交换位置,不影响这两个元素与数组其它元素之间的顺序(逆序)关系,也就是说,每冒泡交换一次,数组逆序对数减一。冒泡交换的总次数就等于逆序对数。
3. 选择排序和冒泡排序相似,插入元素位置移动时不影响的是该元素之后的数组的逆序数对数量,而插入元素之前的数组是已经排好序的,逆序对为0,每次插入逆序对数的改变量等于元素原始位置与插入位置之间的偏移。
4.归并排序,逆序对数的改变发生在合并这一步。数组分成两半A=[A1,A2],A1,A2都是排序好的数组,合并时指向A1的指针为i,指向A2的指针为j,如果A2的元素均比A1的元素大,逆序对数量为0,如果A2的元素全比A1的小,那么逆序数为length(A2)*length(A1)。可以修改归并排序,当j++时,判断当前A1中还有多少个元素没有被归并,当前合并的A2元素产生的逆序对数就是多少,累积求和就可以搞定。
修改冒泡和选择排序求逆序数复杂度是O(n^2),利用归并排序求逆序对数的算法复杂度可以降低到O(nlogn)
数组<2,3,8,6,1>的逆序对数就是5.
求逆序对数最简单的算法就是把数组的元素两两比较下,判断是否逆序对,复杂度为O(n^2).
有没有复杂度更低的算法呢?是否可以从排序算法的角度去思考逆序数?
排序的过程就是要把逆序的数的次序进行纠正的过程,复杂程度很有可能与逆序对的数量是有关的。
1. 逆序对的数量与数组中重复元素的数量是无关系的,因此,不稳定的排序方法如快排等不在考虑范围内
2. 稳定排序算法中,先看冒泡排序,每次进行冒泡交换数组元素时,相邻元素交换位置,不影响这两个元素与数组其它元素之间的顺序(逆序)关系,也就是说,每冒泡交换一次,数组逆序对数减一。冒泡交换的总次数就等于逆序对数。
3. 选择排序和冒泡排序相似,插入元素位置移动时不影响的是该元素之后的数组的逆序数对数量,而插入元素之前的数组是已经排好序的,逆序对为0,每次插入逆序对数的改变量等于元素原始位置与插入位置之间的偏移。
4.归并排序,逆序对数的改变发生在合并这一步。数组分成两半A=[A1,A2],A1,A2都是排序好的数组,合并时指向A1的指针为i,指向A2的指针为j,如果A2的元素均比A1的元素大,逆序对数量为0,如果A2的元素全比A1的小,那么逆序数为length(A2)*length(A1)。可以修改归并排序,当j++时,判断当前A1中还有多少个元素没有被归并,当前合并的A2元素产生的逆序对数就是多少,累积求和就可以搞定。
修改冒泡和选择排序求逆序数复杂度是O(n^2),利用归并排序求逆序对数的算法复杂度可以降低到O(nlogn)
/*
* nixushu.c
*
* Created on: 2012-9-19
* Author: sangerhoo
*/
#include<iostream> using namespace std; /* 归并求逆序对数, arr 存储最终有序结果 * 在函数外申请一个临时数组作为参数传入 * 避免递归不断创建临时数组的开销 */
int Merge(int * arr, int beg, int mid, int end, int * tmp_arr) {
memcpy(tmp_arr + beg, arr + beg, sizeof(int) * (end - beg + 1));
int i = beg;
int j = mid + 1;
int k = beg;
int inversion = 0; // 合并过程中的逆序数
while (i <= mid && j <= end) {
if (tmp_arr[i] <= tmp_arr[j]) {
arr[k++] = tmp_arr[i++];
} else {
arr[k++] = tmp_arr[j++];
inversion += (mid - i + 1);
}
}
while (i <= mid) {
arr[k++] = tmp_arr[i++];
}
while (j <= end) {
arr[k++] = tmp_arr[j++];
}
return inversion;
}
int MergeInversion(int * arr, int beg, int end, int * tmp_arr) {
int inversions = 0; // 记录倒序数 if(beg < end)
{
int mid = (beg + end) >> 1;
inversions += MergeInversion(arr, beg, mid, tmp_arr);
inversions += MergeInversion(arr, mid + 1, end, tmp_arr);
inversions += Merge(arr, beg, mid, end, tmp_arr);
}
return inversions;
}
/* 测试序列 :answer: 10 */
void main() {
int testPoint[10] = { 1, 4, 2, 9, 48, 15, 13, 44, 6, 90 };
int arrcopy[10]; // 临时数组
memcpy(arrcopy,testPoint,sizeof testPoint);
printf("the num of inversions is: %d\n",
MergeInversion(testPoint, 0, 9, arrcopy));
}
}