首页 > 代码库 > 素数打表法。

素数打表法。

打表法1【常规打表法】。

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

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");
} 


素数打表法。