题意:1.p不是素数 2.(a^p)%p=a 输出yes 不满足输出no
思路:
bool is_prime(int n)
{ for(int i=2;i*i<=n;i++) if(n%i==0) return false;
return n!=1;}
2.(a^p)%p快速幂 注意大数快速幂要传入 long long 形参
解决问题的代码:
#include
#include
using namespace std;
typedef long long ll;
bool is_prime(int n)
{
for (int i = ; i*i <= n; i++)
if (n%i == ) //说明不是素数
return false;
return n != ; //忘记判断1
}
ll mod_p(ll a, ll b, ll c) //形参为long long 型不然wa
{
ll ans = 1ll;
a %= c;
while (b)
{
if (b & )
ans = (ans*a) % c;
a = (a*a) % c;
b >>= ;
}
return ans;
}
int main()
{
int p, a;
while (scanf("%d%d", &p, &a) != EOF)
{
if (p == && a == ) break;
if (!is_prime(p) && (mod_p(a, p, p) == a)) printf("yes\n");
else printf("no\n");
}
return ;
}
手机扫一扫
移动阅读更方便
你可能感兴趣的文章