素数的判定

3.素数判定
(1)试除法:
bool  isprime(int n)
{
 int i;
 for(i=2;i<=(int)sqrt(n*1.0);i++)
     if(n%i==0)    return  false;
 return  true;
}
bool isprime(int n)
{
 if(n==2) return true;
 if(n==1||(n&1)==0)  return false;
 for(int i=3;i*i<=n;i+=2)   if(n%i==0) return fals;
 return true;
}
(2)miller-rabin 算法
bool witness(i64 a,i64 n)
{
  i64 x,d=1,i=ceil(log(n-1.0)/log(2.0))-1;
  for(;i>=0;i--)
  {
    x=d; d=(d*d)%n;
    if(d==1&&x!=1&&x!=n-1) return 1;
    if(((n-1)&(1<<i))>0) d=(d*a)%n;
  }
  return  d==1?0:1;
}
bool miller_rabin(i64 n)
{
  if(n==2)   return 1;
  if(n==1||(n&1)==0) r