POJ-3641 Pseudoprime numbers(快速幂）

  题目：POJ-3641 Pseudoprime numbers

  题目大意：给了base-a Pseudoprime numbers一个定义，就是首先某个数p不是素数，然后这个数对于一个数a,满足a^p mod p=a。现在给定两个数p,a让我们判断p是否为base-a Pseudoprime numbers。

  解题思路：按照定义判断即可，p比较大，所以用到了快速幂。

  AC代码：

#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
bool isprime(long long x){
    for (long long i=2;i<=sqrt(x);i++){
        if (x%i==0)return 0;
    }
    return 1;
}
long long qpow(long long a,long long b,long long m){
    long long ans=1;
    long long k=a;
    while (b){
        if (b&1)ans=ans*k%m;
        k=k*k%m;
        b>>=1;
    }
    return ans;
}
int main(){
    long long i,j,a,p;
    while (scanf("%lld%lld",&p,&a)!=EOF&&(p||a)){
        if (isprime(p)){
            printf("no\n");
        }
        else {
            if (qpow(a,p,p)%p==a){
                printf("yes\n");
            }
            else printf("no\n");
        }
    }
}