Miller-Robin 素数测试法 模板

测试单个素数，出错概率比计算机本身出错的概率还要低

算法是基于费马小定理（format），二次探测定理（x*x % p == 1 ，若P为素数，则x的解只能是x = 1或者x = p - 1）加上迭代乘法判断的Miller算法共同构成的

#include <stdio.h>
#include <string.h>
#include <time.h>
#include <iostream>
#include <string>
 
using namespace std;
 
int N;
 
int witness(int a, int n)//随机生成的a,来检测n的素性 
{
    int ans = 1;
    int t = n - 1;//这里需要注意,你如果没有改变乘方的次数的话,最后的判断就是(ans == a) ? 0 : 1;
                 // 并且还要另外开辟空间来存储开始的a,比较麻烦,所以就这样了; 
    int x;
    while (t)
    {
        if (t & 1)
        {
            ans = (long long int)ans * a % n;
        }
        x = a;//从这里开始就是迭代乘法,验证二次验证定理 
        a = (long long int)a * a % n;//这里就相当于 x*x % m = 1 
        if (a == 1 && x != 1 && x != (n - 1))
        {
            return 1; // 这里需要注意,返回一的话就说明,追踪过程中,出现了不是素数的依据. 
        }
        t >>= 1;
    }
    return (ans == 1) ? 0 : 1;
}
        
 
int MillerRobin(int n, int s) // 一般s取50就可以避免所有的偶然性了. 
{
    if (n == 2)
    {
        return 1;
    }
    if (n < 2 || !(n & 1))
    {
        return 0;
    }
    int a;
    for (int i = 0; i < s; i++)
    {
        a = (long long int )rand() * (n - 2) / RAND_MAX + 1; //这样生成的随机数就是真正的随机数了 
        if (witness(a, n))
        {
            return 0;
        } 
    }
    return 1;
}
 
int main()
{
    while (scanf("%d", &N) != EOF)
    {
        if (N == 0)
        {
            break;
        }
        if (MillerRobin(N, 50))
        {
            printf("%d is a prime!
", N);
        }
        else
        {
            printf("%d is not a prime!
", N);
        }
    }
    return 0;
}

参考：https://blog.csdn.net/aledavvv/article/details/8929416