这篇文章是双重的。我是 Python 的新手。
第一部分。
这是我的旧谷歌谜题代码:“e 的连续数字中的前 10 位素数”(https://google-tale.blogspot.ca/2008/07/google-billboard-puzzle.html)
euler = '7182818284590452353602874713526624977572470936999595749669676277240766303535475945713821785251664274274663919320030599218174135966290435729003342952605956307381323286279434907632338298807531952510190115738341879307021540891499348841675092447614606680822648001684774118537423454424371075390777449920695517027618386062613313845830007520449338265602976067371132007093287091274437470472306969772093101416928368190255151086574637721112523897844250569536967707854499699679468644549059879316368892300987931277361782154249992295763514822082698951936680331825288693984964651058209392398294887933203625094431173012381970684161403970198376793206832823764648042953118023287825098194558153017567173613320698112509961818815930416903515988885193458072738667385894228792284998920868058257492796104841984443634632449684875602336248270419786232090021609902353043699418491463140934317381436405462531520961836908887070167683964243781405927145635490613031072085103837505101157477041718986106873969655212671546889570350354'
for i in range(len(euler) - 3):
if (euler[i + 3]== '1' or euler[i + 3]== '3' or euler[i + 3]== '7' or euler[i + 3]== '9'):
number = int(euler[i:i + 3])
for a in range(2, round(number**0.5)):
if number % a == 0:
break
else:
print(number),
break
奇怪的是,使用上面的代码我可以找到正确的答案,即 7427466391(编译需要 4500 万)但是当我测试 3 位数的代码时,它给了我 709,这是不正确的,因为前 3 位数的素数是 281 范围( 3,6)。
我怎样才能解决这个问题?
第二部分。
这是一种算法,它应该是找到素数的最快算法(https://en.wikipedia.org/wiki/Primality_test)(参见页面中间的伪代码)。这是我的代码,但它不起作用,我进入了一个无限循环。
for i in range(len(euler) - 3):
if (euler[i + 3]== '1' or euler[i + 3]== '3' or euler[i + 3]== '7' or euler[i + 3]== '9'):
number = int(euler[i:i + 3])
for y in range(2, number-1):
while y * y <= number:
if (number % y != 0 or number % (y+1) != 0):
break
我该如何解决?