实用算法实践-第1篇 排序 .

1.1    选择排序

许多排序算法比选择排序快多了,但是选择排序能够将排序的交换次数降到最少,这是它比冒泡排序的优点之所在。通过下面这个实例就可以深刻认识这一点。

1.1.1   实例

PKU JudgeOnline, 1674, Sorting by Swapping.

1.2    库函数的应用

快速排序函数qsort()非常好用,下面就是一些应用例子。

据说,STL的sort()函数比stdlib.h 的qsort()函数快

1.2.1   库函数排序整数数组

  1. #include <stdlib.h>   
  2. int cmp(const void*p1, const void*p2)  
  3. {  
  4.     return *(int *)p1 - *(int*)p2;  
  5. }  
  6. qsort(&people[1],N, sizeof(people[0]), cmp);  
#include <stdlib.h>
int cmp(const void*p1, const void*p2)
{
    return *(int *)p1 - *(int*)p2;
}
qsort(&people[1],N, sizeof(people[0]), cmp);

上面的cmp函数中如果p1-p2则是非降序,如果p2-p1则是非升序。

1.2.2   库函数排序字符串数组

  1. char word[1002][21];  
  2. int cmp(const void*p1, const void*p2)  
  3. {  
  4.      return strcmp((char*) p1,(char*)p2);  
  5. }  
  6. qsort(&word[0],num, 21 * sizeof(char),cmp);  
char word[1002][21];
int cmp(const void*p1, const void*p2)
{
     return strcmp((char*) p1,(char*)p2);
}
qsort(&word[0],num, 21 * sizeof(char),cmp);

1.2.3   库函数排序结构体数组

1.2.4   实例

PKU JudgeOnline, 3664, Election Time.

1.2.5   问题描述

N个牛的选举分为两个阶段,第一阶段得票最高的前K个进入第二阶段,第二阶段得票最高的获胜。先输入N、K,然后是N对两阶段得票数。问哪只牛获胜。

1.2.6   输入

53

310

92

56

84

65

 

1.2.7   输出

5

1.2.8   程序

  1. #include<stdio.h>   
  2. #include<string.h>   
  3. #include <stdlib.h>   
  4. #include <iostream.h>   
  5. struct Node  
  6. {  
  7.      int data1;  
  8.      int data2;  
  9.      int pos;  
  10. };  
  11. int cmp1(const void*p1, const void*p2)  
  12. {  
  13.     return(*(Node *)p2).data1 > (*(Node *)p1).data1 ? 1 : -1;  
  14. }  
  15. int cmp2(const void*p1, const void*p2)  
  16. {  
  17.     return(*(Node *)p2).data2 > (*(Node *)p1).data2 ? 1 : -1;  
  18. }  
  19. Nodefirst[50001];  
  20. int main()  
  21. {  
  22.      int i;  
  23.      int N, K;  
  24.      while(cin>> N >> K){  
  25.          for(i =0; i < N; i++){  
  26.               scanf("%d%d",&first[i].data1, &first[i].data2);  
  27.               first[i].pos = i + 1;  
  28.          }  
  29.          qsort(first, N,sizeof(first[0]),cmp1);  
  30. /*       for(i = 0; i< N; i++){ 
  31.               cout<<first[i].pos<<" : "<<first[i].data1<<": "<<first[i].data2<<endl; 
  32.          }*/  
  33.          qsort(first, K,sizeof(first[0]),cmp2);  
  34. /*       for(i = 0; i< N; i++){ 
  35.               cout<<first[i].pos<<" : "<<first[i].data1<<": "<<first[i].data2<<endl; 
  36.          }*/  
  37.          cout << first[0].pos<<endl;  
  38.      }  
  39.      return 1;  
  40. }  
#include<stdio.h>
#include<string.h>
#include <stdlib.h>
#include <iostream.h>
struct Node
{
     int data1;
     int data2;
     int pos;
};
int cmp1(const void*p1, const void*p2)
{
    return(*(Node *)p2).data1 > (*(Node *)p1).data1 ? 1 : -1;
}
int cmp2(const void*p1, const void*p2)
{
    return(*(Node *)p2).data2 > (*(Node *)p1).data2 ? 1 : -1;
}
Nodefirst[50001];
int main()
{
     int i;
     int N, K;
     while(cin>> N >> K){
         for(i =0; i < N; i++){
              scanf("%d%d",&first[i].data1, &first[i].data2);
              first[i].pos = i + 1;
         }
         qsort(first, N,sizeof(first[0]),cmp1);
/*       for(i = 0; i< N; i++){
              cout<<first[i].pos<<" : "<<first[i].data1<<": "<<first[i].data2<<endl;
         }*/
         qsort(first, K,sizeof(first[0]),cmp2);
/*       for(i = 0; i< N; i++){
              cout<<first[i].pos<<" : "<<first[i].data1<<": "<<first[i].data2<<endl;
         }*/
         cout << first[0].pos<<endl;
     }
     return 1;
}

1.3    快速排序的随机化

并不是所有的输入数据的所有排列都是等可能的。为了使得算法能获得较好的平均情况性能,可以加入随机化成分。

1.3.1   程序(未测试)

  1. #include<stdio.h>   
  2. #include<string.h>   
  3. #include <stdlib.h>   
  4. #define maxNum 10001   
  5. #define ONLINE_JUDGE 0   
  6. struct cow{  
  7.     int S;  
  8.     int E;  
  9.     //int sequence;   
  10.     //牛的编号,也就是在输入数据的第几个,   
  11.     //因为在输出的时候需要按照编号大小依次输出值   
  12.     //又因为cow不发生变化,所以可以不存   
  13. };  
  14. cow cow[maxNum];  
  15. int pos[maxNum];  
  16. //在排序之后第i大的数是cow数组的第pos[i]个元素   
  17. //在任何时候,cow都不会变化。   
  18. //所有排序都需要通过调整pos数组的指向来进行。   
  19. int cmp(int i, int j)  
  20. {  
  21.     int Si;  
  22.     int Sj;  
  23.     int Ei;  
  24.     int Ej;  
  25.     Ei = cow[pos[i]].E;  
  26.     Ej = cow[pos[j]].E;  
  27.     Si = cow[pos[i]].S;  
  28.     Sj = cow[pos[j]].S;/* 
  29.     if(Si > Sj){ 
  30.         return 1; 
  31.     }else{ 
  32.         return -1; 
  33.     }*/  
  34.     if(Si <= Sj && Ej <= Ei && Ei - Si > Ej - Sj){  
  35.         return 1;  
  36.     }else{  
  37.         return -1;  
  38.     }  
  39. }  
  40. int sequence[maxNum];  
  41. //pos数组的“逆”,即若pos[i] = j,那么sequence[j] = i;  
  42. //通过sequence数组数组,可以知道cow数组的第i个元素于第sequence[i]名  
  43. void exchange(int i, int j)  
  44. {  
  45.     int temp;  
  46.     sequence[pos[i]] = j;  
  47.     sequence[pos[j]] = i;  
  48.     temp = pos[i];  
  49.     pos[i] = pos[j];  
  50.     pos[j] = temp;  
  51. }  
  52. int Partition(int p, int r)  
  53. {  
  54.     int i;  
  55.     int j;  
  56.     int x;  
  57.     x = r;  
  58.     i = p - 1;  
  59.     for(j = p; j < r; j++){  
  60.         if(cmp(x, j) > 0)//(cow[j] <= x)  
  61.         {  
  62.             i++;  
  63.             exchange(i, j);  
  64.         }  
  65.     }  
  66.     exchange(i + 1, r);  
  67.     return i + 1;  
  68. }  
  69. int randPartition(int p, int r)  
  70. {  
  71.     int i;  
  72.     i = (rand() % (r - p)) + p;  
  73.     exchange(i, r);  
  74.     return Partition(p, r);  
  75. }  
  76. void randQsort(int p, int r)  
  77. {  
  78.     int q;  
  79.     if(p < r)  
  80.     {  
  81.         q = randPartition(p, r);  
  82.         randQsort(p, q - 1);  
  83.         randQsort(q + 1, r);  
  84.     }  
  85. }  
  86. int main()  
  87. {  
  88. #ifndef ONLINE_JUDGE        
  89.     FILE *fin;  
  90.   
  91.     fin = freopen( "test.txt""r", stdin);  
  92.     if( !fin )  
  93.     {  
  94.         printf( "reopen in file failed...\n");     
  95.         while(1){}  
  96.         return 0;  
  97.     }  
  98.     freopen( "out.txt""w", stdout);  
  99. #endif   
  100.     int i;  
  101.     int j;  
  102.     int N;  
  103.     int result;  
  104.     while(scanf("%d", &N)){  
  105.         if(N == 0)  
  106.         {  
  107.             break;  
  108.         }  
  109.         for(i = 0; i < N; i++){  
  110.             scanf("%d%d", &cow[i].S, &cow[i].E);  
  111.             sequence[i] = i;  
  112.             pos[i] = i;  
  113.         }  
  114.         randQsort(0, N - 1);  
  115.         for(i = 0; i < N - 1; i++){  
  116.             //printf("%d %d\n", cow[pos[i]].S, cow[pos[i]].E);  
  117.             result = sequence[i];             
  118.             for(j = result + 1; j < N; j++)  
  119.             {  
  120.                 if(cmp(j, result) > 0)  
  121.                 {  
  122.                     break;  
  123.                 }  
  124.             }  
  125.             result = N - j;  
  126.             printf("%d ", result);  
  127.         }  
  128.         printf("0\n");  
  129.     }  
  130. #ifndef ONLINE_JUDGE      
  131. //  fclose(stdout);   
  132.     fclose(stdin);  
  133. #endif   
  134.     return 1;  
  135. }  
#include<stdio.h>
#include<string.h>
#include <stdlib.h>
#define maxNum 10001
#define ONLINE_JUDGE 0
struct cow{
	int S;
	int E;
	//int sequence;
	//牛的编号,也就是在输入数据的第几个,
	//因为在输出的时候需要按照编号大小依次输出值
	//又因为cow不发生变化,所以可以不存
};
cow cow[maxNum];
int pos[maxNum];
//在排序之后第i大的数是cow数组的第pos[i]个元素
//在任何时候,cow都不会变化。
//所有排序都需要通过调整pos数组的指向来进行。
int cmp(int i, int j)
{
	int Si;
	int Sj;
	int Ei;
	int Ej;
	Ei = cow[pos[i]].E;
	Ej = cow[pos[j]].E;
	Si = cow[pos[i]].S;
	Sj = cow[pos[j]].S;/*
	if(Si > Sj){
		return 1;
	}else{
		return -1;
	}*/
	if(Si <= Sj && Ej <= Ei && Ei - Si > Ej - Sj){
		return 1;
	}else{
		return -1;
	}
}
int sequence[maxNum];
//pos数组的“逆”,即若pos[i] = j,那么sequence[j] = i;
//通过sequence数组数组,可以知道cow数组的第i个元素于第sequence[i]名
void exchange(int i, int j)
{
	int temp;
	sequence[pos[i]] = j;
	sequence[pos[j]] = i;
	temp = pos[i];
	pos[i] = pos[j];
	pos[j] = temp;
}
int Partition(int p, int r)
{
	int i;
	int j;
	int x;
	x = r;
	i = p - 1;
	for(j = p; j < r; j++){
		if(cmp(x, j) > 0)//(cow[j] <= x)
		{
			i++;
			exchange(i, j);
		}
	}
	exchange(i + 1, r);
	return i + 1;
}
int randPartition(int p, int r)
{
	int i;
	i = (rand() % (r - p)) + p;
	exchange(i, r);
	return Partition(p, r);
}
void randQsort(int p, int r)
{
	int q;
	if(p < r)
	{
		q = randPartition(p, r);
		randQsort(p, q - 1);
		randQsort(q + 1, r);
	}
}
int main()
{
#ifndef ONLINE_JUDGE     
	FILE *fin;

	fin = freopen( "test.txt", "r", stdin);
	if( !fin )
	{
		printf( "reopen in file failed...\n");   
		while(1){}
		return 0;
	}
	freopen( "out.txt", "w", stdout);
#endif
	int i;
	int j;
	int N;
	int result;
	while(scanf("%d", &N)){
		if(N == 0)
		{
			break;
		}
		for(i = 0; i < N; i++){
			scanf("%d%d", &cow[i].S, &cow[i].E);
			sequence[i] = i;
			pos[i] = i;
		}
		randQsort(0, N - 1);
		for(i = 0; i < N - 1; i++){
			//printf("%d %d\n", cow[pos[i]].S, cow[pos[i]].E);
			result = sequence[i];			
			for(j = result + 1; j < N; j++)
			{
				if(cmp(j, result) > 0)
				{
					break;
				}
			}
			result = N - j;
			printf("%d ", result);
		}
		printf("0\n");
	}
#ifndef ONLINE_JUDGE   
//	fclose(stdout);
	fclose(stdin);
#endif
	return 1;
}

1.4    计数排序

《算法导论》中计数排序的伪代码如下所示:

注意:10行和11行弄得这么复杂的原因是为了维持计数排序的稳定性。计数排序的稳定性对于使用该排序的基数排序的正确性来说非常重要。
计数排序的算法时间复杂度是Θ(k + n)。可以证明基于比较的排序算法下界为Ω(n lgn)。线性时间的排序算法不是基于比较的排序算法。

1.5    基数排序

虽然基数排序的时间复杂度是Θ(n),但是由于Θ记号中隐藏了常数因子的影响,所以比较基数排序和基于比较的排序算法的好坏需要考虑以下几个因素:

1. 算法的实现特性。排序算法通常可以比基数排序更为有效地利用硬件缓存。

2. 输入数据的特性。

3. 内存的速度。如果基数排序利用计数排序作为中间稳定排序,那么它就不是一个原地排序,而很多的Θ(n lgn)时间的比较排序算法可以做到原地排序。故此,如果内存容量比较宝贵,则原地排序算法比较可取。

一个非常好地利用了基数排序的例子是后缀数组的Skew构造算法。

1.6    实例

PKU JudgeOnline, 3664, Election Time.

PKU JudgeOnline, 3404, Bridge over a roughriver

PKU JudgeOnline, 1700, Crossing River.

PKU JudgeOnline, 3637, Shopaholic.

PKU JudgeOnline, 1674, Sorting by Swapping.

PKU JudgeOnline, 2231, Moo Volume.
本文章欢迎转载,请保留原始博客链接http://blog.csdn.net/fsdev/article

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值