POJ3641 Pseudoprime numbers 快速幂

最新推荐文章于 2024-10-06 14:10:45 发布

LoveKobe_

最新推荐文章于 2024-10-06 14:10:45 发布

阅读量462

点赞数

分类专栏：快速幂运算文章标签： algorithm poj 算法

本文链接：https://blog.csdn.net/u011363774/article/details/37833299

版权

快速幂运算专栏收录该内容

2 篇文章 0 订阅

订阅专栏

这题比较简单，用快速幂进行加速幂运算。

#include <stdio.h>
#include <vector>
#include <math.h>
#include <string.h>
#include <string>
#include <iostream>
#include <queue>
#include <list>
#include <algorithm>
#include <stack>
#include <map>

#include<iostream>  
#include<cstdio>  
using namespace std;

long long qpowmod(long long x, int n, int mod)
{
	long long res = 1;
	while (n > 0)
	{
		if (n & 1)
		{
			res *= x;
			res %= mod;
		}

		x *= x;
		x %= mod;
		n >>= 1;
	}
	return res;
}

bool isprime(int a)
{
	for (int i = 2; i *i <= a;i++)
	{
		if (a % i == 0)
		{
			return false;
		}
	}
	return true;
}

int main()
{
	int p, a;
#ifdef _DEBUG
	freopen("d:\\in.txt", "r", stdin);
#endif
	while (scanf("%d %d\n", &p, &a) != EOF)
	{
		if (p == 0 && a == 0)
		{
			break;
		}
		if (isprime(p))
		{
			printf("no\n");
			continue;
		}
		if (qpowmod(a, p, p) == a)
		{
			printf("yes\n");
		}
		else
			printf("no\n");
	}
	return 1;
}