梅森素数:nefu 120 梅森素数(卢卡斯—莱默判别法+大数乘积取模)

本文介绍了梅森素数的一种判定方法——Lucas-Lehmer判别法,并结合大数乘积取模算法,简化计算过程。通过解题思路和参考代码,阐述了如何利用这种方法来验证梅森数是否为素数,并附带了前47个梅森素数的列表。
摘要由CSDN通过智能技术生成

梅森素数


description

由于梅森学识渊博,才华横溢,为人热情以及最早系统而深入地研究2p-1 型的数(其中p为素数),为了纪念他,数学界就把这种数称为“梅森数”;并以Mp 记之(其中M为梅森姓名的首字母),即Mp=2p-1 。如果梅森数为素数,则称之为“梅森素数”。 比如p=2,3,5,7时,Mp都是素数,但211-1 不是素数 。现在请你求出前N个梅森素数。

input

有多组测试数据。
第一行是一个正整数T,表示测试数据的组数。接下来每组1个数p的值,这里2<= p <= 62。

output

对于每组测试数据,判断Mp 是不是梅森素数,是就输出“yes ”,否就输出“no”,输出后要换行。

sample_input

2
2
7

sample_output

yes
yes

解题思路:

1.梅森素数的判定方法之一——Lucas—Lehmer判别法

构造rk序列判断Mp是否是素数,同时应用大数乘积取模的算法,将乘法变成加法

参考代码+部分解释:

#include 
   
   
    
    
#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
#include 
       
       
         #include 
        
          #include 
         
           using namespace std; const int maxn=50000+10; long long r[100],p; long long multi(long long a,long long b,long long mod)//a * b % mod 乘法换加法 { long long ret=0; while(b){ if(b&1) ret=(ret+a)%mod; b>>=1; a=(a<<1)%mod; } return ret; } int main() { // freopen("input.txt","r",stdin); int T;cin>>T; while(T--){ cin>>p; long long mod=((long long)1< 
           
          
         
       
     
     
    
    
   
   
2.枚举判断——梅森素数表
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
        
        
#include 
        
        
#include 
         
         
           #include 
          
            #include 
           
             using namespace std; const int maxn=50000+10; int p,Mp[9]={2,3,5,7,13,17,19,31,61}; bool solve() { for(int i=0;i<9;i++) if(p==Mp[i]) return true; return false; } int main() { // freopen("input.txt","r",stdin); int T;cin>>T; while(T--){ cin>>p; if(solve()) cout<<"yes"< 
             
            
           
         
       
       
      
      
     
     

附表:梅森素数表(1-47个)

序号

p

Mp=(2^p)-1

Mp的位数

发现时间

发现者

1

2

3

1

古代

古人

2

3

7

1

古代

古人

3

5

31

2

古代

  • 1
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值