問題簡述：參見上述連結。

問題分析：

這個問題是驗證偽素數問題。p是偽素數的條件是，p不是素數並且滿足a^p=a(modp)。

偽素數是數論中與費爾馬小定理有關的一個重要概念。

程式說明：

函式isprime()不是一個真正意義上的素數判斷函式，只進行奇數判定，對於本題條件是沒有問題的。

函式powermod()是模冪計算函式。

AC的C++語言程式如下：

/* POJ3641 UVA11287 HDU1905 Pseudoprime numbers */

#include <iostream>
#include <cmath>

using namespace std;

typedef unsigned long long ULL;

// 試除法判斷一個數是否為素數
bool isprime(ULL n)
{
    ULL end2, i;

    end2 = sqrt(n);
    for(i=3; i<=end2; i+=2) {
        if(n % i == 0)
            break;
    }

    return i > end2;
}

// 模冪計算
ULL powermod(ULL a, ULL n, ULL m)
{
    ULL res = 1LL;
    while(n) {
        if(n & 1LL) {        // n % 2 == 1
            res *= a;
            res %= m;
        }
        a *= a;
        a %= m;
        n >>= 1;
    }
    return res;
}

int main()
{
    ULL p, a;

    while(cin >> p >> a && (p || a)) {
        if(!isprime(p) && powermod(a, p, p) == a % p)
            cout << "yes" << endl;
        else
            cout << "no" << endl;
    }

    return 0;
}