素数判定模板

橙子教师

于 2021-08-20 15:44:59 发布

阅读量309

点赞数

分类专栏：质数（素数）问题

本文链接：https://blog.csdn.net/lvcheng0309/article/details/119824745

版权

质数（素数）问题专栏收录该内容

1 篇文章 2 订阅

订阅专栏

一、素数与合数

素数和合数如果一个整数p只有1和p两个因子，则p为素数，亦称为质数不为素数的其它数为合数。如果n为合数，则n必有一个小于或等于n的平方根的数因子。

二、素数的判别方法

1、试除法

int isprime(int x)
{
    int i,k;
    k=(int)sqrt(x);
    for(i=2;i<=k;i++)
    {
         if(x%i==0)
             return 0;
    }
    return 1;
}

或：

int isprime(int x)
{
    int i;
    for(i=2;i*i<=x;i++)
    {
         if(x%i==0)
             return 0;
    }
    return 1;
}

2、高效判别法

素数出现的规律：

当n≧5时，如果n为素数，那么n mod 6 = 1 或 n mod 6 = 5，即：n一定出现在6x（x≥1）两侧。

当x≥1时，有如下表示方法： ……, 6x, 6x+1, 6x+2, 6x+3, 6x+4, 6x+5, 6(x+1), 6(x+1)+1,…… 不在6x两侧的数为6x+2, 6x+3, 6x+4，即2(3x+1), 3(2x+1), 2(3x+2)，它们一定不是素数，所以素数一定出现在6x的两侧。

int isprime(int x)
{
    int i,k;
    k=(int)sqrt(x);
    if(x<2)
        return 0;
    if(x==2 || x==3)
        return 1;
    if(x%6!=1 && x%6!=5)
         return 0;
    for(i=5;i<=k;i+=6)
         if(x%i==0 || x%(i+2)==0)
             return 0;
    return 1;
}

3、埃氏筛选法

埃拉托斯特尼筛选法：开始假设所有范围内的数都是素数，然后将所有这些素数的倍数（肯定不是素数）筛掉，经过数轮筛选，剩下来的就是素数。

20以内素数的埃氏筛选如下：

幸免于难的数就是素数，有：2 3 5 7 11 13 17 19。

#define MAX 1000010
int table[MAX]={1,1};
void isprime(int x)
{
    int i,j; 
    for(i=2;i<MAX;i++)
    {
         if(table[i]==0)
             for(j=i+i; j<MAX; j+=i)
                 table[j]=1;
    }
}

4、欧拉筛选法

欧拉筛法的基本思想：在埃氏筛法的基础上，让每个合数只被它的最小质因子筛选一次，以达到不重复的目的。即：把 prime 里面记录的素数，升序来当做要消去合数的最小素因子。

#define MAXN 100000010
int pri[MAXN];    //存放素数的数组 
int visit[MAXN];  //访问数组，将素数x的倍数做标记
void getPrime()	
{
    int i,j;
    for(i=2;i<=MAXN;i++)
    {
        if (!visit[i])
            pri[++pri[0]] = i;    //记录素数，pri[0]用于计数
        for(j=1; j<=pri[0] && i*pri[j]<=MAXN; j++)
        {
            visit[i*pri[j]] = 1;
            if(i%pri[j] == 0)
                break;
        }
    }
}