poj 3641 Pseudoprime numbers

#include <iostream>        //二分
using namespace std;
__int64 p;
bool prime(__int64 n)
{
for(__int64 i=2;i*i<=n;++i)
if(n%i==0)
return false;
return true;
}
__int64 query(__int64 a,__int64 n)
{
if(n==1)
return a;
__int64 q
=query(a,n/2),t=(q*q)%p;
if(n%2==0)
return t;
else
return (t*a)%p;
}
int main()
{
__int64 a;
while(cin>>p>>a&&p)
{
if(prime(p)||query(a,p)!=a)
cout
<<"no\n";
else
cout
<<"yes\n";
}
return 0;
}

  

转载于:https://www.cnblogs.com/mjc467621163/archive/2011/07/22/2114364.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值