背景:
这个题目是关于RNA二级结构和完美匹配的。题目给出了一个RNA字符串,要求计算在它的键合图中,有多少种可能的完美匹配。完美匹配是指在图中,每个节点都恰好与另一个节点相连,形成{A, U}或{C, G}的碱基对。
https://rosalind.info/problems/pmch/recent/
题目的要求是:
- 给定:一个长度不超过80的RNA字符串,它的’A’和’U’的个数相同,'C’和’G’的个数相同。
- 返回:在它的键合图中,完美匹配的总数。
题目的一个例子是:
- 给定:AGCUAGUCAU
- 返回:12
分析:
一个A能与一个U相连。三个A和U就有3!次搭配
例如:A1U1,A1U2,A1U3.A2U2,A2U3,A3U3
GC同理。
所以我们只要统计这个字符串中,A或U的个数,G或C的个数。再把他们的阶乘,乘起来,就是他们组合的个数。
from method import fasta #导入处理fasta格式的包
name_list,value_list = fasta("D:/pycharm/Rosalind/Data/rosalind_pmch.txt")
seq = value_list[0] #得到序列
def factorial(num): #阶乘函数
if(num == 0):
return 1
if num == 1:
return 1
return (num*factorial(num-1))
Anums = seq.count('A')
Gnums = seq.count('G')
print(factorial(Anums)*factorial(Gnums))