Pseudoprime numbers伪素数（快速幂+判定素数）

http://acm.sdut.edu.cn/sdutoj/showproblem.php?pid=1469&cid=1142

#include<stdio.h>  
#include<math.h>  
int isprime(long n)  
{  
    long i,k;  
    k=sqrt(n);  
    for(i=2;i<=k;i++)  
    {  
        if(n%i==0)  
         return 0;  
    }  
    return 1;  
}  
long mod(long a,long n,long p)  
{  
    long t;  
    if(n==0) return 1%p;  
    if(n==1) return a%p;  
    t=mod(a,n/2,p);  
    t=t*t%p;  
    if(n%2==1) t=t*a%p;  
    return t;  
}  
int main()  
{  
   long p, a ;  
   while(scanf("%ld %ld",&p, &a),p!=0&&a!=0)  
   {  
      if(isprime(p))  
      printf("no\n") ;  
      else  
      {  
          if(a==mod(a, p, p))  
           printf("yes\n") ;  
          else printf("no\n") ;  
      }  
   }  
  return 0;  
}