UVA 10006 Carmichael Numbers（快速幂）

最新推荐文章于 2022-02-14 17:48:33 发布

yphacker

最新推荐文章于 2022-02-14 17:48:33 发布

阅读量508

点赞数

CC 4.0 BY-SA版权

分类专栏： --------------数论-------------- 文章标签：素数筛选法快速幂

本文链接：https://blog.youkuaiyun.com/piaocoder/article/details/47451109

--------------数论-------------- 专栏收录该内容

39 篇文章

订阅专栏

本文介绍了一种用于检测Carmichael数的算法，并提供了完整的AC代码实现。该算法首先通过素数筛选法判断给定的数是否为素数，然后使用快速幂取模的方法验证该数是否满足Carmichael数的定义。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Carmichael Numbers

题目链接：

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=12&page=show_problem&problem=947

解题思路：

输入n，若满足如下两个条件，则n是Carmichael number

（１）n不是素数

（２）对于所有a(2<=a<n)，有(a^n)%n = a

先判断是否为素数，用素数筛选法，然后用快速幂取模，注意运算过程中的乘法溢出int ，要使用long long...

AC代码：

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

typedef long long ll;
int not_prime[65005];

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

int main(){
    memset(not_prime,0,sizeof(not_prime));
    for(int i = 2; i <= 65000; i++){
        int tmp = 65000/i;
        for(int j = 2; j <= tmp; j++)
            not_prime[i*j] = 1;
    }
    int n;
    while(scanf("%d",&n),n){
        if(!not_prime[n]){
            printf("%d is normal.\n",n);
            continue;
        }
        int flag = 1;
        for(int i = 2; i < n; i++){
            if(mod_pow(i,n,n) != i){
                flag = 0;
                printf("%d is normal.\n",n);
                break;
            }
        }
        if(flag)
            printf("The number %d is a Carmichael number.\n",n);
    }
    return 0;
}