iCeSnaker - Program rhapsody

谱写程序的狂想曲

导航

C#算法设计与分析-寻找素数

Posted on 2004-07-01 18:49  iCeSnaker  阅读(1107)  评论(0编辑  收藏  举报
在这篇文章中,我将使用C#编制两个寻找素数的算法,说明算法设计的重要性以及算法的分析。

 

       素数寻找问题由来已久,一直是一些数学家追求的目的。关于素数的定义及性质,我就不在这里多叙了,相信大家都对此了如指掌。素数的寻找思路比较的简单,根据素数的性质(素数应该不能被除了1和它自身的其他数整除)我们可以从最小的素数2开始,一直到比它小1的数为止,用这些数去整除它,如果它能被整除则它必定不是素数,这是判断单个素数的方法(这个算法思想最简单,时间复杂度最大)。对于寻找比某一个给定的整数值小的所有素数也可以采用这种方法,不过我们会发现,采用这种单个判断的方法所耗的时间比较多。比如查找不大于10的素数,我们必须从2开始一个个判断,共需判断9个数,事实上按照我们后面讲述的方法,只需循环2次就可以了。因此,下面的两种方法都将基于删除法来做。

 

       我们来看看删除法的思想:

 

1.  将小于给定整数值n的所有正整数加到一个数组中;

2.  删除能够被一些整数整除的数;

3.  数组中遗留的元素就是最后要得到的素数序列。

 

对于第二步,我们将给出两种方法来实现。我们先来看看算法:

 

算法一:

class prime

     {

         public static int[] PrimeList;

         public  static void FindPrime(int n)

         {

              int[] IntList;

              IntList=new int[n];             

              for (int p=2;p<=n;p++) IntList[p-1]=p;

              for (int p=2;p<Math.Sqrt(n);p++)

              {

                   int j=p+1;

                   while (j<=n)

                   {

                       if ((IntList[j-1]!=0 ) && ((IntList[j-1]% p)==0) ) IntList[j-1]=0;

                       j=j+1;

                   }

              }

              int i=0;

              for (int p=2;p<=n;p++)

              {

                   if (IntList[p-1]!=0) i=i+1;

              }

              PrimeList=new int[i];

              i=0;

              for (int p=2;p<=n;p++)

              {

                   if (IntList[p-1]!=0)

                   {

                       PrimeList[i]=IntList[p-1];

                       i=i+1;

                   }                 

              }

         }

     }

 

这这个算法中,删除的数是那些被2开始直到n的平方根的整数整除的数。这个算法比起前面介绍的单个素数的寻找方法要好,它的循环次数减少了一多半,但是这个算法还不是最理想的:

 

1.              例如,6既能被2整除,也能被3整除,那么当p=2时,6被删掉了一次;当p=3时,6又被删除了一次,虽然按照我们设定的算法规则,这不会导致冲突(通过判断IntList数组元素是否为0,若为0就不必重复删除),但是这会使得算法的效率低下。

 

2.              还有计算素数序列元素个数时,我们也走了弯路。第一步,我们先计算出了数组元素大小,第二步才开始赋值,事实上这两步我们可以减去计算数组大小这一步,可以把它放在前面完成。

 

3.              已经被删除了的元素,也就是那些不是素数的元素,可以不用拿他们去整除整数,例如4不用拿去整除8,因为能被4整除的数肯定能被2整除,已经在前面循环中被删除了。

 

基于上述考虑,我们得到了一个效率更加高的算法:

 

class primegood

     {

         public static int[] PrimeList;

         public static void FindPrime(int n)

         {

              int[] IntList;

              int len=n-1;

              IntList=new int[n];

              for (int p=2;p<=n;p++) IntList[p-1]=p;

              for (int p=2;p<Math.Sqrt(n);p++)

              {

                   if (IntList[p-1]==0) continue;

                   int j=p*p;

                   while (j<=n)

                   {

                       if (IntList[j-1]!=0 )

                       {

                            IntList[j-1]=0;

                            len=len-1;

                       }

                       j=j+p;

                   }

              }

              PrimeList=new int[len];

              int i=0;

              for (int p=2;p<=n;p++)

              {

                   if (IntList[p-1]!=0)

                   {

                       PrimeList[i]=IntList[p-1];

                       i=i+1;

                   }                 

              }

         }

     }

 

这个算法思想和前面的算法完全一样,不过改正了上面算法中不完善的一些内容。

 

为了说明这两个算法的效率区别,我们编制了如下的主程序来比较一下他们的差异:

 

static void   Main()

         {

              Console.WriteLine("Start!");

              DateTime mytime5=DateTime.Now;

              primegood.FindPrime(100000);

              /*for (int i=0;i<=primegood.PrimeList.Length-1;i++)

              {

                   Console.WriteLine(primegood.PrimeList[i]);

              }*/

              DateTime mytime6=DateTime.Now;

              TimeSpan timeadd3=mytime6-mytime5;

              Console.WriteLine(timeadd3.Ticks);

              DateTime mytime1=DateTime.Now;

              prime.FindPrime(100000);

              DateTime mytime2=DateTime.Now;

              TimeSpan timeadd=mytime2-mytime1;

              DateTime mytime3=DateTime.Now;

              primegood.FindPrime(100000);

              DateTime mytime4=DateTime.Now;

              TimeSpan timeadd2=mytime4-mytime3;

              Console.WriteLine(timeadd.Ticks);

              Console.WriteLine(timeadd2.Ticks);

         }

     }

 

通过运行这个程序,可以发现他们的差别是如此的大(前面的算法所耗时间几乎是后面算法的3060倍),事实上,这两个算法的时间复杂度近似为:n1.5);n);可见,对于同一个问题有着多种不同复杂性的算法实现,算法设计是一门十分重要的学问。