poj 3641

#include <iostream>
#include <cstdio>
#include <iomanip>
#include <string>
#include <cstdlib>
#include <cstring>
#include <queue>
#include <set>
#include <vector>
#include <map>
#include <algorithm>
#include <cmath>
#include <stack>

#define INF 0x3f3f3f3f
#define LINF 0x3f3f3f3f3f3f3f3f
#define ll long long
#define ull unsigned long long
#define uint unsigned int
#define l(x) (x<<1)
#define r(x) (x<<1|1)
#define ms(a,b) memset(a,b,sizeof(a))

using namespace std;



bool isp(ull x) {
	for (ull i = 2; i*i <= x; i++)
		if (x%i == 0)return false;
	return true;
}

ull mpow(ull a, ull p) {
	ull mod = p;
	ull res = 1;
	while (p) {
		if (p & 1) res = (res*a) % mod;
		a = (a*a)%mod;
		p >>= 1;
	}
	return res;
}

ull p, a;

int main() {
	while (cin >> p >> a, p != 0) {
		if (isp(p)) cout << "no" << endl;
		else 
			if (mpow(a, p)==a%p) cout << "yes" << endl;
			else cout << "no" << endl;
	}

	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值