素数算法的优化之路

转载

mob60475707d572 2016-03-01 10:57:00

文章标签 质因子筛选法执行时间数组时间复杂度 文章分类 数据结构与算法人工智能

一、素数的定义

质数又称素数。指在一个大于1的自然数中，除了1和此整数自身外，不能被其它自然数(不包含0)整除的数。由于合数是由若干个质数相乘而得来的，所以。没有质数就没有合数，由此可见质数在数论中有着非常重要的地位。

比方：2。3，5。7，9.....都是素数。

二、构造素数算法

写算法之前。先来说说下面这个东西：

对于随意一个合数n，假设它有两个质因子x，y。显然n = x*y。

所以。由不等式性质可得，x <= sqrt(n)，即 x <= n^(1/2)。

推广一下，对于随意一个合数，假设它有k个质因子x1，x2。x3，.... ，xk。显然n = x1 * x2 * x3 * .... * xk。

所以。由不等式性质可得， x <= n^(1/k)，即每一个质因子必小于或者等于n开k次方，当中x 属于集合{ x1，x2。x3，.... ，xk }。

定义一个全局变量和一个全局数组，要是数组是分配在栈上的话，是有容量限制的。所以我们就定义个全局的。

const unsigned int N = 50000;

bool flag[N+1];

1、算法一：依照定义，就可以高速写出一个简单的构造素数算法

void PrimeCreateCommon()//依照定义
{
  for ( unsigned int i = 2; i <= N; ++i )
  {
    //由前面的铺垫，这里不难理解吧。由于我不知道这个数有几个质因子，
    //我就如果它仅仅有两个。范围拉大一些。这个数的全部质因子都小于或
    //等于sqrt( (double)i ) + 1。其他的约数，必定是这些质因子的倍数
    //而已所以就不是必需继续除下去了。
    bool flag = true;
    for ( unsigned int j = 2; j <= sqrt((double)i)+1; ++j )  
    {                              
      if ( i % j == 0 )         
      {                        
        flag = false;
      }
    }
    //if ( flag )
    //{
    //  cout << i << " "; 
    //}
  }
  cout << endl;
}

分析：该算法的时间复杂度为O（n*sqrt(n) )，在n比較大时（比方n = 50000），我的机器跑了个二十来秒。就也叫算法。呵呵，所以接下来我们来看还有一种算法，那就是古希腊数学家的埃拉托色尼给出的一个比較省力的算法------埃拉托色尼筛法。

2、算法二：筛选法

首先，列出从2開始的数。然后，将2记在素数列表上，再划去全部2的倍数。

依据定义，剩下的最小的数——在这里是3——必然是素数。

将这个数记在素数列表上，再划去全部它的倍数，这样又会剩下一些数，取当中最小的，如此重复操作。

最后剩下的都是素数。

void PrimeCreateExt()//筛选法优化
{
  memset( flag, 0, (N+1)*sizeof(bool) );
  for ( unsigned int i = 2; i <= sqrt(double(N)) + 1; ++i )
  {
    if ( !flag[i] )
    {
      for ( unsigned int j = 2*i; j <= N; j += i )
      {
        flag[j] = true;
      }
    }

  }

//  for ( int i = 2; i <= N; ++i )
//  {
//    if ( !flag[i] )
//    {
//      cout << i << " ";
//    }
//  }
//  cout << endl;
}

三、两种算法的执行时间对照

//測试函数
int main()
{
  clock_t start, finish;
  double duration = 0.0;

  cout << "筛选算法产生素数：" << endl;
  start = clock();
  PrimeCreateExt();
  finish = clock();
  duration = (double)(finish-start);
  cout << "时间消耗:" << duration << "ms" << endl << endl;

  cout << "普通算法产生素数：" << endl; 
  start = clock();
  PrimeCreateCommon();
  finish = clock();
  duration = (double)(finish-start);
  cout << "时间消耗:" << duration << "ms" << endl << endl; 
}