Полиалфавитные подстановки, не получается расшифровать сообщение

#!/usr/bin/python
# -*- coding: utf-8 -*-
 
# Индекс совпадений
def get_Ic(a):
    return float(sum([a.count(c)*(a.count(c)-1) for c in set(list(a))])) / (len(a) * (len(a) - 1))
 
# Определение длины ключа (только для вывода на экран) 
def find_len(a):
    for klen in xrange(1, 20):
        IC = 0
        for i in xrange(klen):
            b = [ a[j] for j in xrange(i, len(a), klen)]
            IC += get_Ic(b)
        IC /= klen
        print "klen:", klen, ":", IC
 
# Индекс взаимных совпадений
def get_MIc(a, b):
    return float(sum([a.count(c)*b.count(c) for c in set(list(a) + list(b))])) / (len(a) * len(b))
 
# Шифрование списка шифром Цезаря
def caesar_enc(a, k, m):
    b = [(x + k) % m for x in a]
    return b
 
# Расшифрование
def caesar_dec(a, k, m):
    b = [(x - k) % m for x in a]
    return b
 
# Определение разности ключей для двух списков (колонок текста)
def get_shift(a, b, m):
    MIC = [get_MIc(a, caesar_enc(b, i, m)) for i in xrange(m)]
    print max(MIC)
    return -MIC.index(max(MIC)) % m 
 
# Чтение шифротекста
f = open("cipher.txt", "r")
xc = f.read()
f.close()
uc = xc.decode('cp1251').strip().upper()
 
# Алфавит
abc = u"АБВГДЕЖЗИЙКЛМНОПРСТУФХЦЧШЩЬЫЭЮЯ "
 
# Кодирование символов
cipher = []
for ch in uc:
    cipher.append(abc.index(ch))
 
# Выводим индексы совпадений для различных длин ключа
find_len(cipher)
 
# Разбиение текста на столбцы в соответствии с выбранной длиной ключа
klen = 6
cols = [[] for i in xrange(klen)]
i = 0
for ch in cipher:
    cols[i].append(ch)
    i += 1
    if i >= klen: i = 0
 
# С помощью индекса взаимных совпадений определяем ключ
# Ключ получается зашифрованным шифром цезаря
key = [get_shift(cols[0], cols[i], len(abc)) for i in xrange(klen)]
print key
 
# Логарифмы частостей наиболее часто встречаемых символов (в произвольном масштабе)
rate = len(abc)*[0]
ch = [u' ', u'О', u'Е', u'А', u'И', u'Н', u'Т', u'С', u'Р', u'В', u'Л']
r  = [1.92, 1.72, 1.61, 1.58, 1.55, 1.50, 1.48, 1.42, 1.35, 1.33, 1.32]
for i in xrange(len(ch)): rate[abc.index(ch[i])] = r[i]
 
# Перебираем все варианты ключа (по количеству символов алфавита)
R = []
for k in xrange(len(abc)):
    work_key = caesar_dec(key, k, len(abc))
    print "KEY:", ''.join([abc[ch] for ch in work_key])
    plain = [(cipher[i] - work_key[i % len(work_key)]) % len(abc) for i in xrange(len(cipher))]
    R.append(sum([rate[p] for p in plain]))
    print ''.join([abc[p] for p in plain])
    print
 
# Выводим наиболее вероятный ключ и текст
k = R.index(max(R))
best_key = caesar_dec(key, k, len(abc))
print "========="
print "BEST KEY:", ''.join([abc[ch] for ch in best_key])
plain = [(cipher[i] - best_key[i % len(best_key)]) % len(abc) for i in xrange(len(cipher))]
print ''.join([abc[p] for p in plain])