您的位置:首页 > 其它

素数打表法。

2015-08-13 19:27 176 查看
打表法1【常规打表法】。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<math.h>
#define N 100000
bool p
;
int prim
;

void is_prim()
{
int i,j,n;
for(i=1;i<=N;i++)
p[i]=true;
n=(int)sqrt(N);
for(i=2;i<=n;i++)
{
for(j=i+i;j<=N;j+=i)
{
p[j]=false;
}
}
j=1;
for(i=1;i<=N;i++)
{
if(p[i])
{
prim[j++]=i;
}
}
}

int main()
{
int i;
is_prim();
for(i=1;i<=1000;i++)
printf("%d ",prim[i]);
system("pause");
return 0;
}


打表法2【压缩路径打表法】

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<iostream>
using namespace std;
int n;
bool vis[10000000];
int prim[10000000];

void is_prim()
{
memset(vis,true,sizeof(vis));
int num=0;
for(int i=2;i<=n;i++)
{
if(vis[i]==true)
{
num++;
prim[num]=i;
}
for(int j=1;((j<=num)&&(i*prim[j]<=n));j++)
{
vis[i*prim[j]]=false;
if(i%prim[j]==0)
break;
}
}
}

int main()
{
memset(prim,0,sizeof(prim));
int cnt=0;
scanf("%d",&n);
is_prim();
for(int i=0;i<=n;i++)
if(prim[i])
{
printf("%d ",prim[i]);
cnt++;
}
printf("\n");
system("pause");
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: