如何使 for 循环使用 2 个字符串？

Question

我有 2 个字符串要比较。问题是我不知道如何同时对两个切片进行迭代，然后在每个切片上运行比较代码。下面是我尝试通过在 for 循环中使用 for 循环来运行的代码（每个字符串一个) 这没有带来任何错误，但只有第一个字符串被迭代切片并打印到屏幕上。第二次尝试是将 2 个 for 循环合并在一起，但是出现了提供的错误消息，它清楚地解释了失败的原因，但是我不知道解决该错误的方法。

for循环中for循环的代码

for i in range(8, len(seq1)):
    seq1i = seq1[0:i]
    print seq1i

    for i in range(8,len(seq2)):
        seq2i = seq2[0:i]
        print seq2i

        for c1, c2 in zip(seq1i, seq2i):
            if c1 == c2:
                counter += 1
            else:
                break
        print "String share the longest common prefix of length %d: %s" % (counter, seq1[:counter])
        if counter >= 8:
            print " Match found, match length %d" % counter

尝试合并 for 循环：

for i1,i2 in range(8,len(seq1, seq2)):
    for c1, c2 in zip(seq1i, seq2i):
        if c1 == c2:
            counter += 1
        else:
            break
    print "String share the longest common prefix of length %d: %s" % (counter, seq1[:counter])
    if counter >= 8:
        print " Match found, match length %d" % counter

这导致了错误：

TypeError                                 Traceback (most recent call last)
<ipython-input-26-9c4822069792> in <module>()
----> 1 for i1,i2 in range(8,len(seq1, seq2)):
      2     for c1, c2 in zip(seq1i, seq2i):
      3         if c1 == c2:
      4             counter += 1
      5         else:

TypeError: len() takes exactly one argument (2 given)

我认为必须有一种方法可以对字符串进行每次迭代切片，并且在制作它们时可以运行匹配的（如果 c1 == c2）代码。

想要的输出是：

如果 seq1 和 seq2：

seq1 = 'ATCGATCGATCGAAAA'
seq2 = 'ATCGATCGATCCCCCC'

A
A
AT
AT
ATC
ATC
ATCG
ATCG
ATCGA
ATCGA
ATCGAT
ATCGAT
ATCGATC
ATCGATC
ATCGATCG
ATCGATCG
ATCGATCGA
ATCGATCGA
ATCGATCGAT
ATCGATCGAT
ATCGATCGATC
ATCGATCGATC
String share the longest common prefix of length 10:ATCGATCGATC

Answer 1

我确信使用高级切片和/或 itertools 有更好的解决方案，但老实说，我认为一个简单的循环是这里最简单的事情。

for i, c in enumerate(seq1):
    try:
        if c != seq2[i]:
            break
        print seq1[:i + 1] + '\n' + seq1[:i + 1]
    except IndexError:
        break
print 'Longest common prefix %s of length %d' % (seq1[:i], i)

Answer 2

尽管来自os.path 模块，但您应该为此使用os.path.commonprefix。

>> s1 = 'GATTACA'
>> s2 = 'GATTCCC'
>> os.path.commonprefix([s1,s2])
'GATT'

Answer 3

打印中间结果有点低效，但这似乎可行：

def lcp(seq1, seq2):
    ret = []
    for i in xrange(min(len(seq1), len(seq1))):
        if seq1[i] != seq2[i]:
            break
        else:
            ret.append(seq1[i])
            print ''.join(ret)
    print "strings share a longest common prefix of length:", len(ret)

seq1 = 'ATCGATCGATCGAAAA'
seq2 = 'ATCGATCGATCCCCCC'
lcp(seq1, seq2)

输出：

A
AT
ATC
ATCG
ATCGA
ATCGAT
ATCGATC
ATCGATCG
ATCGATCGA
ATCGATCGAT
ATCGATCGATC
strings share a longest common prefix of length: 11

Answer 4

您实际上可以zip 字符串，这似乎比这里的大多数答案更pythonic：

str1 = 'ACTGGAG'
str2 = 'ACTA'

prefix = ''
for a,b in zip(str1,str2):
    if a==b:
        prefix += a
    else:
        break

print prefix # prints ACT