我想可能有比我在这里提出的方法更有效的方法,但是当你意识到你不需要计算分子和分母在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))
这输出:在2/7 = 0.2857142857142857左边的分数是285713/99996 = 0.28571414285657143
请注意,引用提及d≤1,000,000,因此您需要range(1000001)包含该限制。