Miller Robin大素数判定
Miller Robin算法
当要判断的数过大,以至于根n的算法不可行时,可以采用这种方法来判定素数。
用于判断大于2的奇数(2和偶数需要手动判断),是概率意义上的判定,因此需要做多次来减少出错概率。
Template:
```C++
typedef long long ll;
ll kmul(ll a,ll b,ll mod)
{
	ll res=0;
	while (b)
	{
		if (b&1)
			res=(res+a)%mod;
		a=(a+a)%mod;
		b>>=1;
	}
	return res;
}
ll kpow(ll a,ll b,ll mod)
{
	ll res=1;
	while (b)
	{
		if (b&1)
			res=kmul(res,a,mod)%mod;
		a=kmul(a,a,mod)%mod;
		b>>=1;
	}
	return res;
}
bool Mil_Rb(ll n,ll a)
{
	ll d=n-1,s=0,i;
	while (!(d&1))
	{
		d>>=1;
		s++;
	}
	ll t=kpow(a,d,n);
	if (t==1||t==-1)
		return 1;
	for (i=0;ia[i]&&!Mil_Rb(n,a[i]))
			return 0;
	}
	return 1;
}
```
 
                     
                    
                 
                    
                
 
 
                
            
         
         浙公网安备 33010602011771号
浙公网安备 33010602011771号