UVA 10006(素数判断+快速幂)

My_day_Hyt

于 2016-08-27 22:25:06 发布

阅读量521

点赞数

CC 4.0 BY-SA版权

分类专栏： ACM-数论文章标签： uva 算法

本文链接：https://blog.youkuaiyun.com/My_day_Hyt/article/details/52337423

ACM-数论专栏收录该内容

1 篇文章

订阅专栏

本文介绍了一种用于检测Carmichael数的算法实现。通过筛选法预先标记合数，并利用快速幂模运算验证给定数值是否为Carmichael数。文章提供了完整的C++源代码，展示了如何判断一个数是否正常或属于特殊的Carmichael数。

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

/*过两天写这两个算法的详解*/
#include <iostream>
#include<stdio.h>
#include<algorithm>
#include<math.h>
#include<string.h>
using namespace std;
#define ll long long
ll ans = 1;
int n;
const int maxn = 65005;
int a[maxn];
void judge()
{
    int m = sqrt(65000+0.5);
    memset(a,0,sizeof(a));
    for(int i=2;i<=m;i++)
    {
        if(!a[i])
        {
            for(int j=i*i;j<65010;j+=i)
                a[j]=1;
        }
    }
}
int powmod(int i,int n,int m)
{
    int res = 1;
    int tmp = i;
    while(n>0)
    {
        if(n&0x1)   res *= tmp;
        tmp *= tmp;
        n>>=1;
    }
    res = res%m;
    return res;
}
int main()
{
    judge();
    while(scanf("%d",&n)&&n!=0)
    {
        if(a[n] != 1)   printf("%d is normal.\n",n);
        else{
            bool flag = true;
            for(int i=2;i<n;i++)
            {
                ans = powmod(i,n,n);
                if(ans != i)    {flag = false;break;}
            }
            if(flag == 0)    printf("The number %d is a Carmichael number.\n",n);
                else    printf("%d is normal.\n",n);
        }
    }
    return 0;
}