附: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)和质数表一样快,且省内存
继续努力~~