素数的检测

常用的素数检测算法

（一）Eratosthenes 筛法

        Eratosthenes 筛法，也称为埃拉托斯特尼筛法（Eratosthenes Sieve），是一个古老的算法，用于找到小于或等于给定整数的所有质数。这个算法以古希腊数学家埃拉托斯特尼（Eratosthenes of Cyrene）的名字命名，他是第一个使用这种方法来找到质数的人。

算法原理：

        Eratosthenes 筛法的基本思想是从2开始，依次筛选掉2的倍数、3的倍数、4的倍数，以此类推，直到筛选到不超过所给整数的平方根为止。剩下的未被筛选掉的数就是质数。

代码：

#include <stdio.h>
#define TRUE 1
#define FALSE 0
int main()
{
    int SIZE;
    printf("Please input SIZE: ");
    scanf("%d",&SIZE);
    int i,j,k;
    int a[SIZE];
    int* p;
    for(p=a; p<a+SIZE; ++p)
    {
        *p=TRUE;
    }
    a[0]=a[1]=FALSE;
    i=2;
    while(i<SIZE)
    {
        while(a[i++]==TRUE)
        {
            j=i-1;
            break;
        }
        for(k=2; (j*k<SIZE)&&(i<SIZE); ++k)
        {
            a[j*k]=FALSE;
        }
    }
    for(p=a; p<a+SIZE; ++p)
    {
        if(*p==TRUE)
        {
            printf("%8d",p-a);
        }
    }
    printf("\n");
    getchar();
    return 0;
}

（二）Miller-Rabin算法

        Miller-Rabin算法是一种概率性的素数测试算法，用于确定一个给定的数是否是质数。尽管它不像确定性的素数测试算法那样绝对确定一个数是否为质数，但它在实际应用中被广泛使用，因为它具有较高的效率和可扩展性。

算法原理

1. 选择一个要测试的整数n：首先，选择一个要测试的整数n，它应该是一个大于2的奇数。如果n是偶数，那么它肯定不是质数。

2. 将n-1表示为2^r * d：对于奇数n，可以将n-1表示为2的幂r乘以一个奇数d，即n-1 = 2^r * d，其中r和d都是正整数。这是因为n-1是偶数，可以将其分解为2的幂和一个奇数。

3. 选择随机整数a：随机选择一个介于2和n-2之间的整数a。

4. 计算a^d mod n：计算a的d次方模n的结果。可以使用模幂运算（快速幂算法）来高效地计算。

5. 检查一些条件：

   • 如果a^d mod n等于1，那么n可能是质数，但我们需要进行更多的测试。
   • 如果存在一个正整数i（i从0到r-1），使得a^(2^i * d) mod n等于n-1，那么n可能是质数，但我们需要进行更多的测试。
   • 如果以上两个条件都不满足，那么n肯定不是质数。

6. 重复测试：重复步骤3至步骤5，通常重复k次，其中k是一个参数，控制了测试的精度。较大的k值会增加测试的准确性，但也会增加计算时间。

7. 结论：如果n通过了所有测试，那么它被认为是质数的概率非常高。如果n没有通过所有测试，那么它肯定不是质数。

代码：

#include <stdio.h>
#include <stdlib.h>
#include <stdbool.h>
#include <math.h>

// 定义快速幂取模函数
unsigned long long mod_exp(unsigned long long base, unsigned long long exp, unsigned long long mod) {
    unsigned long long result = 1;
    base %= mod;
    while (exp > 0) {
        if (exp % 2 == 1) {
            result = (result * base) % mod;
        }
        base = (base * base) % mod;
        exp /= 2;
    }
    return result;
}

// Miller-Rabin算法主函数
bool is_prime_miller_rabin(unsigned long long n, int k) {
    if (n <= 1 || n == 4) {
        return false;
    }
    if (n <= 3) {
        return true;
    }

    // 寻找d和r，使得n-1 = 2^r * d
    unsigned long long d = n - 1;
    int r = 0;
    while (d % 2 == 0) {
        r++;
        d /= 2;
    }

    // 进行Miller-Rabin测试
    for (int i = 0; i < k; i++) {
        unsigned long long a = 2 + rand() % (n - 3); // 生成2到n-2之间的随机整数
        unsigned long long x = mod_exp(a, d, n);
        if (x == 1 || x == n - 1) {
            continue;
        }
        for (int j = 0; j < r - 1; j++) {
            x = mod_exp(x, 2, n);
            if (x == n - 1) {
                break;
            }
        }
        if (x != n - 1) {
            return false;
        }
    }
    return true;
}

int main() {
    unsigned long long n;
    int k = 10; // 进行Miller-Rabin测试的次数

    printf("输入一个正整数 n：");
    scanf("%llu", &n);

    if (is_prime_miller_rabin(n, k)) {
        printf("%llu 是素数。\n", n);
    } else {
        printf("%llu 不是素数。\n", n);
    }

    return 0;
}