#include <stdio.h>
#include <stdlib.h>
#include <time.h>
int prime(int N,unsigned long long int* sum)
{
unsigned int i,j,k,count_p=0;
if(N>100000000)
i=17;
else if(N>70000)
i=10;
else if(N>10000)
i=5;
else if(N>=100)
i=4;
else if(N>=2)
i=1;
else
return 0;
//unsigned int* Prime=(int*)malloc(N/i*(sizeof(unsigned int)));
char* flag=(char*)malloc(N+1);
for(flag[0]=flag[1]=-1,i=2;i<=N;i++)
flag[i]=0; //0代表质数,1代表合数。未筛选前全默认是质数。
for(i=2;i<=N;i++)
{
if(1==flag[i])
continue; //如果i已经是被标记为合数,则退出
else //否则i就是质数,就要把i的倍数都标记为合数
{
*sum+=i;
count_p++;
for(j=i+i;j<=N;j+=i)
flag[j]=1;
}
}
free(flag);//free(Prime);
return count_p;
} //int prime
int main(int argc,char* argv[])
{
int N=1000000000;
unsigned long long int sum=0;
clock_t after,before;
printf("别动,请耐性等待程序执行...\n");
before=clock();
int num=prime(N,&sum);
after=clock();
printf("%d 以内有 %d 个质数。质数之和为%llu。\n计算用时为 %d 毫秒。\n",N,num,sum,after-before);
} //int main
Eraosthenes质数数筛选法
最新推荐文章于 2024-05-30 16:17:47 发布