蒙特卡罗法判断素数：

蒙特卡罗法判断素数的原理是：对于给定的正整数n，首先随机生成0到n-1之间的两个整数a和b，然后计算a^2-b^2的值。如果a^2-b^2能被n整除，则n不是素数；如果a^2-b^2不能被n整除，则n可能是素数。随着蒙特卡罗测试的次数增加，判断结果的可靠性也会提高。存在一定的误差概率，使用于一些较小的数。并不实用，即做分享。

#include <iostream>
#include <cstdlib>
#include <ctime>
using namespace std;

bool isPrime(int n, int k)
{
    // 排除n为1和2的情况
    if (n <= 2)
    {
        return (n == 2);
    }

    srand(time(NULL)); // 初始化随机数种子

    for (int i = 0; i < k; i++)
    { // 进行k次蒙特卡罗测试
        int a = rand() % (n - 3) + 2; // 生成2到n-1之间的随机数a
        int x = 1;
        int power = a;

        while (power > 0)
        { // 使用快速幂算法计算 a^n mod n 的值
            if (power & 1)
            {
                x = (long long)x * a % n;
            }
            a = (long long)a * a % n;
            power >>= 1;
        }

        if (x != 1)
        { // 如果 a^(n-1) mod n 不等于 1，则 n 可能不是素数
            return false;
        }
    }

    return true; // 经过 k 次蒙特卡罗测试，n 可能是素数
}

int main()
{
    int n, k;
    cout << "请输入一个正整数：";
    cin >> n;
    cout << "请输入进行蒙特卡罗测试的次数：";
    cin >> k;

    if (isPrime(n, k))
    {
        cout << n << "是质数" << endl;
    }
    else
    {
        cout << n << "不是质数" << endl;
    }

    return 0;
}