算法系列之一 ：Google方程式

算法系列之一 ： Google方程式

 

    有一个字符组成的等式：WWWDOT - GOOGLE = DOTCOM，每个字符代表一个0－9之间的数字，WWWDOT、GOOGLE和DOTCOM都是合法的数字，不能以0开头。请找出一组字符和数字的对应关系，使它们互相替换，并且替换后的数字能够满足等式。这个字符等式是Google公司能力倾向测试实验室的一道题目，这种题目主要考察人的逻辑推导能力和短期记忆能力，通常棋下的好的人解决这类问题会更得心应手一些（飞行棋例外）。此类型的题目有很多变种，各种编程比赛中常常能见到它们的身影。比如2005年的GOOGLE中国编程挑战赛第二轮淘汰赛有一道名为“SecretSum”的500分的竞赛题，与本题如出一辙，只不过字母都是三个，而且用的是加法计算。现在言归正传，先看看如何分析这个问题。

 

以人的思维方式分析问题

 

将横式改成竖式可能更直观一些：

根据以上竖式减法，从左向右依次可以得到6个算式，分别是：

 W – G = D                     （算式 1）

W – O = O                     （算式 2）

W – O = T                     （算式 3）

D – G = C                     （算式 4）

O – L = O                      （算式 5）

T – E = M                      （算式 6）

 根据以上6个算式可以分析出两个关键信息：一个是W要足够大，因为考虑到它可能被借位的情况还要等于G和D的和；另一个则是本问题的突破口，就是算式2和算式3两次出现的W – O计算。现在分析算式2和算式3，根据是否需要借位，算式2和算式3一共有四种借位组合结果，下面分别对这四种借位组合结果进行分析。

  

1.      W – O = T不需要借位，W – O = O也不需要借位

 

     由于W – O = T和W – O = O都不需要借位，则可由算式2变形得到算式1.1：

 W = 2O                      （算式1.1）

将算式1.1带入算式3，又可以得到算式1.2:

 O = T                        （算式 1.2）

 

 根据算式1.2，O和T代表的数字是同一个数字，这与题目要求不符，因此，这种借位组合不能得到正确的结果。

 

2.      W – O = T需要借位，W – O = O不需要借位

 

     根据借位情况，对算式2和算式3进行借位修正，得到两个修正算式：

 W – 1 – O = O                （算式2.1）

W + 10 – O = T               （算式2.2）

 由算式2.1变形得到算式2.3：

 W = 2O + 1                  （算式2.3）

 将算式2.3带入算式2.2可以得到算式2.4：

 T = O + 11                    (算式2.4)

 

 对算式2.3分析，由于W是个位数，最大值是9，所以O的取值只能是1－4，但是无论如何，由算式2.4计算出的T都超过9，这与题目要求不符，因此，这种情况也是无解的情况。

 

3.      W – O = T不需要借位，W – O = O需要借位

 

    根据借位情况，对算式2和算式3进行借位修正，得到两个修正算式：

 W + 10 – O = O               （算式3.1）

W – O = T                   （算式3.2）

 由算式3.1变形得到算式3.3：

 W = 2O - 10                 （算式3.3）

 将算式3.3带入算式3.2由可得到算式3.4：

 O – 10 = T                  （算式3.4）

 

 O显然是不能比10大的个位数，因此，这种情况也是无解的情况

 

 4.      W – O = T需要借位，W – O = O也需要借位

 

    根据借位情况，对算式2和算式3进行借位修正，得到两个修正算式：

 W – 1 + 10 – O = O            （算式4.1）

W + 10 – O = T               （算式4.2）

 由算式4.1变形得到算式4.3：

 W = 2O - 9                   （算式4.3）

 将算式4.3代入算式4.2得到算式4.4：

 O + 1 = T                    （算式4.4）

 

由于W不能小于0，因此，根据算式4.3，O的取值最小为5。根据算式4.4继续分析，因为T不能大于9，因此O的最大值只能取值为8。根据O的取值区间[5,8]，可依次计算出W和T的值，如下表所示：

 

O	W	T
5	1	6
6	3	7
7	5	8
8	7	9

 

已知O、W、T的取值，可以进一步推算其他字符代表的数字，上表中得到了四组目前合法的取值，但是并不是四组取值都能最终推算出正确的结果，本题的答案只有一个，也就是说只有一组O、W、T的取值是正确的，下面就分别进行分析。

 

            O = 5, W = 1, T = 7

 

在这种情况下，考察算式1：W – G = D，W = 1显然无法满足此种情况，更何况算式2： W – O = O还要从它这里借位，因此，这种情况无解。

 

            O = 6, W = 3, T = 7

 

在这种情况下，算式2： W – O = O还要从它这里借位，因此算式1：W – G = D对应的实际情况是2 – G = D，G和D不能同时为1，而且G和D都是第一位数字，不能是0，因此无法满足算式1，这种情况也是无解。

 

            O = 7, W = 5, T = 8

 

在这种情况下，需要考察另另外两个关键算式，分别是算式5和算式6。根据这两个算式是否需要借位进行不同的假设，根据组合，仍然有四种假设，下面分别分析这四种假设：

 

假设一：算式5需要借位，算式6不需要借位。则此时算式5可修正为O + 10 – L = O，推算出L = 10，显然不符合题意，假设一不成立；

 

假设二：算式5需要借位，算式6需要借位，则算式5和算式6应该修正为算式4.3.1和算式4.3.2：

 

O -1 + 10 – L = O                   （算式4.3.1）

T + 10 – E = M                     （算式4.3.2）

 

因为已知T=8，带入4.3.2可得E+M=18，显然对于两个不相同的个位数无法满足这个等式，因此假设二也不成立；

 

假设三：算式5不需要借位，算式6不需要借位，此时根据算式6可知E和M的和是8（T=8），排除E=M=4的情况后，E和M的组合可以是(1,7)、(2,6)和(3,5)，又因为数字5和7分别被W和O使用，因此E和M只能是2或6。再回头来看算式1，因为算式2需要借位，算式1实际相当于G + D = 4，G和D只能取值1和3，若G=1，D=3，则根据算式4计算出C＝2，这与E或M矛盾。若G=3，D=1，则算式4需要借位，这又与算式3的假设矛盾。由此看来，假设三也不能得到正确的结果；

 

假设四：算式5不需要借位，算式6需要借位，此时根据算式5被修正为O – 1 – L = O，这种情况下也是无解的。

 

            O = 8, W = 7, T = 9

 

在这种情况下，根据算式5和算式6是否借位的假设，仍然有四种假设，下面分别分析这四种假设：

 

假设一：算式5需要借位，算式6不需要借位。则此时算式5可修正为O + 10 – L = O，推算出L = 10，显然不符合题意，假设一不成立；

 

假设二：算式5需要借位，算式6需要借位，则算式5和算式6应该修正为算式4.4.1和算式4.4.2：

 

O -1 + 10 – L = O                   （算式4.4.1）

T + 10 – E = M                     （算式4.4.2）

 

因为已知T=9，带入4.4.2可得E+M=19，两个不同的个位数的和不可能大于18，因此假设二也不成立；

 

假设三：算式5不需要借位，算式6不需要借位，此时根据算式6可知E和M的和是9（T=9），E和M的组合可以是(1,8)、(2,7) 、(4,5)和(3,6)，又因为数字8和7分别被O和W使用，因此E和M只能是(4,5)和(3,6)。进一步假设E＝4，M=5（反过来E＝5，M=4是一样的，不影响分析）。再看算式1，因为算式2需要借位，算式1实际相当于G + D = 6，由于M或E是5，所以G和D只能取值2和4。若G=2，D=4，则根据算式4计算出C＝2，这与G＝2矛盾。若G=4，D=2，则算式4需要借位，这又与算式3的假设矛盾，因此E＝4，M=5的情况无解。再次进一步假设E＝3，M=6（反过来E＝6，M=3是一样的，不影响分析）。同样再看算式1，G和D的值可取是(2,4)和(1,5)，G和D取值1和2的情况刚刚分析过无解，因此G和D的取值只能是1和5，前面分析过，算式4没有借位，也就是说要保证D > G，因此，D＝5，G=1，根据算式4计算出C＝4，这样就得到了一组解：O = 8，W = 7，T = 9，D = 5，L = 0， G = 1，C = 4，E = 3/6，M = 6/3。最终的等式是：

 

777589 - 188103 = 589486

或

777589 - 188106 = 589483

 

假设四：算式5不需要借位，算式6需要借位，此时根据算式5被修正为O – 1 – L = O，这种情况下也是无解的。

 

完整的分析过程结束，得到了一组答案，事实上通过计算机穷举算法也只能得到这一组结果，下面就看看如何用计算机算法求解本题的答案。

 

用计算机穷举所有的解

 

    以上是用人的思维方式的解题过程，如果方法正确，加上运气好（三次假设都是正确的，避免在错误分支上浪费时间），两分钟内就可得到结果。但是考虑到更通用的情况，字母数字没有规律，也没有可供分析的入手点和线索，比如：

    AAB – BBC = CCD

这样的问题，该什么方法解决呢？只能“猜想”，用穷举的方法试探每一种猜想，对每个字母和数字穷举所有可能的组合，直到得到正确的结果。当然，这样的力气活交给计算机做是最合适不过了。

 

1.      建立数学模型

 

    要想让计算机解决问题，就要让计算机能够理解题目，这就需要建立一个计算机能够识别、处理的数学模型，首先要解决的问题就是建立字母和数字的映射关系的数学模型。本题的数学模型很简单，就是一个字母二元组：{char, number}。考察等式：

 

WWWDOT - GOOGLE = DOTCOM

 

共出现了9个不同的字母：W、D、O、T、G、L、E、C和M，因此，最终的解应该是9个字母对应的字母二元组向量：[ {'W', 7}, {'D', 5}, {'O', 8}, {'T', 9}, {'G', 1}, {'L', 0}, {'E', 3}, {'C', 4}, {'M', 6} ]。穷举算法就是对这个字母二元组向量中每个字母二元组的number元素进行穷举，number的穷举范围就是0－9共10个数字，当然，根据题目要求，有一些字符不能为0，比如W、G和D。排列组合问题的穷举多使用多重循环，看样子这个穷举算法应该是9重循环了，在每层循环中对一个字母进行从0到9遍历。问题是，必须这样吗，对于更通用的情况，不是9个字母的问题怎么办？首先思考一下是否每次都要遍历0－9。题目要求每个字母代表一个数字，而且不重复，很显然，对每个字母进行的并不是排列，而是某种形式的组合，举个例子，就是如果W字母占用了数字7，那么其它字母就肯定不是7，所以对D字母遍历是就可以跳过7。进一步，假设某次遍历的字母二元组向量中除M字母外其它8个字母已经有对应的数字了，比如：

 

[ {'W', 7}, {'D', 5}, {'O', 8}, {'T', 9}, {'G', 1}, {'L', 0}, {'E', 3}, {'C', 4}, {'M', ?} ] （序列－1）

 

那么M的可选范围就只有2和6，显然没必要使用9重循环。

 

现在换一种想法，对9个二元组的向量进行遍历，可以分解为两个步骤，首先确定第一个二元组的值，然后对剩下的8个二元组进行遍历。显然这是一种递归的思想（分治），算法很简单，但是要对10个数字的使用情况进行标识，对剩下的二元组进行遍历时只使用没有占用标识的数字。因此还需要一个标识数字占用情况的数字二元组定义，这个二元组可以这样定义：{number, using}，0－9共有10个数字，因此需要维护一个长度为10的数字二元组向量。数字二元组向量的初始值是：

 

[{0, false}, {1, false},{2, false},{3, false},{4, false},{5, false},{6, false},{7, false},{8, false},{9, false}] （序列－2）

 

每进行一重递归就有一个数字的using标志被置为true，当字母二元组向量得到（序列－1）的结果时，对应的数字二元组向量的值应该是：

 

[{0, true}, {1, true},{2, false},{3, true},{4, true},{5, true},{6, false},{7, true},{8, true},{9, true}] （序列－3）

 

此时遍历这个数字二元组向量就可以知道M字母的可选值只能是2或6。

 

穷举遍历的结束条件是每层递归中遍历完所有using标志是false的数字，最外一层遍历完所有using标志是false的数字就结束了算法。

 

根据题目要求，开始位置的数字不能是0，也就是W、G和D这三个字母不能是0，这是一个“剪枝”条件，要利用起来，因此，对字母二元组进行扩充成字母三元组，添加一个leading标志：{char, number, leading}。下面就是这个数学模型的C语言定义：

4 typedef struct     5 {     6   char c;     7   int value;     8     bool leading;     9 }CharItem;    10     11 typedef struct    12 {    13     bool used;    14   int value;    15 }CharValue;

 根据此数学模型初始化字母三元组和数字二元组向量：

29     30     CharItem char_item[max_char_count] =    31         {    32           { 'W', -1, true  }, { 'D', -1, true  }, { 'O', -1, false },    33           { 'T', -1, false }, { 'G', -1, true  }, { 'L', -1, false },    34           { 'E', -1, false }, { 'C', -1, false }, { 'M', -1, false }    35         };    36     37     CharValue char_val[max_number_count] =    38         {    39           {false, 0}, {false, 1}, {false, 2}, {false, 3},    40           {false, 4}, {false, 5}, {false, 6}, {false, 7},    41           {false, 8}, {false, 9}    42         };    43

2.      穷举算法

     建立数学模型，其实就是为了让计算机理解题目并处理相关的数据，算法就是告诉计算机如何使用这些模型中的数据。本文介绍的是穷举算法，算法的核心其实前面已经提到了，就是穷举所有的字母和数字的组合，对每种组合进行合法性判断，如果是合法的组合，就输出结果。

整个算法的核心是SearchingResult()函数，其实这个函数非常简单：

54     55 void SearchingResult(CharItem ci[max_char_count],    56                      CharValue cv[max_number_count],    57                      int index, CharListReadyFuncPtr callback)    58 {    59     if(index == max_char_count)    60     {    61         callback(ci);    62         return;    63     }    64     65     for(int i = 0; i < max_number_count; ++i)    66     {    67         if(IsValueValid(ci[index], cv[i]))    68         {    69             cv[i].used = true;/*set used sign*/    70             ci[index].value = cv[i].value;    71             SearchingResult(ci, cv, index + 1, callback);    72             cv[i].used = false;/*clear used sign*/    73         }    74     }    75 }    76

 SearchingResult()函数有四个参数，ci就是存储遍历结果的字母三元组向量，cv是存储遍历过程中数字占用情况的数字二元组向量，index是当前处理的字母三元组在字母三元组向量中的位置索引，0表示第一个字母三元组。callback是一个回调函数，当ci中所有三元组都分配了数字，就调用callback对这组解进行判断，如果满足算式就输出结果。SearchingResult()函数的代码分两部分，前一部分是结束条件判断和结果输出，后一部分是算法的关键。算法就是遍历cv中的所有数字二元组，对于每一个可用的数字（当前没有被占用，并且满足第一个数字不是0的要求），首先设置占用标志，然后将当前字母三元组的值与这个数字的值绑定，最后递归处理下一个字母三元组。

SearchingResult()函数是一个通用过程，负责字母和数字的组合，回调函数(callback)负责根据题目要求对SearchingResult()函数得到的字母和数字的组合进行筛选，只输出正确的组合。对于本题，回调函数可以这样实现：

6      7 void OnCharListReady(CharItem ci[max_char_count])     8 {     9     char *minuend    = "WWWDOT";    10     char *subtrahend = "GOOGLE";    11     char *diff       = "DOTCOM";    12     13     int m = MakeIntegerValue(ci, minuend);    14     int s = MakeIntegerValue(ci, subtrahend);    15     int d = MakeIntegerValue(ci, diff);    16     if((m - s) == d)    17     {    18         std::cout << m << " - " << s << " = " << d << std::endl;    19     }    20 }    21

 

3.      结果验证

 根据char_item和char_val的初始数据，求解本题的Google方程式：

 

SearchingResult(char_item, char_val, 0, OnCharListReady);

 

穷举算法可以得到两个结果（M和E可以互换）：

 

777589 - 188103 = 589486

777589 - 188106 = 589483

 

由于算法具有通用性，对于前文例子中的等式：

 

AAB – BBC = CCD

 

只需要构造新的字母三元组向量，并修改回调函数的过滤数据即可。新的字母三元组可按照如下方式构造：

33     34     CharItem char_item[max_char_count] = { {'A', -1, true}, {'B', -1, true},{'C', -1, true},    35                                            {'D', -1, false} };    36

 回调函数与前文的OnCharListReady()函数类似，此处不再列出。根据新的字符三元组和回调函数运行算法，可以得到13组结果：

 

443 - 331 = 112

553 - 332 = 221

554 - 441 = 113

665 - 551 = 114

774 - 443 = 331

775 - 552 = 223

776 - 661 = 115

885 - 553 = 332

886 - 662 = 224

887 - 771 = 116

995 - 554 = 441

997 - 772 = 225

998 - 881 = 117

 

对于加法、乘法和除法算式，同样只要使用不用的回调函数进行结果判断即可，不需要修改SearchingResult()函数，例如加法算式：

 

ABC + ABC = BCE

 

可以得到5组结果：

124 + 124 = 248

125 + 125 = 250

249 + 249 = 498

374 + 374 = 748

375 + 375 = 750