poj3641

最新推荐文章于 2020-12-16 07:44:56 发布

幽涧恋桥

最新推荐文章于 2020-12-16 07:44:56 发布

阅读量468

点赞数

分类专栏： ACM 数论文章标签： poj

本文链接：https://blog.youkuaiyun.com/u013614723/article/details/44964147

版权

ACM 同时被 2 个专栏收录

37 篇文章

订阅专栏

数论

8 篇文章

订阅专栏

题目链接：http://poj.org/problem?id=3641

题意：a^p = a (mod p)

解题思路：快速幂运算

#include <iostream>
#include <math.h>
using namespace std;

long long mod_pow(long long x,long long n,long long mod)
{
    if(n==0)
        return 1;
    long long res=mod_pow(x*x%mod,n/2,mod);
    if(n&1)
        res=res*x%mod;
    return res;
}

bool isPrime(long long p)
{
    if(p==2)
        return true;
    if(p%2==0||p==1)
        return false;
    bool flag=true;
    for(int i=3;i<=sqrt(p);i=i+2)
    {
        if(p%i==0)
        {
            flag=false;
            return flag;
        }

    }
    return flag;
}

int main()
{
    long long p,a,ans;
    while(cin>>p>>a,p!=0||a!=0)
    {
        ans=mod_pow(a,p,p);
        if(ans==a&&!isPrime(p))
            cout<<"yes"<<endl;
        else
            cout<<"no"<<endl;
    }
    return 0;
}