考虑以下问题。给定的是长度为n的数组p,表示双射。即元素0假设p是一个置换,它可以用循环符号来表示,例如here。
例如,如果
P = [16, 3, 10, 6, 5, 9, 1, 19, 13, 14, 7, 0, 2, 8, 12, 4, 17, 15, 11, 18]
那么循环符号的结果是
[(16, 17, 15, 4, 5, 9, 14, 12, 2, 10, 7, 19, 18, 11, 0), (3, 6, 1), (13, 8)]
下面是完成此任务的Python方法
def toCycle(p):
covered = cur = 0
perm = []
n = len(p)
done = [0]*n
while covered < n:
while cur < n and done[cur] == -1:
cur+=1
cycle = [p[cur]]
sec = p[p[cur]]
done[p[cur]] = -1
done[cur] = -1
covered+=1
while sec != cycle[0]:
cycle.append(sec)
done[sec] = -1
sec = p[sec]
covered+=1
perm+=[tuple(cycle)]
return perm
该算法显然是在线性时间内运行的(done/p的每一个元素被访问的次数都是固定的),因此没有多少可以渐近地完成的工作。
由于我必须对大量的大排列使用这种方法,我想知道
你能快点吗你对表演有什么建议吗
进步?
最佳答案
def cycling(p, start, done):
while not done[e]:
done[e] = True
e = p[e]
yield e
def toCycle(p):
done = [False]*len(p)
cycles = [tuple(cycling(p, 0, done))]
while not all(done):
start = done.index(False)
cycles.append(tuple(cycling(p, start, done)))
return cycles
以你的例子来说,我的代码比你的快30%。