POJ2299-Ultra-QuickSort

 

转载请注明出处：優YoU http://user.qzone.qq.com/289065406/blog/1304259927

题目大意：

给出长度为n的序列，每次只能交换相邻的两个元素，问至少要交换几次才使得该序列为递增序列。

 

解题思路：

一看就是冒泡，交换一次记录一次就可以了

但是n的范围达到50W，冒泡O(n^2)的复杂度铁定超时（即使有7000ms，其实这是一个陷阱）

直接用快排又不符合题目的要求（相邻元素交换），快排是建立在二分的基础上的，操作次数肯定比在所要求的规则下的交换次数要更少

 

那么该怎么处理？

 

其实这题题目已经给出提示了：Ultra-QuickSort

特殊的快排，能和快排Quicksort相媲美的就是归并排序Mergesort了，O(nlogn)

 

但是用归并排序并不是为了求交换次数，而是为了求序列的 逆序数（学过《线代》的同学应该很熟悉了）

一个乱序序列的 逆序数 = 在只允许相邻两个元素交换的条件下,得到有序序列的交换次数

 

例如例子的

9 1 0 5 4

由于要把它排列为上升序列，上升序列的有序就是  后面的元素比前面的元素大

而对于序列9 1 0 5 4

9后面却有4个比9小的元素，因此9的逆序数为4

1后面只有1个比1小的元素0，因此1的逆序数为1

0后面不存在比他小的元素，因此0的逆序数为0

5后面存在1个比他小的元素4, 因此5的逆序数为1

4是序列的最后元素，逆序数为0

 

因此序列9 1 0 5 4的逆序数 t=4+1+0+1+0 = 6  ，恰恰就是冒泡的交换次数

 

PS：注意保存逆序数的变量t，必须要用__int64定义，int和long都是无法保存的。。。。会导致WA。 以前的long long 在现在的VC编译器已经无法编译了。

注意__int64类型的输出必须使用指定的c格式输出，printf(“%I64d”,t);

cout是无法输出__int64类型的

 

序列数组s[]用int就足够了，每个元素都是小于10E而已

 

 

//Memory Time
//3768K 2422MS 

#include<iostream>
using namespace std;

const int inf=1000000000;  //10E

__int64 t; //s[]数列逆序数

void compute_t(int* s,int top,int mid,int end)
{
	int len_L=mid-top+1;
	int len_R=end-mid;

	int* left=new int[len_L+2];
	int* right=new int[len_R+2];

	int i,j;
	for(i=1;i<=len_L;i++)
		left[i]=s[top+i-1];
	left[len_L+1]=inf;   //设置无穷上界，避免比较大小时越界

	for(j=1;j<=len_R;j++)
		right[j]=s[mid+j];
	right[len_R+1]=inf;   //设置无穷上界，避免比较大小时越界

	i=j=1;
	for(int k=top;k<=end;)
		if(left[i]<=right[j])
			s[k++]=left[i++];
		else
		{
			s[k++]=right[j++];
			t+=len_L-i+1;    //计算逆序数
		}

	delete left;
	delete right;

	return;
}

void mergesort(int* s,int top,int end)
{
	if(top<end)
	{
		int mid=(top+end)/2;
		mergesort(s,top,mid);
		mergesort(s,mid+1,end);
		compute_t(s,top,mid,end);
	}
	return;
}

int main(void)
{
	int n;  //序列长度
	while(cin>>n)
	{
		if(!n)
			break;

		/*Input*/

		int* s=new int[n+1];
		for(int i=1;i<=n;i++)
			cin>>s[i];

		/*Merge-Sort*/

		t=0;  //initial
		mergesort(s,1,n);

		/*Output*/

		printf("%I64d\n",t);

		/*Relax*/

		delete s;

	}
	return 0;
}