我一直在阅读一本关于算法的书(重点是Python)的练习,当我尝试解决某个问题时,我会遇到奇怪的行为。练习相当简单,目标是从k
生成随机整数range(n)
,然后通过一系列“是/否”问题找到此数字。只要n
不是太大,下面的代码就可以工作,但是当它得到无限循环时。但这不应该发生,因为在最坏的情况下,step_size
最终会等于1,而position
会爬到k
的值,直到position == k
。
当我同时打印step_size
和position
时,似乎position
在它变得太大时没有更新,所以我猜它与获取值有关近似和截断(这就是为什么我试图明确地将其强制转换,实际上可能没有做任何事情)。代码的某些小部分是不必要的,但我希望它适用于一般n
。关于发生了什么的任何想法(你可能需要增加n
才能得到我正在谈论的行为)?
from random import randrange
from math import ceil
n = 10**15
k = randrange(n)
success = False
step_size = ceil(n / 2.0)
position = long(ceil(n / 2.0))
print(k)
while success == False:
if position == k:
success = True
print(str(int(position)) + ' is the value of k.')
elif k < position:
step_size = ceil(step_size / 2.0)
position -= step_size
print(step_size, position)
else:
step_size = ceil(step_size / 2.0)
position += step_size
print(step_size, position)
答案 0 :(得分:3)
您需要使用position += long(step_size)
和position -= long(step_size)
if position == k:
success = True
print(str((position)) + ' is the value of k.')
elif k < position:
step_size = ceil(step_size / 2.0)
position -= long(step_size)
print(step_size, position,k)
else:
step_size = ceil(step_size / 2.0)
position += long(step_size)
print(step_size, position)
如果你不position
将是一个不长的浮动。
对于n = 10 ** 200
的更改,我们得到的值为:
11008769984569663730780658772914869218303604494435631537943840906720756507580926318660187453313745419228469341648307070867052432521026422402953409000922062458195678772749579263632191432518529106302726
11008769984569663730780658772914869218303604494435631537943840906720756507580926318660187453313745419228469341648307070867052432521026422402953409000922062458195678772749579263632191432518529106302726 is the value of k.
答案 1 :(得分:2)
你遇到了精确问题。确保您的值始终表示为整数(而不是浮点数)。 Python整数具有任意精度,但浮点数不具有。