UVA-10006 Carmichael Numbers

题目传送门

题目的意思就是给出n ，看它是否满足下面两个条件：

1. n是合数 。
2. 对于区间 [ 2 , n ) 的任意数a ，都使得 pow (a , n) % n == a成立。

首先判断n是不是合数，我使用素数筛进行一次性记录判断。
然后用快速幂遍历所有小于n的a（2<=a< n）,判断它是不是满足第二个条件.
注意使用long long 。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=65000+10;
bool isprime[maxn];
int prime[maxn];
int tot,mod;
ll mypow(ll a,ll n)
{
	ll res=(ll)1;
	while(n)
	{
		if(n&1)
		res=res*a%mod;
		a=a*a%mod;
		n>>=1;
	}
	return res;
}
void prime_match()
{
	tot=0;
	fill(isprime,isprime+maxn,true);
	isprime[0]=isprime[1]=false;
	for(int i=2;i<maxn;i++)
	{
		if(isprime[i])
		{
			prime[tot++]=i;
		}
		for(int j=0;j<tot&&i*prime[j]<maxn;j++)
		{
			isprime[i*prime[j]]=false;
			if(!(i%prime[j]))
			break;
		}
	}
}
int main()
{
	int n,i;
	prime_match();
	while(scanf("%d",&n)==1&&n)
	{
		if(isprime[n])
		{
			printf("%d is normal.\n",n);
			continue;
		}
		else
		{
			mod=n;
			for(i=2;i<n;i++)
			{
				int cnt=mypow(i,n);
				if(cnt!=i)
				break;
			}
			if(i>=n)
			{
				printf("The number %d is a Carmichael number.\n",n);
			}
			else
			{
				printf("%d is normal.\n",n);
			}
		}
	}
	return 0;
}