python解答欧拉计划第59题

原题如下：

大概翻译了一下，这个题目是关于异或加密的，前四段是对异或加密的简单介绍，第五段是问题描述，cipher.txt里面是使用异或加密后的密文对应的ascii码，已知明文是常见的英文单词，密钥是由三位小写英文字母构成，要我们找到密钥并求出明文对应ascii码值的和。

要解决这个问题，首先得了解一下什么是异或加密。我们知道，字符在计算机中有对应的ascii码值，异或加密就是将明文字符对应的ascii码值与密钥字符对应的ascii码值相异或把得到的ascii码值对应的字符作为密文，从而实现加密。对于明文长度长于密钥的情况，则循环使用密钥对明文进行加密，如明文为I am happy，密钥为key，加密时应将明文与keykeykeyk，需注意这里空格也有对应的ascii码，加密的时候如果是对整个明文直接进行加密，空格也会被加密（这里是个坑，这个题目就有涉及）。

异或运算具有可逆性，如65 XOR 42 = 107 则 107 XOR 42 = 65，因此对于异或运算a XOR b = c，我们只需知道abc中的任意两个，将这两个相异或便可求出第三个，利用这个性质我们在知道明文和密文的情况下便可求得密钥，知道密钥和密文的情况下可求得明文，这也是这个题目的解题思路——先求得密钥，再求出明文。

一开始我看到这个题，我是想先构造所有三位小写字母可能的组合，一共有26*26*26=17576种组合，然后尝试把每一种组合当作密钥去和密文进行异或，观察哪一种可以得到正确的英文单词但这样子无疑计算量很大。然后我又想到了字母频率攻击，因为异或加密使用密钥加密的过程和维吉尼亚密码很像，所以密文中肯定也符合字母频率规律，在英文中字母e出现的频率最高，我的想法是统计出使用同一个密钥字母加密的密文中出现频率最高的数字，这个数字对应的明文应该是字母e的ascii码，将这两个数字相异或便可得到对应密钥，因为密钥长度为三，所以我将密文的第1 4 7 10...划为一组，第2 5 8 11...划为一组，第3 6 9 12..划为一组，分别找到每组中出现次数最高的数字，将这个数字和e的ascii码101相异或，得到密钥，源代码如下：

输出的结果是71 79 68，然而。。。这三个数字和101异或得到的密钥不是小写字母。。。坑，思路应该没错，问题出在哪里呢？通过查找资料我发现，问题出在异或加密将空格也加密了，题目中的密文又比较短，出现频率最高的不是e应该是空格，所以应该将得到的三个数字分别与空格的ascii码32相异或（坑。。。），尝试了一下果然没问题，顺利得到了密钥，剩下的问题就很简单了，将密钥和密文异或，得到明文然后求和，顺利得到答案，107359。

下面是完整的源代码：

import re
import collections
a = '79,59,12,2,79,35,8,28,20,2,3,68,8,9,68,45,0,12,9,67,68,4,7,5,23,27,1,21,79,85,78,79,85,71,38,10,71,27,12,2,79,6,2,8,13,9,1,13,9,8,68,19,7,1,71,56,11,21,11,68,6,3,22,2,14,0,30,79,1,31,6,23,19,10,0,73,79,44,2,79,19,6,28,68,16,6,16,15,79,35,8,11,72,71,14,10,3,79,12,2,79,19,6,28,68,32,0,0,73,79,86,71,39,1,71,24,5,20,79,13,9,79,16,15,10,68,5,10,3,14,1,10,14,1,3,71,24,13,19,7,68,32,0,0,73,79,87,71,39,1,71,12,22,2,14,16,2,11,68,2,25,1,21,22,16,15,6,10,0,79,16,15,10,22,2,79,13,20,65,68,41,0,16,15,6,10,0,79,1,31,6,23,19,28,68,19,7,5,19,79,12,2,79,0,14,11,10,64,27,68,10,14,15,2,65,68,83,79,40,14,9,1,71,6,16,20,10,8,1,79,19,6,28,68,14,1,68,15,6,9,75,79,5,9,11,68,19,7,13,20,79,8,14,9,1,71,8,13,17,10,23,71,3,13,0,7,16,71,27,11,71,10,18,2,29,29,8,1,1,73,79,81,71,59,12,2,79,8,14,8,12,19,79,23,15,6,10,2,28,68,19,7,22,8,26,3,15,79,16,15,10,68,3,14,22,12,1,1,20,28,72,71,14,10,3,79,16,15,10,68,3,14,22,12,1,1,20,28,68,4,14,10,71,1,1,17,10,22,71,10,28,19,6,10,0,26,13,20,7,68,14,27,74,71,89,68,32,0,0,71,28,1,9,27,68,45,0,12,9,79,16,15,10,68,37,14,20,19,6,23,19,79,83,71,27,11,71,27,1,11,3,68,2,25,1,21,22,11,9,10,68,6,13,11,18,27,68,19,7,1,71,3,13,0,7,16,71,28,11,71,27,12,6,27,68,2,25,1,21,22,11,9,10,68,10,6,3,15,27,68,5,10,8,14,10,18,2,79,6,2,12,5,18,28,1,71,0,2,71,7,13,20,79,16,2,28,16,14,2,11,9,22,74,71,87,68,45,0,12,9,79,12,14,2,23,2,3,2,71,24,5,20,79,10,8,27,68,19,7,1,71,3,13,0,7,16,92,79,12,2,79,19,6,28,68,8,1,8,30,79,5,71,24,13,19,1,1,20,28,68,19,0,68,19,7,1,71,3,13,0,7,16,73,79,93,71,59,12,2,79,11,9,10,68,16,7,11,71,6,23,71,27,12,2,79,16,21,26,1,71,3,13,0,7,16,75,79,19,15,0,68,0,6,18,2,28,68,11,6,3,15,27,68,19,0,68,2,25,1,21,22,11,9,10,72,71,24,5,20,79,3,8,6,10,0,79,16,8,79,7,8,2,1,71,6,10,19,0,68,19,7,1,71,24,11,21,3,0,73,79,85,87,79,38,18,27,68,6,3,16,15,0,17,0,7,68,19,7,1,71,24,11,21,3,0,71,24,5,20,79,9,6,11,1,71,27,12,21,0,17,0,7,68,15,6,9,75,79,16,15,10,68,16,0,22,11,11,68,3,6,0,9,72,16,71,29,1,4,0,3,9,6,30,2,79,12,14,2,68,16,7,1,9,79,12,2,79,7,6,2,1,73,79,85,86,79,33,17,10,10,71,6,10,71,7,13,20,79,11,16,1,68,11,14,10,3,79,5,9,11,68,6,2,11,9,8,68,15,6,23,71,0,19,9,79,20,2,0,20,11,10,72,71,7,1,71,24,5,20,79,10,8,27,68,6,12,7,2,31,16,2,11,74,71,94,86,71,45,17,19,79,16,8,79,5,11,3,68,16,7,11,71,13,1,11,6,1,17,10,0,71,7,13,10,79,5,9,11,68,6,12,7,2,31,16,2,11,68,15,6,9,75,79,12,2,79,3,6,25,1,71,27,12,2,79,22,14,8,12,19,79,16,8,79,6,2,12,11,10,10,68,4,7,13,11,11,22,2,1,68,8,9,68,32,0,0,73,79,85,84,79,48,15,10,29,71,14,22,2,79,22,2,13,11,21,1,69,71,59,12,14,28,68,14,28,68,9,0,16,71,14,68,23,7,29,20,6,7,6,3,68,5,6,22,19,7,68,21,10,23,18,3,16,14,1,3,71,9,22,8,2,68,15,26,9,6,1,68,23,14,23,20,6,11,9,79,11,21,79,20,11,14,10,75,79,16,15,6,23,71,29,1,5,6,22,19,7,68,4,0,9,2,28,68,1,29,11,10,79,35,8,11,74,86,91,68,52,0,68,19,7,1,71,56,11,21,11,68,5,10,7,6,2,1,71,7,17,10,14,10,71,14,10,3,79,8,14,25,1,3,79,12,2,29,1,71,0,10,71,10,5,21,27,12,71,14,9,8,1,3,71,26,23,73,79,44,2,79,19,6,28,68,1,26,8,11,79,11,1,79,17,9,9,5,14,3,13,9,8,68,11,0,18,2,79,5,9,11,68,1,14,13,19,7,2,18,3,10,2,28,23,73,79,37,9,11,68,16,10,68,15,14,18,2,79,23,2,10,10,71,7,13,20,79,3,11,0,22,30,67,68,19,7,1,71,8,8,8,29,29,71,0,2,71,27,12,2,79,11,9,3,29,71,60,11,9,79,11,1,79,16,15,10,68,33,14,16,15,10,22,73'
b = a.split(',')
c1 = b[::3]             #切片取出密文中的第1 4 7 10....个数字，这是使用第一个密钥字母加密的
c2 = b[1::3]            #切片取出密文中的第2 5 8 11....个数字，这是使用第二个密钥字母加密的
c3 = b[2::3]            #切片取出密文中的第3 6 9 12....个数字，这是使用第三个密约字母加密的
                                                                                                #分别找出C1 C2 C3 中出现频率最高的数字，
k = [int(max(c1, key=c1.count))^32,int(max(c2, key=c2.count))^32,int(max(c3, key=c3.count))^32] #并将这个数字与32异或，得到的数字即为密钥，
                                                                                                #并按顺序存入列表k中
sum0 = 0
i = 0
while(i<3):
    sum0 += sum([int(b[j]) ^ k[i] for j in range(i, len(b), 3)])#按顺序使用密钥异或相应的密文，从b中第i个数开始，每隔两个数和k[i]对应的密钥异或一次
                                                                #并将得到的数字求和，
    i += 1
print (sum0)