2

これは、私が正しい答えに出くわした質問の 1 つですが、なぜそれが正しい答えなのか理解できず、ウィキペディアは役に立ちませんでした。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 シーケンスの総数を計算し、最終結果のモジュロを取ります (またはそう思います)。これを試す前に、2回間違った答えを得ました:

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

これにより、明らかに正しい答えが得られました。x*y ごとにモジュロを実行する必要があるのはなぜですか? 私はモジュロを理解していませんか、それとも Python を理解していませんか?

編集:申し訳ありませんが、タイプミス。

4

1 に答える 1