暑假开始啃算法导论了,先解决堆的问题。首先堆排序使用的是最大堆,特点是堆顶元素大于其他所有数,堆排序就是利用了这一点,需要三个函数:
1.维护最大堆性质:MaxHeapify(int a[], int i),基本思路是比较节点i和它的左右孩子,选出最大的交换位置,注意交换后可能导致它的子树不符合最大堆原则,所以要递归再求一遍,它的复杂度为O(lgn)
2.建堆:BuildMaxHeap(int a[]),基本思路是:叶子节点不用进行维护,所以只要从i=heapsize/2开始到i遍历维护一遍堆即可得到一个最大堆。时间复杂度为O(n)
3.堆排序算法:Heapsort(int a[]),基本思路是,先用建堆函数建立一个最大堆,这样最大值就是a[1],我们将a[i]与a[n]交换位置,并将堆内元素减一,之后再对堆顶元素也就是现在的a[1]进行维护最大堆性质操作,就可以确保当前的堆顶是最大值,这样不断交换,直到堆中只用2个元素,总的时间复杂度为O(nlgn)。
堆排序具有归并排序的时间复杂度和插入排序的常数级空间复杂度,是非常优越的排序算法。
#include<stdio.h>
#define maxsize 500
int heapsize;
void exchange(int *a, int *b)
{
int t;
t = *a;
*a = *b;
*b = t;
}
void MaxHeapify(int a[], int i)//在节点和该节点的左右孩子中选出最大的
{
int left = 2 * i;
int right = 2 * i + 1;
int largest;
if (left<=heapsize&&a[left]>a[i])
largest = left;
else
largest = i;
if (right <= heapsize&&a[right] > a[largest])
largest = right;
if (largest != i)
{
exchange(&a[i], &a[largest]);
MaxHeapify(a, largest);//更换后子节点可能违背最大堆的性质,所以再进行一次
}
}
void BuildMaxHeap(int a[])//在线性时间内将数组构造为一个最大堆
{
for (int i = heapsize / 2; i >= 1; i--)
MaxHeapify(a, i);
}
void Heapsort(int a[])//不断进行堆顶的更改,确保堆顶一定是当前最大值
{
BuildMaxHeap(a);
for (int i = heapsize; i >= 2; i--)
{
exchange(&a[i], &a[1]);//将堆顶放到后面
heapsize--;//堆慢慢减少
MaxHeapify(a, 1);
}
}
int main()
{
int a[maxsize],n,i;
while (scanf("%d", &n) != EOF)
{
for (i = 1; i <= n; i++)
scanf("%d", &a[i]);
heapsize = n;
Heapsort(a);
for (i = 1; i <= n; i++)
printf("%d ", a[i]);
printf("\n");
}
return 0;
}