关于解决快速判断longlong型的整数是否为素数

VIP文章 llqxs

于 2017-04-13 21:27:43 发布

阅读量1.5k

点赞数 2

分类专栏：数论文章标签：米勒罗宾定理

本文链接：https://blog.csdn.net/qq2442438699/article/details/70162294

版权

核心做法是米勒罗宾定理。

测试n（n很大）

在（1，n-1）中选出一个a，1<=a<=n-1；

如果满足a^(n-1)%n==1;那么n有1/2的概率是素数；

如果不满足，那么n一定不是素数；

只需选取多个a，那么测试的正确率会无限接近100%

// 18位素数：154590409516822759  
#include<stdio.h>
typedef  unsigned long long ll;
ll mul(ll a,ll b,ll n)
{
    ll ans=0;
    while(b)
    {
        if(b&1) ans=(ans+a)%n;
        a=(a+a)%n;
        b=b>>1;
    }
    return ans;
}
ll Pow(ll a,ll b,ll n)
{
    ll result=1;
    ll base=a%n;
    while(b)
    {
        if(b&1) result=mul(result,base,n)%n;
        base=mul(base,base,n)%n;
        b=b>>1;
    }
    return result;
}
main()
{
    ll n;
    scanf("%llu",&n);
    if(n<=73)
    {
        int i;
        for(i=2; i*i<=n; i++)
            if(n%i==0)