UVA 10006 - Carmichael Numbers (模平方)

题意:对任意的2<=x<n都有x^n=x mod n成立的合数n称为Carmichael number，给定一个整数n，判断它是不是Carmichael number。

注意：前提条件是合数，用模平方法枚举所有的x，乘法会超int。

#include <iostream>
#include<cstdio>
using namespace std;
#define N 65005
int prime[N]={0};
int mod_pow(long long  x,long long  n,long long  mod);
int main()
{
    prime[1]=1;
    for(int i=2;i<N;i++)
    {
        if(prime[i]==0){
        for(int j=i+i;j<N;j=j+i)
        {
            prime[j]=1;
        }
        }
    }
    int n;
    while(~scanf("%d",&n)&&n)
    {
        if(prime[n])
        {
            int i;
            for(i=2;i<n;i++)
            {
                if(mod_pow(i,n,n)!=i)
                break;
            }
            if(i==n)
                printf("The number %d is a Carmichael number.\n",n);
            else
                printf("%d is normal.\n",n);
        }
        else
            printf("%d is normal.\n",n);
    }

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