证明见《算法分析与设计——c++语言描述(第2版)》的第232页,写的很详细。
其中用到了快速幂来判定是否满足费马小定理,同时判定x*x%n==1的解是否只有x=1和x=n-1,这样Carmichael误判为素数的几率会更小
#include<bits/stdc++.h>
using namespace std;
int fast(int a,int n,int p){
//就是个快速幂,true则必为合数,false则以高概率为素数
int ret=1;
while(n>0){
if(n%