数据结构上机实践第二周项目3——体验复杂度

【项目3 - 体验复杂度】

在数据结构与算法中，对于算法的选择，要考虑到时间复杂度的重要性，在小规模的计算中，或许不同时间复杂度的程序所用时间并无多少影响，但在实际应用中，大数据时代，我们会明白复杂度不同对于算法的差异，下面将用实际例子体验复杂度。

排序是计算机科学中的一个基本问题，产生了很多种适合不同情况下适用的算法，也一直作为算法研究的热点。本项目提供两种排序算法，复杂度为O(n2)的选择排序selectsort，和复杂度为O(nlogn)的快速排序quicksort，在main函数中加入了对运行时间的统计。 

我们将以近十万条数据作为输入测试体验。

选择排序的源程序  （复杂度是O(n2)） 

[cpp]  view plain  copy

1. #include <stdio.h>  
2. #include <time.h>  
3. #include <stdlib.h>  
4. #define MAXNUM 100000  
5. void selectsort(int a[], int n)  
6. {  
7.         int i, j, k, tmp;  
8.         for(i = 0; i < n-1; i++)  
9.         {  
10.                 k = i;  
11.                 for(j = i+1; j < n; j++)  
12.                 {  
13.                         if(a[j] < a[k])  
14.                                 k = j;  
15.                 }  
16.                 if(k != j)  
17.                 {  
18.                         tmp = a[i];  
19.                         a[i] = a[k];  
20.                         a[k] = tmp;  
21.                 }  
22.         }  
23. }  
24.   
25. int main()  
26. {  
27.     int x[MAXNUM];  
28.     int n = 0;  
29.     double t1,t2;  
30.     FILE *fp;  
31.     fp = fopen("numbers.txt", "r");  
32.     if(fp==NULL)  
33.     {  
34.         printf("打开文件错！请下载文件，并将之复制到与源程序文件同一文件夹下。\n");  
35.         exit(1);  
36.     }  
37.     while(fscanf(fp, "%d", &x[n])!=EOF)  
38.         n++;  
39.     printf("数据量：%d, 开始排序....", n);  
40.     t1=time(0);  
41.     selectsort(x, n);  
42.     t2=time(0);  
43.     printf("用时 %d 秒!", (int)(t2-t1));  
44.     fclose(fp);  
45.     return 0;  
46. }  

快速排序源代码 （ 复杂度为 O(nlogn)）     

[cpp]  view plain  copy

1. #include <stdio.h>  
2. #include <time.h>  
3. #include <stdlib.h>  
4. #define MAXNUM 100000  
5. void quicksort(int data[],int first,int last)  
6. {  
7.     int i, j, t, base;  
8.     if (first>last)  
9.         return;  
10.     base=data[first];   
11.     i=first;  
12.     j=last;  
13.     while(i!=j)   
14.     {  
15.         while(data[j]>=base && i<j)   
16.             j--;  
17.         while(data[i]<=base && i<j)   
18.             i++;  
19.         /*交换两个数*/  
20.         if(i<j)  
21.         {  
22.             t=data[i];  
23.             data[i]=data[j];  
24.             data[j]=t;  
25.         }  
26.     }  
27.     data[first]=data[i];   
28.     data[i]=base;   
29.     quicksort(data,first,i-1);  
30.     quicksort(data,i+1,last);   
31. }  
32.   
33. int main()  
34. {  
35.     int x[MAXNUM];  
36.     int n = 0;  
37.     double t1,t2;  
38.     FILE *fp;  
39.     fp = fopen("numbers.txt", "r");  
40.     if(fp==NULL)  
41.     {  
42.         printf("打开文件错！请下载文件，并将之复制到与源程序文件同一文件夹下。\n");  
43.         exit(1);  
44.     }      
45.     while(fscanf(fp, "%d", &x[n])!=EOF)  
46.         n++;  
47.     printf("数据量：%d, 开始排序....", n);  
48.     t1=time(0);  
49.     quicksort(x, 0, n-1);  
50.     t2=time(0);  
51.     printf("用时 %d 秒!", (int)(t2-t1));  
52.     fclose(fp);  
53.     return 0;  
54. }  

第一种排序的运行结果如下：

用时达到12秒。

第二种排序方法运行结果如下：

同样的数据，第二种排序方法竟然在瞬间完成。

这样的结果让我们切实的体验到了复杂度在实际应用中的差异。

（2）汉诺塔 
　　有一个印度的古老传说：在世界中心贝拿勒斯（在印度北部）的圣庙里，一块黄铜板上插着三根宝石针。印度教的主神梵天在创造世界的时候，在其中一根针上从下到上地穿好了由大到小的64片金片，这就是所谓的汉诺塔。不论白天黑夜，总有一个僧侣在按照下面的法则移动这些金片：一次只移动一片，不管在哪根针上，小片必须在大片上面。僧侣们预言，当所有的金片都从梵天穿好的那根针上移到另外一根针上时，世界就将在一声霹雳中消灭，而梵塔、庙宇和众生也都将同归于尽。 
　　可以算法出，当盘子数为n个时，需要移动的次数是f(n)=2n−1。n=64时，假如每秒钟移一次，共需要18446744073709551615秒。一个平年365天有31536000秒，闰年366天有31622400秒，平均每年31556952秒，移完这些金片需要5845.54亿年以上，而地球存在至今不过45亿年，太阳系的预期寿命据说也就是数百亿年。真的过了5845.54亿年，不说太阳系和银河系，至少地球上的一切生命，连同梵塔、庙宇等，都早已经灰飞烟灭。据此，2n从数量级上看大得不得了。

用递归算法求解汉诺塔问题，其复杂度可以求得为O(2n)，是指数级的算法。体验盘子数discCount为4、8、16、20、24时在时间耗费上的差异。

源代码如下：  

[cpp]  view plain  copy

1. #include <stdio.h>  
2. #define discCount 4  
3. long move(int, char, char,char);  
4. int main()  
5. {  
6.     long count;  
7.     count=move(discCount,'A','B','C');  
8.     printf("%d个盘子需要移动%ld次\n", discCount, count);  
9.     return 0;  
10. }  
11.   
12. long move(int n, char A, char B,char C)  
13. {  
14.     long c1,c2;  
15.     if(n==1)  
16.         return 1;  
17.     else  
18.     {  
19.         c1=move(n-1,A,C,B);  
20.         c2=move(n-1,B,A,C);  
21.         return c1+c2+1;  
22.     }  
23. }  

当discount为4时，运行结果如下：

当discount为8时，运行结果如下：

似乎运行时间并没有多大差异

当discount为16时，运行结果如下：

当discount为20时，运行结果如下：

当discount为24时，运行结果如下：

这些次数，望而生畏，连程序计算都稍微卡顿了一下

当discount为26时

运行已经明显有将近2秒的延迟了

当discount为30时

程序整整运行了14秒之久，可见复杂度在实际应用中的差异。