Rabin-Miller素性测试算法

算法原理：

Th1 p为素数, 0

证明：

a mod p =1或者 a mod p =-1

Th2 任何一个素数p都能用表示

根据费马小定理 对于素数p和1

所以

由此可以得到数列中

（1）所有值都为1 或者

（2）有的值不为1 但是其平方后模p为1 所以该值为 -1mod p =p-1

如果一个数不满足上诉条件，则这个数一定为合数；

如果一个数满足上诉条件，则它可能是素数，也可能是合数，但它是素数的概率（0.75）大于

它是合数的概率。为了降低误判（将合数判为合数），多次运行该算法，降低错误率

算法实现：

bool millerRabinPrimeTest(unsigned long p,int cnt){//对于素数p进行检验 ，检验cnt次
	unsigned long  k=0;
	unsigned long  q= p-1;
	unsigned long r=0;

	if(q!=1&&q%2==1)return false;
	for (;q%2!=1;k++)q/=2;
	for (int i=0;i>=1;
		t=(t*t)%n;
	}

如果要求算法的错误率小于 t 则