Rosalind第28题——ros_bio28

如果第一次阅读，请查看写在前面

import math#读取文件
with open('../examples/ros_bio28_PROB.txt') as f:file = f.readlines()dna = file[0].rstrip()
probability = file[1]
probability = probability.split(' ')#根据提供的各个概率计算序列相加的对数
all_lo = []
for i in probability:i = float(i)gpro = i/2apro = (1-i)/2loga = 0for j in dna:if j == 'A' or j == 'T':loga += math.log(apro, 10)else:loga += math.log(gpro, 10)all_lo.append(loga)#格式化为三位小数
result = ''
for line in all_lo:line = format(line, '.3f')result += line + ' '
print(result)

Rosalind第28题——ros_bio28_PROB相关推荐

LeetCode集锦（十） - 第28题 Implement StrStr
LeetCode集锦(十) - 第28题 Implement StrStr 问题 Implement strStr().Return the index of the first occurrence ...
【12c】OCP 062近期新出现的考试原题-第28题
第28题.choose one Unified auditing is enabled in your database. The HR_ADMIN and OE_ADMIN roles exist ...
UESTC 28题-补刀
好吧,我又来不务正业了,UESTC-28题链接: http://acm.uestc.edu.cn/#/problem/show/28 这道题没涉及什么算法,自己写了个暴力破解的版本如下,但是运行时间太 ...
史上最强Dubbo面试28题答案详解：核心功能+服务治理+架构设计等
1.Dubbo是什么? Dubbo 是一个分布式.高性能.透明化的 RPC 服务框架,提供服务自动注册.自动发现等高效服务治理方案, 可以和 Spring 框架无缝集成. RPC 指的是远程调用协议, ...
Rosalind第五题：计算GC内容
问题 DNA字符串的GC含量由字符串中" C"或" G"的符号百分比给出.例如," AGCTATAG"的GC含量为37.5%.请注意,任何D ...
LeetCode第28题实现strStr()之KMP算法（C++）【代码已提交成功】
目录初步思路朴素匹配算法 KMP算法 NEXT数组利用NEXT数组改进朴素匹配算法初步思路这是一道难度为简单的题,所以不熟悉的话可能第一反应就是朴素匹配的算法.但因为考研的时候学过数据结构, ...
Rosalind第四题：兔子与递归关系
http://rosalind.info/problems/fib/ 问题序列是对象(通常是数字)的有序集合,允许重复.序列可以是有限的也可以是无限的.两个例子是有限奇数数列和无限奇数数列.我们用这 ...
[剑指offer]面试题第[28]题[Leedcode][JAVA][第101题][对称二叉树][队列][递归]
[问题描述][简单] 给定一个二叉树,检查它是否是镜像对称的.例如,二叉树 [1,2,2,3,4,4,3] 是对称的.1/ \2 2/ \ / \ 3 4 4 3 [解答思路] 1. 递归时间复杂度 ...
程序员面试题100题第28题——全排列
题目: 输入一字符串,输出字符串的所有排列分析: 我们以三个字符abc为例.首先我们固定第一个字符a,求后面两个字符bc的全排列.当bc的排列求好之后,这就是一部分解: 现在我们把第一个字符a和后面 ...

Rosalind第28题——ros_bio28_PROB

Rosalind第28题——ros_bio28_PROB相关推荐

最新文章

热门文章