poj3641(快速幂取模)
来源:互联网 发布:浙江c语言二级考试时间 编辑:程序博客网 时间:2024/06/05 02:36
Description
Fermat's theorem states that for any prime number p and for any integer a > 1, ap = a (mod p). That is, if we raise a to the pth power and divide by p, the remainder is a. Some (but not very many) non-prime values of p, known as base-a pseudoprimes, have this property for some a. (And some, known as Carmichael Numbers, are base-a pseudoprimes for all a.)
Given 2 < p ≤ 1000000000 and 1 < a < p, determine whether or not p is a base-a pseudoprime.
Input
Input contains several test cases followed by a line containing "0 0". Each test case consists of a line containing p and a.
Output
For each test case, output "yes" if p is a base-a pseudoprime; otherwise output "no".
Sample Input
3 210 3341 2341 31105 21105 30 0
Sample Output
nonoyesnoyesyes
#include<iostream>#include<cstdio>#include<algorithm>#include<cmath>using namespace std;bool isPrime(int a){int num = (int)sqrt(a);for(int i = 2; i<=num; i++){if(a%i==0){return false;}}return true;}int quick_mod(int a, int p){int ret = 1;int mod = p;for( ; p; p>>=1, a = (long long)a*a%mod){if(p&1){ret = (long long)ret*a%mod;}}return ret;}int main(){int p, a;while(cin>>p>>a){if(p==0){break;}bool ok = false;if(!isPrime(p)){int num = quick_mod(a, p);if(num == a){ok = true;}}if(ok){cout<<"yes"<<endl;}elsecout<<"no"<<endl;}return 0;}
- poj3641(快速幂取模)
- POJ3641 Pseudoprime numbers(快速幂取模)
- 快速幂-poj3641
- POJ3641 快速幂
- POJ3641 Pseudoprime numbers 快速幂
- POJ3641 Pseudoprime numbers 【快速幂】
- Pseudoprime numbers(poj3641)快速幂
- poj3641
- poj3641
- poj3641
- poj3641
- poj3641
- POJ3641
- poj3641
- poj3641
- Pseudoprime numbers poj3641(快速幂+素数判定)
- poj3641(快速幂,以及素性测试方法的总结)
- poj3641 简单数论 快速幂 取模,素数
- wpa_supplicant--忘记网络后,wpa自动连接规则
- service能否实现不被Kill
- 错误记录--更改tomcat端口号方法,Several ports (8005, 8080, 8009)
- [LeetCode] Maximum Subarray
- Redis源码解析—源码目录介绍
- poj3641(快速幂取模)
- The architecture of REDIS
- redis设计思想
- thinkphp3.2开发排错一例——当输入数据长度超过表里设置的长度后,提示操作失败,没有提示超过长度
- java equals函数
- Cocos2dx 3.0正式版本UI解析错乱解决方案
- vim程序编辑器
- wikioi1066 引水入城 DFS
- codeforces A. Cakeminator 题解