求逆序对个数（inversions）

逆序对：设a是一个有n个不同元素的数组，，如果在i<j的情况下，a[i]>a[j]，则(i, j)称为a的一个逆序对。

给定一个数组，求逆序对的个数，可以通过做循环，在n^2的时间内求出逆序对数目，但通常用归并排序在nlog(n)的时间内解决，只是当求出逆序对数目后，数组也已经排好序了，如果不想让原数组发生变化，可以把原数组拷贝到另一个数组，对新数组操作就好。整个过程一边排序一边数，递归到最深的地方如果有一个数，那么返回0，如果有两个数，看做有两个元素的子数组，两个子数组都已经排好序，并且左边只有一个元素的子数组的逆序对是0，右边只有一个元素的子数组的逆序对是0，那么两个元素的子数组的逆序对的个数就是1,若这两个数逆序，否则是0对。所以，当递归进行到最后一步，也就是这样的：a[1]到a[n/2]已经排好序，并且逆序对的个数为Inv1，a[n/2+1]到a[n]已经排好序，并且逆序对的个数为Inv2，那么整个数组的逆序对总数就是Inv1+Inv2+最后一步归并过程中的逆序对数。那么在merge这一步的逆序对个数又是怎么计算的呢？若i是左边子数组的下标（i的范围是1到mid），j是右边子数组的下标（j的范围是mid+1到n），假设某对(i, j)，有a[i]>a[j]，那么a[i]到a[mid]都大于a[j]，贡献的逆序对数为mid-i+1。

#include<stdio.h> #include<stdlib.h> int merge(int *a, int low, int mid, int high) { int n1 = mid-low+1; int n2 = high-mid; int *tmp1 = malloc(n1*sizeof(int)); int *tmp2 = malloc(n2*sizeof(int)); int i, j, k; for(i=0; i<n1; i++) tmp1[i] = a[low+i]; for(i=0; i<n2; i++) tmp2[i] = a[mid+i+1]; k = low; int cross_inversion = 0; for(i=0, j=0; i<n1 && j<n2;) { if(tmp1[i]<=tmp2[j]) a[k++] = tmp1[i++]; else { //从i到n1-1，所以是n1-1-i+1 = n1-i个逆序对 a[k++] = tmp2[j++]; cross_inversion += (n1-i); } } while(i<n1) a[k++] = tmp1[i++]; while(j<n2) a[k++] = tmp2[j++]; free(tmp1); free(tmp2); return cross_inversion; } int count_inversions(int *a, int low, int high) { int inversions = 0; if(low < high) { int mid = (low+high)/2; inversions += count_inversions(a, low, mid); inversions += count_inversions(a, mid+1, high); inversions += merge(a, low, mid, high); } return inversions; } int main(void) { int a[] = {8, 7, 6, 5, 4, 3, 2, 1}; int n = sizeof(a)/sizeof(a[0]); printf("inversions:%d/n", count_inversions(a, 0, n-1)); return 0; }