数论（米勒测试）

a^(p-1)%p==1

 

 

#include  < iostream >
#include  < time.h >
using   namespace  std;


inline __int64 GetMod(__int64  base ,__int64 power,__int64 Mod)
{
    __int64 t = 1 ;
    __int64 res = base ;
    
     while  (power > 1 )
    {
         if  (power & 1 )
        {
            t = (t * res) % Mod;
        }
        
        res = (res * res) % Mod;
        power = power >> 1 ;
        
    }
    
     return  (res * t) % Mod;
}

inline  bool  slove(__int64 n)
{
     if  (n == 1 )
    {
         return   false ;
    }
    
     if  (n == 2 )
    {
         return   true ;
    }
    
     if  ((n & 1 ) == 0 )
    {
         return   false ;
    }
    
    
     long  a;
    
    a = 2 ;
    if(n%a==0) return false;
     if  (n > a && GetMod(a,n - 1 ,n) != 1 )
    {
         return   false ;
    }
    
    a = 3 ;
    if(n%a==0) return false;
     if  (n > a && GetMod(a,n - 1 ,n) != 1 )
    {
         return   false ;
    }
    
    
    a = 5 ;
    if(n%a==0) return false;
     if  (n > a && GetMod(a,n - 1 ,n) != 1 )
    {
         return   false ;
    }
    
    a = 7 ;
    if(n%a==0) return false;
     if  (n > a && GetMod(a,n - 1 ,n) != 1 )
    {
         return   false ;
    }
    
     return   true ;
}

int  main()
{
    __int64 n;
     long  T;
    scanf( " %ld " , & T);
     while  (T -- )
    {
        scanf( " %I64d " , & n);
         if (slove(n))
        {
            puts( " yes " );
        }
         else
        {
            puts( " no " );
        }    
    }
    
     return   0 ;
}