素数生成算法

 

 

质数的定义

一个数，如果只有1和它本身两个因数，这样的数叫做质数，又称素数。 

 

试除判断法

算法描述：从上述定义可知，素数不能被1和它本身之外的数整除，所以，判断一个数x是否素数只要看它是否能被2~sqrt(x)间的数整除即可；而求N内所有素数则是循环重复上述过程。

优化：判断一个数x是否素数只要看它是否能被2~sqrt(x)间的素数整除即可。

 

#include <math.h> #include <list> // 试除法求maxNum内的所有素数,返回个数 unsigned int PrimeNum_1(unsigned int maxNum) { std::list<int> primes; primes.push_back(2); for(unsigned int i = 3; i < maxNum; ++i) { int k = (int)sqrt((double)i); bool bPrime = true; // 试除因子选已经筛选出的素数，且不大于k的素数 for(std::list<int>::const_iterator iter = primes.begin(); (*iter) <= k && iter != primes.begin(); ++iter) { if(i%(*iter) == 0) { bPrime = false; break; } } if(bPrime) primes.push_back(i); } return primes.size(); } 

 

 

合数过滤筛选法

算法描述：我们知道，素数N不能被2~(N-1)间的任何数整除；反过来看，只要能被2~(N-1)间的任何数整除的N，都不是素数。所以我们可以采用一个简单的排除法：就是对N以内的所有数，只要逐个去除值为2~(N-1)的倍数的数，剩下的就是素数。

 

// 合数过滤法求maxNum内的所有素数，返回个数 unsigned int PrimeNum_2(unsigned int maxNum) { unsigned int count = 0; // 分配标记数组，浪费flag[0] char* flag = (char*)malloc(maxNum+1); unsigned int i, j; // 初始化标记 flag[2] = 1; // 步长为2，maxNum为奇数时，没有初始化flag[maxNum] for(i = 3; i < maxNum; ) { flag[i++] = 1; // 所有大于2的偶数，均是非素数 flag[i++] = 0; } if(maxNum%2) { flag[maxNum] = 1; } for(i = 3; i < maxNum/2; ++i) { // 已经判定是和数的，不再判定 if(!flag[i]) continue; // 用加法代替乘法，加法运算更快 for(j = i+i; j <= NULL; j += i) { flag[j] = 0; } } // 剩下的凡是flag标识为1的索引，即是输出 std::list<unsigned int> primes; for(i = 2; i <= maxNum; i++) { if(flag[i]) primes.push_back(i); } return primes.size(); } 

 

PS：上面的合数过滤法在初始化时仅仅过滤了偶数。对于maxNum比较大时，考虑一次性过滤更多的呢？

 

// 合数过滤法求maxNum内的所有素数，返回个数 unsigned int PrimeNum_3(unsigned int maxNum) { unsigned int i; std::list<unsigned int> primes; // 忽略flag[0] char* flag = (char*)malloc(maxNum+1); unsigned int fbLen = 2*3*5*7*11*13; // 忽略filterBlock[0] char filterBlock[fbLen+1]; for(i = 1; i <= fbLen; ) { filterBlock[i++] = 1; // 6*i+1 filterBlock[i++] = 0; // 6*i+2，合数 filterBlock[i++] = 0; // 6*i+3，合数 filterBlock[i++] = 0; // 6*i+4，合数 filterBlock[i++] = 1; // 6*i+5 filterBlock[i++] = 0; // 6*i+6，合数 } // 初始化5,7,11,13的倍数 for(i = 5; i <= fbLen; i += 5) filterBlock[i] = 0; for(i = 7; i <= fbLen; i += 7) filterBlock[i] = 0; for(i = 11; i <= fbLen; i += 11) filterBlock[i] = 0; for(i = 13; i <= fbLen; i += 13) filterBlock[i] = 0; // 初始化，过滤所有2,3,5,7,11,13的倍数，并且采用块拷贝的形式 int remainderLen = maxNum%fbLen; char* p = flag+1; char* pstop = p+maxNum-remainderLen; while(p < pstop) { // 块拷贝，要从filterBlock[1]开始，因为忽略了filterBlock[0] memcpy(p, filterBlock+1, fbLen); p += fbLen; } if(remainderLen > 0) { memcpy(p, filterBlock+1, remainderLen); } // 重新置2,3,5,7,11,13号元素的标识为1 flag[2] = flag[3] = flag[5] = flag[7] = flag[11] = flag[13] = 1; // 开始过滤，从17开始，因为2,3,5,7,11,13的倍数已经过滤，剩下的是13以后的素数的倍数且倍数为13以后的素数 // 倍数存在于[17, maxNum/13]中 for(i = 17; i <= maxNum/13; ++i) { if(0 == flag[i]) continue; // 还需要判定不是素数的合数是那些具有大于13的素数因子，且一定是奇数倍 unsigned int step = i+i; for(unsigned int j = i*17; j <= maxNum; j += step) { flag[j] = 0; } } // 剩下的凡是flag标识为1的索引，即是素数 for(i = 2; i <= maxNum; i++) { if(flag[i]) primes.push_back(i); } free(flag); return primes.size(); }