Question

这是我偶然发现正确答案的问题之一，但我不明白为什么它是正确的，维基百科没有帮助。对于Rosalind，我写了一个简单的脚本，用于从蛋白质字符串中获取所有可能的RNA序列的数量（模数为1,000,000）。我知道这不是最有效的代码（部分原因是因为我回收了以前所做的事情中的比特），但这里是：

protein = """<large protein string>"""
protein = ''.join(protein.split('\n'))

translate = {'UUU' : 'F','CUU' : 'L','AUU' : 'I','GUU' : 'V','UUC' : 'F','CUC' : 'L','AUC' : 'I','GUC' : 'V','UUA' : 'L','CUA' : 'L','AUA' : 'I','GUA' : 'V','UUG' : 'L','CUG' : 'L','AUG' : 'M','GUG' : 'V','UCU' : 'S','CCU' : 'P','ACU' : 'T','GCU' : 'A','UCC' : 'S','CCC' : 'P','ACC' : 'T','GCC' : 'A','UCA' : 'S','CCA' : 'P','ACA' : 'T','GCA' : 'A','UCG' : 'S','CCG' : 'P','ACG' : 'T','GCG' : 'A','UAU' : 'Y','CAU' : 'H','AAU' : 'N','GAU' : 'D','UAC' : 'Y','CAC' : 'H','AAC' : 'N','GAC' : 'D','UAA' : 'Stop','CAA' : 'Q','AAA' : 'K','GAA' : 'E','UAG' : 'Stop','CAG' : 'Q','AAG' : 'K','GAG' : 'E','UGU' : 'C','CGU' : 'R','AGU' : 'S','GGU' : 'G','UGC' : 'C','CGC' : 'R','AGC' : 'S',
'GGC' : 'G','UGA' : 'Stop','CGA' : 'R','AGA' : 'R','GGA' : 'G','UGG' : 'W','CGG' : 'R','AGG' : 'R','GGG' : 'G',
}
aminos = translate.values()
sample = [l for l in protein] + ['Stop']

score = []
for s in sample:
    c = aminos.count(s)
    score.append(c)

import math
result = reduce(lambda x, y: x*y, score) % 1000000
print result

这计算RNA序列的总数并取最终结果的模数（或者我认为）。在我决定尝试之前，我得到了错误的答案两次：

import math
result = reduce(lambda x, y: x*y % 1000000, score)
print result

这显然产生了正确的答案。为什么必须在每个x * y执行模数？我不理解模数或我不理解Python吗？

编辑：对不起，错字。

Answer 1

之间的区别

reduce(lambda x, y: x*y, score) % 1000000

和

reduce(lambda x, y: x*y % 1000000, score)

第一个必须与long一起使用score中的所有值的乘积，而第二个将使用不大于max(score) * 999999的值。

任意大整数不能存储在有限的内存中，也不能在恒定时间内计算它们的产品，所以你更有可能达到OutOfMemoryError或者花很长时间使用第一个选项

为什么必须在每次迭代期间执行模数？

1 个答案: