转自:http://www.slyar.com/blog/poj-2299-c.html
Slyar:题目本质就是求逆序对了,简单介绍一下。逆序对是指在序列{a0,a1,a2...an}中,若ai<aj(i>j),则(ai,aj)上一对逆序对。而逆序数顾名思义就是序列中逆序对的个数。例如: 1 2 3是顺序,则逆序数是0;1 3 2中(2,3)满足逆序对的条件,所以逆序数只有1; 3 2 1中(1,2)(1,3)(2,3)满足逆序对,所以逆序是3。由定义不能想象,序列n的逆序数范围在[0,n*(n-1)/2],其中顺序时逆序数为 0,完全逆序时逆序数是n*(n-1)/2。
可以利用归并排序时计算逆序个数,时间复杂度是nlog2n,而空间复杂度2n。 利用归并求逆序是指在对子序列s1和s2在归并时,若s1[i]>s2[j](逆序状况),则逆序数加上s1.length-i,因为s1中i后面的数字对于s2[j]都是逆序的。具体看注释吧。
PS.归并算法基本思路
设两个有序的子文件放在同一向量中相邻的位置上:R[low..m],R[m+1..high],先将它们合并到一个局部的暂存向量R1中,待合并完成后将R1复制回R[low..high]中。
合并过程中,设置i,j和p三个指针,其初值分别指向这三个记录区的起始位置。合并时依次比较R[i]和R[j]的关键字,取关键字较小的记录复制到R1[p]中,然后将被复制记录的指针i或j加1,以及指向复制位置的指针p加1。
重复这一过程直至两个输入的子文件有一个已全部复制完毕(不妨称其为空),此时将另一非空的子文件中剩余记录依次复制到R1中即可。
AC代码:
1 #include <stdio.h> 2 #include <stdlib.h> 3 4 #define MAX 500001 5 6 int n, a[MAX], t[MAX]; 7 __int64 sum; 8 9 /* 归并 */ 10 void Merge(int l, int m, int r) 11 { 12 /* p指向输出区间 */ 13 int p = 0; 14 /* i、j指向2个输入区间 */ 15 int i = l, j = m + 1; 16 /* 2个输入区间都不为空时 */ 17 while(i <= m && j <= r) 18 { 19 /* 取关键字小的记录转移至输出区间 */ 20 if (a[i] > a[j]) 21 { 22 t[p++] = a[j++]; 23 /* a[i]后面的数字对于a[j]都是逆序的 */ 24 sum += m - i + 1; 25 } 26 else 27 { 28 t[p++] = a[i++]; 29 } 30 } 31 /* 将非空的输入区间转移至输出区间 */ 32 while(i <= m) t[p++] = a[i++]; 33 while(j <= r) t[p++] = a[j++]; 34 /* 归并完成后将结果复制到原输入数组 */ 35 for (i = 0; i < p; i++) 36 { 37 a[l + i] = t[i]; 38 } 39 } 40 41 /* 归并排序 */ 42 void MergeSort(int l, int r) 43 { 44 int m; 45 if (l < r) 46 { 47 /* 将长度为n的输入序列分成两个长度为n/2的子序列 */ 48 m = (l + r) / 2; 49 /* 对两个子序列分别进行归并排序 */ 50 MergeSort(l, m); 51 MergeSort(m + 1, r); 52 /* 将2个排好的子序列合并成最终有序序列 */ 53 Merge(l, m, r); 54 } 55 } 56 57 int main() 58 { 59 int i; 60 while(1) 61 { 62 scanf("%d", &n); 63 if (n == 0) break; 64 sum=0; 65 for(i = 0; i < n; i++) 66 { 67 scanf("%d", &a[i]); 68 } 69 MergeSort(0, n - 1); 70 printf("%I64d\n", sum); 71 } 72 //system("pause"); 73 return 0; 74 }