poj 2299（利用合并排序求逆序数）

题目描述：

已知：给定一任意互异整数序列S，只能对S中的相邻元素交换以达到排序的目的

求：将S排序为升序所需最少交换次数

输入：第1行为序列长度N，第2至第N+1行为各元素。有多个case，以N=0为输入结束

输出：各case对应序列排为升序所需最少交换次数

Sample Input：

5

9

1

5

0

4

3

1

2

3

0

Sample Output：

6

0

思路：

先看逆序数的定义：在一个序列中，如果一对元素的前后位置与大小顺序相反，则称它们为一个逆序对。整个序列中的逆序对个数称为该序列的逆序数

本题所求的最少交换次数即输入序列的逆序数，因为每交换一次，序列的逆序数就减一，交换过程结束后，序列的逆序数为0，即排序完毕

另一方面，采用合并排序求逆序数的原因在于：本题数据量很大，若采用时间复杂度为O(n^2)的简单两层循环遍历会导致超时；而采用合并排序，时间复杂度仅为O(n*logn)。不仅如此，合并排序本身的算法思想也很适合求逆序数。在合并排序的”合并“过程中，两个输入序列S1 = (a[low], ..., a[mid])和S2 = (a[mid+1], ..., a[high])是各自有序的，而空间上S2整体在S1之后，因此，若当前插入的是S2中的元素而S1中仍剩余k个元素未插入，则逆序数加k，整个合并排序结束后自然得到了逆序数

代码：

#include <stdio.h>
#include <stdlib.h>
#include <limits.h>
long a[500100];
long long cnt;  //统计逆序数的cnt最大可能为(n-1)!（n为输入序列长度），因此设为long long型；另外，cnt设为全局变量免去了递归过程中的很多麻烦
void merge(long low, long mid, long high){
	long len1 = mid-low+1;
	long len2 = high-mid;
	long * tmp1 = (long *)malloc((len1+1)*sizeof(long));
	long * tmp2 = (long *)malloc((len2+1)*sizeof(long));
	for(long i = 0; i < len1; i++)
		tmp1[i] = a[low+i];
	for(long i = 0; i < len2; i++)
		tmp2[i] = a[mid+1+i];
	tmp1[len1] = LONG_MAX;
	tmp2[len2] = LONG_MAX;
	long p1 = 0;
	long p2 = 0;
	for(long i = 0; i < len1+len2; i++){
		if(tmp1[p1] <= tmp2[p2]){
			a[low+i] = tmp1[p1];
			p1++;
		}
		else{
			a[low+i] = tmp2[p2];
			p2++;
			cnt += len1-p1;  //统计逆序数
		}
	}
	free(tmp1);
	free(tmp2);
}

void merge_sort(long low, long high){
	if(low < high){
		long mid = (low+high)/2;
		merge_sort(low, mid);
		merge_sort(mid+1, high);
		merge(low, mid, high);
	}
}

int main(){
	long n;
	while(1){
		scanf("%ld", &n);
		if(n == 0) break;
		for(long i = 0; i < n; i++)
			scanf("%ld", &a[i]);
		cnt = 0;
		merge_sort(0, n-1);
		printf("%lld\n", cnt);
	}
	return 0;
}