HDU1905

附：HDU1905

题意：对a,p, 求解非质数p能否满足a^p%p==a。

思路：其实也就是质数判断和快速幂问题。

程序如下：

#include<stdio.h>
#include<math.h>
typedef __int64 ll;
/* 1：筛法
#define MAX 50000
#define MAXNUM 10000
int isprime[MAX],prime[MAXNUM];
void creatprimetable()
{
	for(ll k=4;k<MAX;k+=2)
		isprime[k]=1;
	for(ll i=3;i<=sqrt(MAX);++i)
	{
		if(isprime[i]==0)
		{
			for(ll j=i*i;j<MAX;j+=2*i)
				isprime[j]=1;
		}
	}
	for(ll m=2;m<MAX;++m)
		if(isprime[m]==0)
			prime[++prime[0]]=m;
}
int whetherprime(ll x)
{
	int flag=1;
	for(ll i=1;i<=prime[0]&&prime[i]<=sqrt(x);++i)
		if(x%prime[i]==0)
		{
			flag=0;
			break;
		}
	return flag;
}
*/

/* 2：一般方法：
ll whetherprime(ll x)
{
	if(x==2)
		return 1;
	if(x>2&&x%2==0)
		return 0;
    for(ll i=3;i<=sqrt(x);i+=2)
        if(x%i==0)
            return 0;
    return 1;
}
*/
// 最快0MS 米勒-罗宾算法依据费马小定理 对0<a<p,a^(n-1) mod n!=1,n必是合数。找几十次即可.
int whetherprime(ll x)
{
    ll randvalue;
    int n=TIMES;
    while(n--)
    {
        randvalue=rand()%x;
        if(randvalue!=0&&powermod(randvalue,x-1,x)!=1)
            return 0;
    }
    return 1;
}
ll powermod(ll a,ll b,ll c)
{
	a%=c;
	ll ans=1;
	while(b)
	{
		if(b&1)
			ans=(ans*a)%c;
		a=(a*a)%c;
		b=b>>1;
	}
	return ans;
}
int main()
{
	ll a,p;
//	creatprimetable();
	while(scanf("%I64d%I64d",&p,&a)!=EOF&&p|a)
	{
		if(whetherprime(p)==0&&powermod(a,p,p)==a)
			printf("yes\n");
		else
			printf("no\n");
	}
	return 0;
}
//注意:用64位__int64(两个_)才能过，int不一定能过.
//用直接2到sqrt(x)和质数表一样快，且省内存

继续努力~~