埃拉托色尼筛选法求素数

算法的思想很简单，找一张表，先把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;//1061109567
typedef long long ll;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
using 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;//1061109567
typedef long long ll;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
using 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;
}