Расстояние Хэмминга между двумя двоичными строками не работает

На этом сайте я нашел интересный алгоритм для расчета расстояния Хэмминга:

def hamming2(x,y):
    """Calculate the Hamming distance between two bit strings"""
    assert len(x) == len(y)
    count,z = 0,x^y
    while z:
        count += 1
        z &= z-1 # magic!
    return count

Дело в том, что этот алгоритм работает только с битовыми строками, и я пытаюсь сравнить две строки, которые являются двоичными, но они в строковом формате, как

'100010'
'101000'

Как я могу заставить их работать с этим алгоритмом?

4 ответа

Решение

Реализуйте это:

def hamming2(s1, s2):
    """Calculate the Hamming distance between two bit strings"""
    assert len(s1) == len(s2)
    return sum(c1 != c2 for c1, c2 in zip(s1, s2))

И проверить это:

assert hamming2("1010", "1111") == 2
assert hamming2("1111", "0000") == 4
assert hamming2("1111", "1111") == 0

Это то, что я использую для расчета расстояния Хэмминга.
Он считает количество различий между строками одинаковой длины.

def hamdist(str1, str2):
    diffs = 0
    for ch1, ch2 in zip(str1, str2):
        if ch1 != ch2:
            diffs += 1
    return diffs

Если мы придерживаемся исходного алгоритма, нам нужно преобразовать строки в целые числа, чтобы иметь возможность использовать побитовые операторы.

def hamming2(x_str, y_str):
    """Calculate the Hamming distance between two bit strings"""
    assert len(x_str) == len(y_str)
    x, y = int(x_str, 2), int(y_str, 2)  # '2' specifies we are reading a binary number
    count, z = 0, x ^ y
    while z:
        count += 1
        z &= z - 1  # magic!
    return count

Тогда мы можем назвать это следующим образом:

print(hamming2('100010', '101000'))

Хотя этот алгоритм хорош как новинка, преобразование в строку, скорее всего, сводит на нет любое преимущество в скорости, которое он может иметь. Ответ @dlask опубликован гораздо более кратко.

Я думаю, это хорошо объясняет The Hamming distance

def hammingDist(s1, s2):
    bytesS1=bytes(s1, encoding="ascii")
    bytesS2=bytes(s2, encoding="ascii")
    diff=0;
    for i in range(min(len(bytesS1),len(bytesS2))):
        if(bytesS1[i]^bytesS2[i]!=0):
            diff+=1;
    return(diff)
Другие вопросы по тегам