python-2.7 - 对于马尔可夫链蒙特卡罗的循环问题

Question

所以这是我的下一个问题。我试图遍历并找出 State_Space 中有多少条目的第 25 个条目为 1，但它一直告诉我答案为 0。这是代码。

import random
import numpy

M=numpy.zeros((52),dtype=int)


z=0
State_Space=[]

for i in range(1,100500):
    x=random.randint(1,50)


    T=M.copy()
    if T[x]==1:
        T[x]=0
    if T[x]==0:
        T[x]=1


    if not any(numpy.array_equal(T, X) for X in State_Space):
        if T[x+1]==0 and T[x-1]==0:
            State_Space.append(T)
            if T[25]==1:
                z=z+1
            M=T


    else:
         if T[x+1]==0 and T[x-1]==0:
            M=T


print z
print len(State_Space)

输出：

0, 22