首先，显然的，$a^{xy} =(a^x)^y$ 。
然后可以假设其中 $a^x ≡1(modp)$，则有 $a^{xy} =(a^x)^ y ≡1^y ≡1(modp)$。
所以我们可以 $O(\sqrt{p}​)$分解 $p−1$ 的质因数，然后检查 $p−1$ 的每一个因数 $f$ 是否满足 $a^ f ≡1(modp)$。
如果有任意一个因数满足上述条件，直接输出 No。
如果所有因数均不满足，输出 Yes。总复杂度 $O(T\sqrt{p}logp)$。