博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
python的数学函数(1)-python组合函数模块itertools
阅读量:4553 次
发布时间:2019-06-08

本文共 1264 字,大约阅读时间需要 4 分钟。

要解决的问题:

输出n个 ['A','T','C','G'] 所有的排列组合。
比如n=2 时,输出为 AA,AT,AC,AG,TA,TT,TC,TG,………………
n=3时,输出为 AAA,AAT,AAC,AAG,ATA,ATT,ATC,ATG,……………………
解法1 :传统的map reduce函数
map(func,seq1[,seq2...]):将函数func作用于给定序列的每个元素,并用一个列表来提供返回值;如果func为None,func表现为身份函数,返回一个含有每个序列中元素集合的n个元组的列表。
reduce(func,seq[,init]):func为二元函数,将func作用于seq序列的元素,每次携带一对(先前的结果以及下一个序列的元素),连续的将现有的结果和下一个值作用在获得的随后的结果上,最后减少我们的序列为一个单一的返回值:如果初始值init给定,第一个比较会是init和第一个序列元素而不是序列的头两个元素。
filter(bool_func,seq):此函数的功能相当于过滤器。调用一个布尔函数bool_func来迭代遍历每个seq中的元素;返回一个使bool_seq返回值为true的元素的序列。
用map reduce解此题:

def sequence(n): base=['A','T','C','G'] print reduce(lambda x,y:[(a+b) for a in x for b in y],map(lambda x:[x]*n,[base])[0]) print sequence(4)

这是发帖人自己的一个实现
[(a+b) for a in x for b in y] 产生一个笛卡尔积
map(lambda x:[x]*n,[base])[0] 其实等价于 [base]*4,写法太费解
reduce源码:

def reduce(bin_func,seq,initial=None): lseq = list(seq) if initial is None: res = lseq.pop(0) else: res = initial for eachItem in lseq: res = bin_func(res,eachItem) return res

解法2 :itertools组合函数
import itertools
# python 3 or use __future__
print([''.join(x) for x in itertools.product('ATCG', repeat=4)]
product(A, B) 等价于 ((x,y) for x in A for y in B).
product(A, repeat=4) 等价于 product(A, A, A, A).

转载于:https://www.cnblogs.com/shenguanpu/archive/2011/12/23/2299742.html

你可能感兴趣的文章
禁止键盘上的刷新键F5等
查看>>
SAP中对于获取订单的状态
查看>>
oracle PL/SQL块
查看>>
CentOS7集群环境Elastic配置
查看>>
【EXCEL】指定の項目の内容一覧を表示
查看>>
Head first java chapter 4 对象的行为
查看>>
luogu_4503【题解】企鹅QQ 哈希
查看>>
linux 面试
查看>>
Linux:Gentoo系统的安装笔记(三)
查看>>
打开IE窗口并最大化显示
查看>>
Windows API SendMessage 和 PostMessage 内部实现
查看>>
服务器发送邮件出现Could not connect to SMTP host错误 解决办法
查看>>
sklearn.preprocessing.LabelBinarizer
查看>>
C teaching
查看>>
分隔指定内容,提取章节数
查看>>
this point
查看>>
leetcode 30 Substring with Concatenation of All Words
查看>>
验证登录信息是否合法
查看>>
线程池
查看>>
git版本控制器的基本使用
查看>>