poj3641

题意是输入p,m如果 p是素数输出no否则看m的p次方%p是否等于m，要用到快速幂取模，复习一下

#include <stdio.h>
#include <math.h>
__int64 modular(__int64 a,__int64 b,__int64 c)//快速幂取模(a的b次方取余c)
{
    __int64 ans=1;
    a=a%c;
    while(b>0)
    {
        if ((b&1)==1)//如果是奇数
 		ans=(ans*a)%c;//先将这个数乘进去
 		b>>=1;//b=b/2;
 		a=(a*a)%c;
    }
   return ans;
}
int isprime(__int64 x)
{
    __int64 i,j,k=(int)sqrt(x);
    if ((x==2)||(x==3))
    return 1;
    if (x==0)
    return 0;
    for (i=2;i<=k;i++)
    if (x%i==0)
    return 0;
    return 1;
}
 int main()
 {
     __int64 p,m,result;
     while(scanf("%I64d%I64d",&p,&m)==2)
     {
         if ((p==0)&&(m==0))
            break;
         if (isprime(p))
            printf("no\n");
         else
            {
            if (modular(m,p,p)==m)
              printf("yes\n");
            else
              printf("no\n");
            }
      }
      return 0;
}

//快速幂取模的模板
__int64 modular(__int64 a,__int64 b,__int64 c)//快速幂取模(a的b次方取余c)
{
    __int64 ans=1;
    a=a%c;
    while(b>0)
    {
        if ((b&1)==1)//如果是奇数
 		ans=(ans*a)%c;//先将这个数乘进去
 		b>>=1;//b=b/2;
 		a=(a*a)%c;
    }
   return ans;
}