POJ 3641 Pseudoprime numbers G++ 快速幂运算巧妙 检测素数 背

 

#include <iostream>
#include <cstdio>
using namespace std;
//谢谢博友好文章  抄博友程序  英语  快速求幂巧妙   检测素数   背 
__int64 mi(__int64 a,__int64 b)//抄博友程序   背 
{
	__int64 ba=b;
	__int64 jg=1;
	while(b>0)
	{
		if(b&1)//抄博友程序 
		{
			jg=(jg*a)%ba;
		}
		a=(a*a)%ba;
		b=b>>1;
	}
	return jg;
} 
int fun(__int64 a)//抄博友程序 
{
	for(int i=2;i*i<a;i++)
	{
		if(a%i==0)
		{
			return 0;//不是素数 
		}
	}
	return 1;
}
int main()
{
	while(1)
	{
		__int64 p, a;
		cin>>p>>a;
		if(p==0&&a==0)
		{
			break;
		}
		__int64 da=mi(a,p); 
		int flag=fun(p);
		if(da==a&&flag==0)
		{
			cout<<"yes"<<endl;
		}else
		{
			cout<<"no"<<endl;
		}
	} 
	return 0;
}