判断伪素数~ 先判断是不是素数,然后用快速幂 判断。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <cmath>
#include <cstdlib>
#include <queue>
#include <stack>
using namespace std;
long long pp;
long long pow1(long long p,long long a)
{
    if(p==0)
        return 1;
    long long x=pow1(p/2,a);
    long long ans =x*x%pp;
    if(p%2==1)
        ans=ans*a%pp;


    return ans;
}


int main()
{
    bool has=false;
    long long p,a,b;
    while(cin>>p>>a)
    {


        has=true ;
       if(p==0&&a==0)
        break;




      for(int i=2;i*i<p;i++)
            {
                if(p%i==0)
                {
                     has=false;
                     break;
                }
            }


        if(has==false )
         {
             pp=p;
            b=pow1(p,a);
            if(a==b)
            {
                printf("yes\n");
            }
            else
            {
                printf("no\n");


            }
         }
         else
         {
             printf("no\n");
         }


    }
return 0;
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值