题目内容:找第n个默尼森数。
P是素数且M也是素数,并且满足等式M=2^P-1,则称M为默尼森数。
例如,P=5,M=2^P-1=31,5和31都是素数,因此31是默尼森数。
输入格式:
按提示用input()函数输入
输出格式:
int类型
输入样例:
4
输出样例:
127
题目分析:
输入的是第几个默尼森数,故先找P满足的素数,对应公式求出M。
再看M是否为素数,若M也是素数,则M是默尼森数。
P:2,3,5,7,11…
公式求出数:3,7,31,127,2047(不是素数),8191…
M:3,7,31,127,8191…
这里我犯了个错,以为素数P对应的数M都是默尼森数,结果错了。。。。。。2^7-1就不是
概要设计
一、判断素数函数
def prime(num):
if num<=1:
return False
for i in range(2,int(sqrt(num))+1):
if num%i==0:
return False
return True
二、默尼森函数
def monisen(no):
countM=0
for p in range(2,100):
if prime(p):
if prime(pow(2,p)-1):
countM=countM+1
if countM==no:
break
return pow(2,p)-1
参数no是你要找的第几个默尼森数
找100内的素数(可以改成更大范围),若P是素数,则判断M是否为素数。
若是,则countM增1。
若countM就是你要找的第几个数(no),则返回M
三、主函数
def main():
print(monisen(int(input())))
main()
完整代码
from math import sqrt
def prime(num):
if num<=1:
return False
for i in range(2,int(sqrt(num))+1):
if num%i==0:
return False
return True
def monisen(no):
countM=0
for p in range(2,100):
if prime(p):
if prime(pow(2,p)-1):
countM=countM+1
if countM==no:
break
return pow(2,p)-1
def main():
print(monisen(int(input())))
main()
第一次写CSDN。。。。。。做个记录