读书人

uva 10006 Carmichael Numbers(高速幂

发布时间: 2013-10-22 16:16:51 作者: rapoo

uva 10006 Carmichael Numbers(快速幂)

题目链接:uva 10006 Carmichael Numbers


题目大意:判断一个数是否为Carmichael数, (非素数, 并且满足a^n % n == a, a 的取值为2 ~ n - 1)。


解题思路:Eratosthenes筛选法求出素数,然后对应n如果为非素数,就对每个a进行判断,中间用到快速幂。


#include <stdio.h>#include <string.h>const int N = 65005;int n, prime[N];void init() {memset(prime, 0, sizeof(prime));prime[0] = prime[1] = 1;for (int i = 2; i < N; i++)for (int j = i * 2; j < N; j += i)prime[j] = 1;}long long count(int a, int c) {if (c == 1) return a;long long p = count(a, c / 2);p = p * p;if (c % 2) p *= a;return p % n;}bool judge(int n) {if (!prime[n]) return false;for (int i = 2; i < n; i++)if (count(i, n) != i) return false;return true;}int main () {init();while (scanf("%d", &n) && n) {if (judge(n))printf("The number %d is a Carmichael number.\n", n);elseprintf("%d is normal.\n", n);}return 0;}


读书人网 >编程

热点推荐