uvaoj 10006 Carmichael Numbers 快速幂取模

uvaoj 10006 Carmichael Numbers 快速幂取模

这个题目是让判断给定的一个数是不是 Carmichael Numbers，即卡米切尔数，这种数满足费马小定理，但是不是素数。

满足下面条件的数n是卡米切尔数，对于任意的a(2<=a<n)有：a^n % n = a。因为要求幂取模，所以，可以使用递归来二分，或者使用快速幂来写。

代码如下：

/*************************************************************************
	> File Name: 10006.cpp
	> Author: gwq
	> Mail: gwq5210@qq.com 
	> Created Time: 2015年01月05日 星期一 23时05分56秒
 ************************************************************************/

#include <cmath>
#include <ctime>
#include <cctype>
#include <climits>
#include <cstdio>
#include <cstdlib>
#include <cstring>

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <string>
#include <vector>
#include <sstream>
#include <iostream>
#include <algorithm>

#define INF (INT_MAX / 10)
#define clr(arr, val) memset(arr, val, sizeof(arr))
#define pb push_back
#define sz(a) ((int)(a).size())

using namespace std;
typedef set<int> si;
typedef vector<int> vi;
typedef map<int, int> mii;
typedef long long ll;

const double esp = 1e-5;

#define N 65010

int isnpri[N];

ll pow_mod(ll a, ll n, ll b)
{
	if (n == 1) {
		return a % b;
	} else {
		ll tmp = pow_mod(a, n / 2, b) % b;
		ll ret = tmp * tmp % b;
		if (n % 2 == 0) {
			return ret;
		} else {
			return (ret * a) % b;
		}
	}
}

int judge(ll n)
{
	if (!isnpri[n]) {
		return 0;
	}

	for (int i = 2; i < n; ++i) {
		if (pow_mod(i, n, n) != i) {
			return 0;
		}
	}

	return 1;
}

int main(int argc, char *argv[])
{
	int n;

	clr(isnpri, 0);
	isnpri[0] = isnpri[1] = 1;
	for (int i = 2; i * i < N; ++i) {
		if (!isnpri[i]) {
			for (int j = i * i; j < N; j += i) {
				isnpri[j] = 1;
			}
		}
	}
	while (scanf("%d", &n) != EOF) {
		if (n == 0) {
			break;
		}
		if (judge(n)) {
			printf("The number %d is a Carmichael number.\n", n);
		} else {
			printf("%d is normal.\n", n);
		}
	}

	return 0;
}