埃拉托色尼筛选法求素数

来源：互联网发布：夏普zd140wj网络设置编辑：程序博客网时间：2024/05/22 19:26

算法的思想很简单，找一张表，先把1划掉，然后在把2的倍数划掉，接着找下一个没有被划掉的数，是3，接着把3的倍数划掉，接着找下一个没有被划掉的数，是5，接着把5的倍数划掉，接着找下一个没有被划掉的数，依次类推，划掉的数不是素数，没有划掉的数是素数

#include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <queue>#include <stack>#include <map>#include <cstring>#include <climits>#include <cmath>#include <cctype>const int inf = 0x3f3f3f3f;//1061109567typedef long long ll;#define lson l,m,rt<<1#define rson m+1,r,rt<<1|1using namespace std;int prime[1000];int isprime[10000];int k=0;void doprime(int n){    for(int i=2; i<=n; i++)    {        if(isprime[i] == 0)        {            prime[k++] = i;            for(int j=i*i; j<=n; j+=i)//2到i-1之前和i的乘积已经被划出来了，所以直接从i*i开始就行，很多题都是在这个模板上搞，只是第二个循环有所不同            {                isprime[j] = 1;            }        }    }}int main(){    memset(isprime,0,sizeof(isprime));    doprime(1000);    int n;    while(scanf("%d",&n) != EOF)    {        if(isprime[n] == 0)        {            printf("yes\n");        }        else        {            printf("no\n");        }    }    return 0;}

上面这个模板不仅能判断1000以内的一个数是否是素数，prime数组里存的是1000以内的所有素数

如果只要求1000内的某个数是否是素数，下面这个就比较省时间了

#include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <queue>#include <stack>#include <map>#include <cstring>#include <climits>#include <cmath>#include <cctype>const int inf = 0x3f3f3f3f;//1061109567typedef long long ll;#define lson l,m,rt<<1#define rson m+1,r,rt<<1|1using namespace std;int isprime[10000];void doprime(int n){    for(int i=2; i*i<=n; i++)//注意这里是小于等于,举个例子啊n是25的时候，没有等于就会造成错误    {        if(isprime[i] == 0)        {            for(int j=i*i; j<=n; j+=i)            {                isprime[j] = 1;            }        }    }}int main(){    memset(isprime,0,sizeof(isprime));    doprime(1000);    int n;    while(scanf("%d",&n) != EOF)    {        if(isprime[n] == 0)        {            printf("yes\n");        }        else        {            printf("no\n");        }    }    return 0;}

1 0