uva 10006

一题如此基础的题目，搞了那么久，晕。。。不过有很多细节错误自己改过来了，给自己鼓鼓掌^_^

求幂的时候要用longlong型的，可是为什么啊？x的值根本就不可能比65000大啊，每次都取模啊，求反例。。。。。。

还有判断素数要在前，刘汝佳先生说过&&是短路来着。。。。。。

#include <cstdio>
#include <cstring>
#include <cmath>

int prime[65002];
int n;

void is_prime()
{
    int m = sqrt(65001 + 0.5);
    memset(prime,0,sizeof(prime));
    prime[0] = prime[1] = -1;
    for(int i = 2;i <= m;i ++)//这里判断素数也挺快滴，也是刘汝佳先生的白皮书上滴。
    {
        if(!prime[i])
        {
            prime[i] = 1;
            for(int j = i*i;j <= 65001;j += i)
            {
                prime[j] = -1;
            }
        }
    }
}

long long pow_mod(int a,int k,int m)//这个求幂挺快的logk,刘汝佳先生的那本白皮书上滴
{
    if(k == 1)
        return (long long)a;

    long long x = pow_mod(a,k/2,m);
    long long ans = x * x % m;
    if(k % 2 == 1)
    {
        ans = ans * a % m;
    }

    return ans;
}

bool fermat(int k)
{
    for(int i = 2;i < k;i ++)
    {
        if(pow_mod(i,k,k) != i)
            return false;
    }

    return true;
}
int main()
{
    is_prime();
    while(scanf("%d",&n) == 1)
    {
        if(!n)
            break;
        if((prime[n] == -1)&&fermat(n))
        {
            printf("The number %d is a Carmichael number.\n",n);
        }
        else 
        {
            printf("%d is normal.\n",n);
        }
    }

    return 0;
}