nefu 120 Lucas-Lehmer 梅森素数判别法二分-大数乘法换加法

最新推荐文章于 2022-08-24 11:28:25 发布

weixin_34096182

最新推荐文章于 2022-08-24 11:28:25 发布

阅读量85

点赞数

CC 4.0 BY-SA版权

文章标签： php

原文链接：http://www.cnblogs.com/814jingqi/p/3310440.html

本文介绍了一种利用卢卡斯莱默判别法判断梅森素数的有效方法，并通过具体的C++代码实现展示了该算法的工作原理。文章详细解释了如何通过位运算提高算法效率，避免数据溢出等问题。

题目地址：http://acm.nefu.edu.cn/test/problemshow.php?problem_id=120

#include<iostream>
#include<cmath>
#include<cstdlib>
using  namespace std;
typedef long long inta;



inta multi(inta a,inta b, inta m)
{
  inta rel=0;
  while(b>0)
  {
     if(b&1)
       {
           rel=(rel+a)%m;
       }
     b>>=1;
     a=(a<<1)%m;
  }
  return rel;
}
int test(int n)
{
    int count=0;
   while(n>0)
    {n>>=1;
      count++;
    }
 return count;
}
int main()
{

  int size=0;

   cin>>size;
   int p;
   while(cin>>p)
   {

      if(p<=3)
      cout<<"yes"<<endl;

      else
      {
          inta mp=1;
          mp<<=p;
          mp-=1;


        inta r=4;
        for(int i=0;i<p-2;i++)

         {

           r=(multi(r,r,mp)-2)%mp;


         }
        if(r==0) cout<<"yes"<<endl;
        else cout<<"no"<<endl;
      }
   }
}

注意要点

1 梅森素数判断用卢卡斯莱默判别法

2 long long 进行乘法时直接乘会溢出需要转化成加法边加边取模

3 使用位运算加法这样算法复杂度为O（log2（n））才不会超时

转载于:https://www.cnblogs.com/814jingqi/p/3310440.html