poj-3641 Pseudoprime numbers(费马小定理)

最新推荐文章于 2018-12-12 13:41:41 发布

双之城

最新推荐文章于 2018-12-12 13:41:41 发布

阅读量548

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/u011026037/article/details/23697039

版权

数论专栏收录该内容

46 篇文章 0 订阅

订阅专栏

判断a ^ p = a % p (p是合数) true?yes：no

#include<stdio.h>
#include<math.h>
typedef long long  LL;
bool is_prime(LL p)
{
    if(p == 1|| p == 2) return true;
    for(LL i = 2;i <= sqrt(p);i++)
        if(p%i==0) return false;
    return true;
}
LL fast_pow(LL x,LL p,LL MOD)//快速幂
{
    LL res = 1;
    while(p)
    {
        if(p&1) res = (res * x) % MOD;
        x = (x*x)%MOD;
        p>>=1;
    }
    return res;
}
int main()
{
    LL a,p;
    while(scanf("%lld%lld",&p,&a),p|a)
    {
        if(is_prime(p)) {printf("no\n");continue;}
        LL r = fast_pow(a,p,p);
        if(r == a%p) printf("yes\n");//此处不能用r == 1作为判断
        else printf("no\n");
    }
}