素数判定 - 素数打表

最新推荐文章于 2022-04-30 19:47:16 发布

有情剑客无情剑

最新推荐文章于 2022-04-30 19:47:16 发布

阅读量523

点赞数

分类专栏：数据结构与算法文章标签： CC++ 数据结构与算法

本文链接：https://blog.csdn.net/zngto/article/details/20164443

版权

数据结构与算法专栏收录该内容

1 篇文章 0 订阅

订阅专栏

素数：只能被1和其本身整除

普通求一个数是不是素数，可能就是：方法1

bool JuNum(int t_v)
{
    for (int i = 2; i <= sqrt(t_v + 0.5); (i==2)?(i++):(i+=2))
    {
        if (!(t_v%i)) return false;
    }
    return true;
}

然而方法1效率太低，当数字很大时就太过耗时。

下面介绍一种更高效的方法，一种素数打表法，它判断并统计小于一定值内的素数时可以说是线性的：方法2

#define N 100000
#define M 10000
bool isPrime[N+5];
int prime[M];
int cnt;
void JuPr()
{
    int i, j;

    memset(isPrime, 0, sizeof(isPrime));
    cnt = 0;
    for (i = 2; i <= N; i++)
    {
        if (!isPrime[i]) prime[cnt++] = i;
        for (j = 0; j < cnt && (i*prime[j] <= N); j++)
        {
            isPrime[i*prime[j]] = 1;
            //非素数都可以分解成两个数的乘积的形式或者若干个素数的乘积
            //假如非素数n，那么n=A*B（A是素数，B是任意整数）成立
            //那么有1，B=A^K（K是某一整数）；2，A和B互质
            //如果i%prime[j] == 0，那么i = prime[j] * J，且J>prime[j]
            //所以一定存在i*prime[j+1] = I(>i) * prime[k]
            if (!(i%prime[j])) break; //避免对一个数重复多次判定
        }
    }
}