我想可能有比我这里介绍的更有效的方法,但至少当你意识到你不需要计算分子和分母都在0..1000000范围内的所有因子时,你可以避免双重循环。在
你可以只在分母上循环(从1开始,而不是0),然后增加分子(从0开始),直到它超过给定的分数。此时,将分子减量一次,所以它代表一个潜在的候选解。然后在下一次迭代中,当分母大一倍时,不必再次将分子重置为0,它可以在它离开的地方继续,因为可以确定新的分数将小于给定的分数:这是您真正获得时间的地方。在
这意味着您可以使用线性方法而不是二次法:def get_prev_fraction(num, denom, limit = 1000000):
bestnum = 0
bestdenom = 1
a = 0
for b in range(1,limit+1):
while a*denom < b*num:
a += 1
a -= 1
if a*bestdenom > b*bestnum:
bestnum, bestdenom = a, b
return bestnum, bestdenom
num, denom = get_prev_fraction(2, 7)
print("the fraction which is to the left of {}/{}={} is {}/{}={}".format(2, 7, 2.0/7, num, denom, num/denom))
该输出:the fraction which is to the left of 2/7=0.2857142857142857 is 285713/999996=0.28571414285657143
注意,引用中提到了d≤1000000,所以您需要做range(1000001)来包含这个限制。在