模拟中的无限循环

Question

我从python开始..我在下面写的详细信息..当我尝试在其内部调用函数时，它进入无限循环并给我一个错误..这种递归不是吗允许吗？

下面的发布代码.. 感谢您的所有帮助:)

该程序假设我们有 100 名乘客登机。假设如果第一个丢失了登机牌，他会随机找到一个座位并坐在那里。然后其他进入的乘客如果空着就坐在他们的位置上，如果有人的话，就坐在其他随机座位上。 最终目的是找出最后一位乘客不会坐在他/她自己座位上的概率。我还没有添加循环部分 将使其成为适当的模拟。上面的问题实际上是一个概率之谜。我正在尝试验证答案，因为我并没有真正遵循推理。

import random
from numpy import zeros

rand = zeros((100,3))
# The rows are : Passenger number , The seat he is occupying and if his designated     seat is occupied. I am assuming that the passengers have seats which are same as the order in which they enter. so the 1st passenger enter has a designated seat number 1, 2nd to enter has no. 2 etc.

def cio(r):  # Says if the seat is occupied ( 1 if occupied, 0 if not)
    if rand[r][2]==1:
        return 1
    if rand[r][2]==0:
        return 0

def assign(ini,mov):    # The first is passenger no. and the second is the final seat he gets. So I keep on chaning the mov variable if the seat that he randomly picked was occupied too. 
    if cio(rand[mov][2])== 0 :
        rand[mov][2] = 1
        rand[mov][1] = ini
    elif cio(rand[mov][2])== 1 :
        mov2 = random.randint(0,99)
 #       print(mov2)            Was used to debug.. didn't really help
        assign(ini,mov2)        # I get the error pointing to this line :(

# Defining the first passenger's stats.
rand[0][0] = 1
rand[0][1] = random.randint(1,100)
m = rand[0][1]
rand[m][2]= 1

for x in range(99):
    rand[x+1][0] = x + 2

for x in range(99):
    assign(x+1,x+1)

if rand[99][0]==rand[99][1] :
    print(1);
else :
    print(0);

如果你们都遇到同样的错误，请告诉我。如果我违反任何规则，因为这是我发布的第一个问题。对不起，如果看起来太长了。

本该如此…… 在这种情况下，代码确实可以使用以下模块正常工作：

def assign(ini,mov):
if cio(mov)== 0 :     """Changed here"""
    rand[mov][2] = 1
    rand[mov][1] = ini
elif cio(mov)== 1 :    """And here"""
    mov2 = random.randint(0,99)
    assign(ini,mov2)  

我在 Windows 7 上使用 Python 2.6.6，使用的是 Python 的 Enthought Academic Version 的软件。 http://www.enthought.com/products/getepd.php

这个谜题的答案也是 0.5，这实际上是我通过运行 10000 次得到的（几乎）。

我没有在这里看到它，但它必须在线提供.. http://www.brightbubble.net/2010/07/10/100-passengers-and-plane-seats/

Answer 1

虽然允许递归，但它不是您最好的首选。

Python 对递归函数实施上限。您的循环似乎超出了上限。

你真的想要某种 while 循环分配。

def assign(ini,mov):    
   """The first is passenger no. and the second is the final seat he gets. So I keep on chaning the mov variable if the seat that he randomly picked was occupied too. 
   """
   while cio(rand[mov][2])== 1:
      mov = random.randint(0,99)

   assert cio(rand[mov][2])== 0
   rand[mov][2] = 1
   rand[mov][1] = ini

这可能是你想要做的更多。

注意对您的 cmets 的更改。在 def 之后的三引号字符串。

Answer 2

您也许可以使用动态规划找到确切的解决方案 http://en.wikipedia.org/wiki/Dynamic_programming 为此，您需要在递归函数中添加记忆： What is memoization and how can I use it in Python?

如果你只是想使用随机数模拟来估计概率，那么我建议你在概率变得非常小时时在一定深度后打破你的递归函数，因为这只会改变一些较小的小数位（大多数可能......您可能希望在更改深度时绘制结果变化）。

要测量深度，您可以在参数中添加一个整数： f（深度）： 如果深度>10： 还东西 否则：f(深度+1)

默认情况下允许的最大递归深度为 1000，尽管您可以更改此设置，但在得到答案之前您只会耗尽内存