Carmichael Numbers (Uva No.10006) -- 快速幂运算_埃氏筛法_打表

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>
#include <cstring>
using namespace std;

typedef long long LL;
const int maxn =  + ;

//int prime[maxn + 1];  //第i个素数,保存区间内素数
bool is_prime[maxn];  //is_prime[i]为true表示i是素数
bool judge[maxn];     //能随机访问某数是否为素数 

void sieve(int n) {
    memset(judge, , sizeof(judge));
//    int p = 0;
    for (int i = ; i <= n; i++) is_prime[i] = true;
    is_prime[] = is_prime[] = false;
    for (int i = ; i <= n; i++) {
        if (is_prime[i]) {
//            prime[p++] = i;
            judge[i] = true;
            for (int j = *i; j <= n; j+=i) is_prime[j] = false;
        }
    }
}

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

bool solve(int n)
{
    for (int i = ; i < n; i++) {
        if (mod_pow(i, n, n*) != i)
            return false;
    }
    return true;
}

int main()
{
    int n;
    sieve(maxn);
    while (scanf("%d", &n) != EOF && n)
    {
        if (!judge[n] && solve(n)) printf("The number %d is a Carmichael number.\n", n);
        else printf("%d is normal.\n", n);
    }
    return ;
}