#include<string.h>
#include<algorithm>
#include<math.h>
#include<iostream>
using namespace std;
const int N = 10000001;
bool PRIME[N]={0};//素数表
int PrimeTable[664579 + 10];//这里的664579是1 - 10000000里面的素数的最大个数
//int MAX = (int )sqrt((double)N);
//设置素数表,不管偶数,因为所有的偶数都不是素数,只检查奇数
void getPrimetable()
{
for(int i=3;i<=3333;i+=2)
{
if (PRIME[i]==0)
{
for(int j=i*i;j<N;j+=2*i)//检查9,15(3*5),21(3*7)。。。。奇数乘以偶数得到一偶数,所以这里j+=2*i
{
PRIME[j]=1;//不和要求,不是素数
}
}
}
cout<<"checking prime number..."<<endl;
//将所有的素数都挑拣出来放到primeTable【】中
int total_prime=1;//考虑2是素数
PrimeTable[1]=2;
for (int i=3;i<N;i++)
{
if (i%2 != 0 && PRIME[i]==0)//这里排除掉所有的偶数
{
PrimeTable[++total_prime]=i;
}
}
cout<<"Set prime Table DONE!"<<endl;
}
int main()
{
getPrimetable();
cout<<"input a number between 2 and 10000000"<<endl;
int n,i;
while (cin>>n)
{
if (n<2 || n>1000000
求一千万以内的素数的比较高效的方法,
最新推荐文章于 2021-03-20 14:09:51 发布