Category: 经济、IT观察与思考

这真的不是什么可以引以为豪的事情....我一直认为我是懂动态规划的，直到这两天重新看到动态规划的代码发现自己看不懂，然后恍然间意识到上次看懂都是7年前的事情了。google了一番搜到自己的blog真的是欲哭无泪，然后痛定思痛，觉得这次把它搞懂，重新写一篇笔记，这样万一若干年之后再回头看这个，至少保证这次的笔记有更多的含金量自己可以看懂。（更惭愧的是，高级宏观的时候天天在手动解动态规划，最多的就是无穷期动态规划，现在居然不怎么记得当年是怎么解的了...）

动态规划的用处还真多。很多例子都是斐波那次数列的，但是其实我感觉这样的例子并没有很明显的感觉。倒是今天看到一个文本排列的例子觉得很有意思，原来latex计算每行放多少词是用动态规划算的。想想word打字的时候也是不时会重排一下，所以大概也是在后面不停的算动态规划的最优结果吧。

动态规划的想法其实并不麻烦，大致就是一种以空间换时间的交换。今天耐着性子把youtube上mit公开课关于动态规划的两节看完了（19节和20节），顺便拿笔在旁边记了一堆笔记。然后找了两个例子用python练习了一番，又看了一下其他人的答案，开车回家的路上顺便又想了一下，这才觉得这次好像是想明白动态规划了。所以在这里记一下。

最短路径的动态规划解法

先来个简单的例子？路径问题好了。这个好像是最经典的动态规划例子了。我这里随便画了一个图（神器在此）。

假设如同上图所示，我希望从A走到D，其中各条路径的长度已经标注在图上。那么最短的路径是哪条呢？

最笨的办法，我们可以把每条路径都列出来，一个一个走一遍呗。这里的可能性就是

• A -> B1 -> C1 -> D： 3＋3+10 ＝ 16
• A -> B1 -> C2 -> D：3+5+2 ＝ 10
• A -> B2 -> C1 -> D：4+6+10 ＝ 20
• A -> B2 -> C2 -> D：4+8+2=14

所以很显然，最短路径是第二条：A -> B1 -> C2 -> D。

那么如果问题再复杂一点呢？

这里我们还是可以继续采用笨办法，只是可能的路径多了一点：

• A -> B1 -> C1 -> D1 -> E： 3＋3+10+6 ＝ 22
• A -> B1 -> C1 -> D2 -> E： 3＋3+4+2 ＝ 12
• A -> B1 -> C2 -> D1 -> E：3+5+2+6 ＝ 16
• A -> B1 -> C2 -> D2 -> E：3+5+6+2＝ 16
• A -> B2 -> C1 -> D1 -> E：4+6+10+6 ＝ 26
• A -> B2 -> C1 -> D2 -> E：4+6+6+2 ＝ 18
• A -> B2 -> C2 -> D1 -> E：4+8+2+6 = 8
• A -> B2 -> C2 -> D2 -> E：4+8+6+2 = 20

所以很显然，第二条胜出。在这个过程中我们一共计算了8条路径、24次加法。有没有发现什么规律呢？我们其实有很多重复的计算：比如C1 -> D1 -> E我们计算了两遍。总结一下：

• 不管我们是怎么走到C1的，从C1到最后终点E的最短路径一定是C1 -> D2 -> E，距离为6。同理，不管我们怎么走到C2的，从C2到E怎么走都是8的距离。
• 这样，继续往前推，不管我们是怎么走到B1的，若是从B1到C1再到E，最短距离就是3+6＝ 9 （C1 -> D1 -> E）；若是从B1到C2再到E，最短距离就是5+8=13，所以从B1到E的最短距离就是9。同理，不管我们怎么到B2的，从B2到E的最短距离是 6+6<8+8，故为12。
• 那么回到最初的起点，A到B1再到E，最短距离就是3+9=12；A到B2再到E，最短距离就是4+12=16。

所以在这个过程中，每一步其实我们只进行局部计算就好了，不需要把各种可能性都列出来。下面是我们在每一步可以排除的路径。

• A -> B1 -> C1 -> D1 -> E：
• A -> B1 -> C1 -> D2 -> E： 3＋3+4+2 ＝ 12
• A -> B1 -> C2 -> D1 -> E：
• A -> B1 -> C2 -> D2 -> E：
• A -> B2 -> C1 -> D1 -> E：
• A -> B2 -> C1 -> D2 -> E：4+6+6+2 ＝ 18
• A -> B2 -> C2 -> D1 -> E：
• A -> B2 -> C2 -> D2 -> E：

这就是动态规划的力量了：在我们这个例子中，我们可以倒推出来，给定某一步的各种可能性、后面的最优走法。然后这样每次前面的只需要对比怎么走下一步、后面的最优路径就已经计算好了。

这大概是我自己对动态规划最直观的理解了：每一步都是相对独立的，所以可以先不管前面的、针对后面的进行优化、然后慢慢往前推。因为只要到了某步、后面的结果其实并不取决于是怎么到当前步的。

斐波那次数列的动态规划解法

然后看一下最著名的斐波那次数列好了，就是那个著名的数兔子数列：第一年没有兔子，第二年一只兔子，第三年1一只兔子，第四年2只兔子，第五年3只兔子，第六年5只兔子...每次都是前两年的兔子的和。

这个例子是绝佳的演示为什么动态规划体现了用空间换时间。比如我们要算第10年几只兔子，然后上面已经算好了直到第6年的兔子。我们先来看一种最笨的算法。括号里面的数字是倒推出来写的。

第十年的兔子 ＝ 第八年的兔子（5＋5+3）＋第九年的兔子（5+5+3+5+3）

第九年的兔子 ＝ 第八年的兔子 （5＋5+3）＋ 第七年的兔子（5+3）

第八年的兔子 ＝ 第六年的兔子 （5）＋ 第七年的兔子 （5+3）

第七年的兔子 ＝ 第六年的兔子（5） ＋ 第五年的兔子（3）

.....

一直算下去的话，为了算第10年的兔子，我们要算7次加法。这个过程中可以看出来，除了第五年和第六年的兔子是已知的之外，我们算了两遍第七年的兔子、两遍第八年的兔子、一遍第九年的兔子，然后才算出来第十年的兔子，显然是重复计算。

那动态规划这里怎么用呢？很简单啊，算过的我们就存起来，然后下一次再问到就不用算了呗；没算过的就现算呗。同样，我们已知的是｛第一年，0｝，｛第二年，1｝，｛第三年，1｝，｛第四年，2｝，｛第五年，3｝，｛第六年，5｝。

所以这个过程就是：

• 为了算第十年的，我需要知道第八年的和第九年的，然后这俩都要算。
• 为了算第九年的，我需要知道第七年的和第八年的，然后这俩都要算。
• 为了算第八年的，我需要知道第七年的和第六年的，然后第七年要算。
• 为了算第七年的，我需要知道第六年的和第五年的，这俩都知道，所以第七年是8。
• 这样，第八年的就是8+5=13。
• 这样，第九年就是8+13 ＝ 21。
• 这样，第十年就是，13+21=34。

于是，每一年我只算了一遍，算好了就存起来了，下次备用就好了。

于是有童鞋说，另外一种办法就是从头到尾算就好了嘛、一个一个往后算、到了第10年停就是了。其实这样从前往后和刚才倒推＋存储是一摸一样的计算过程、每一年只算了一遍。因为有存储的存在、动态规划会极大降低时间复杂度。不过显然最省内存的就是从头往后算了，因为我只需要记住n－1和n－2两年的兔子就可以了，不需要知道再往前的年份的。这又体现了一种相对独立的感觉：给定n－1和n－2，n跟n－3...等等就完全没关系了，想想这不就是类似时间序列中的AR(2)过程嘛！

<不过话说最高效的算法，还是通项公式吧，直接就出结果了。但那个就跟这里没关系了呢。>

强盗问题的动态规划解法

最后再来俩个比较好玩的问题吧。House Robber problem，直接复制一下别人的翻译：

你是一名专业强盗，计划沿着一条街打家劫舍。每间房屋都储存有一定数量的金钱，唯一能阻止你打劫的约束条件就是：由于房屋之间有安全系统相连，如果同一个晚上有两间相邻的房屋被闯入，它们就会自动联络警察，因此不可以打劫相邻的房屋。

给定一列非负整数，代表每间房屋的金钱数，计算出在不惊动警察的前提下一晚上最多可以打劫到的金钱数。

我们先来看一下简单的情况。

如图，这条街上一共有6栋房子，我们假设它们依次有1-6块钱。然后可行的打劫策略有：

• 打劫1
  • 打劫3
    • 打劫5：  1+3+5=9
    • 打劫6： 1+3+6 ＝10
  • 打劫4
    • 打劫6： 1+4+6=11
• 打劫2
  • 打劫4
    • 打劫6： 2+4+6 ＝ 12
  • 打劫5： 2+5 ＝ 7

所以简单计算可知，2，4，6是最佳打劫策略。在这个过程中有没有什么熟悉的感觉？比如，给定打劫4，最优的策略一定是打劫6；给定打劫3、最优的策略已经是打劫6。所以我们可以一步步倒推出来、给定某个点、往后最优策略是什么，然后往前慢慢比较前一步怎么走到这个点就可以了。其实无非就是另外一个最短（长）路径问题。

那大致思路就是：已知第i步最优策略，那么只需要比较从i-2走过来和i-3走过来哪个更优就可以了。而我们又知道，这个过程可以借助存储表来降低时间复杂度，而借助存储和从头到尾算又是等价的，所以如果采取从头到尾写的话，就是上面链接给出的代码了：

状态转移方程：

dp[i] = max(dp[i - 1], dp[i - 2] + num[i - 1])

其中，dp[i]表示打劫到第i间房屋时累计取得的金钱最大值。

时间复杂度O(n)，空间复杂度O(n)

直白的讲，就是第i步的累积最大值是比较 第i－1步的累积最大值（此时不打劫i） 和 第i-2步累积最大值＋第i步金钱（此时打劫i）。

类似的还有一个好玩的问题：

如果这些房子是相互连城圈的，就是第1间和最后一间连在一起，那么就不能同时打劫第1间和最后一间了，此时的最优策略是什么？

答案无非就是，把一个list分别去掉头保留尾、去掉尾保留头，然后分别算一遍动态规划，看看哪个是最优的就可以了。

今天就写到这里吧，希望日后回头看自己还能看得懂...

想起来差不多十年前从图书馆里面一本又一本的借出来各种算法和数据结构的书籍，却从来没读完过...

今天看到一个东西，格雷码。看了半天硬是没怎么看懂（除了第一种递归的办法）...一看到二进制我就懵掉了，脑子里面一点线索都没有的悲催感。想想这块儿的知识应该也会挺有意思的吧？什么二叉树啊，红黑树啊，好像我都完全不了解是怎么个玩法。唉。

抄一点百度百科的笔记。争取过段时间再来看，能多看懂一些...现在脑子里面完全没有一点感觉。

递归生成码表

这种方法基于格雷码是反射码的事实，利用递归的如下规则来构造：

1. 1位格雷码有两个码字
2. (n+1)位格雷码中的前2ⁿ个码字等于n位格雷码的码字，按顺序书写，加前缀0
3. (n+1)位格雷码中的后2ⁿ个码字等于n位格雷码的码字，按逆序书写，加前缀1^[3] 

2位格雷码	3位格雷码	4位格雷码	4位自然二进制码
00 01 11 10	000 001 011 010 110 111 101 100	0000 0001 0011 0010 0110 0111 0101 0100 1100 1101 1111 1110 1010 1011 1001 1000	0000 0001 0010 0011 0100 0101 0110 0111 1000 1001 1010 1011 1100 1101 1110 1111

异或转换

二进制码→格雷码（编码）：

此方法从对应的n位二进制码字中直接得到n位格雷码码字，步骤如下：

1. 对n位二进制的码字，从右到左，以0到n-1编号
2. 如果二进制码字的第i位和i+1位相同，则对应的格雷码的第i位为0，否则为1（当i+1=n时，二进制码字的第n位被认为是0，即第n-1位不变）^[3] 

公式表示：

（G：格雷码，B：二进制码）

例如：二进制码0101，为4位数，所以其所转为之格雷码也必为4位数，因此可取转成之二进位码第五位为0，即0 b3 b2 b1 b0。 0 xor 0=0，所以g3=0 0 xor 1=1，所以g2=1 1 xor 0=1，所以g1=1 0 xor 1=1，所以g0=1 因此所转换为之格雷码为0111

格雷码→二进制码（解码）：

从左边第二位起，将每位与左边一位解码后的值异或，作为该位解码后的值（最左边一位依然不变）。依次异或，直到最低位。依次异或转换后的值（二进制数）就是格雷码转换后二进制码的值。

公式表示：

（G：格雷码，B：二进制码）

原码：p[n:0]；格雷码：c[n:0](n∈N）；编码：c=G(p）；解码：p=F(c）；

书写时按从左向右标号依次减小，即MSB->LSB，编解码也按此顺序进行

举例： 如果采集器器采到了格雷码：1010 就要将它变为自然二进制： 0 与第四位 1 进行异或结果为 1 上面结果1与第三位0异或结果为 1 上面结果1与第二位1异或结果为 0 上面结果0与第一位0异或结果为 0 因此最终结果为：1100 这就是二进制码即十进制 12 当然人看时只需对照表1一下子就知道是12

...................c[n]=p[n]，

解码：

利用卡诺图

利用卡诺图相邻两格只有一位变化以及卡诺图的变量取值以低阶格雷码的顺序排布的特征，可以递归得到高阶格雷码。由于此方法相对繁琐，使用较少。生成格雷码的步骤如下：

1. 将卡诺图变量分为两组，变量数目相近（最好相等）
2. 以逻辑变量高位在左低位在右建立卡诺图
3. 从卡诺图的左上角以之字形到右上角最后到左下角遍历卡诺图，依次经过格子的变量取值即为典型格雷码的顺序

三位格雷码（三位格雷码由建立在二位基础上）

AB╲ C	0	1
00	0→	1↓
01	↓2	←3
11	6→	7↓
10	4	←5

格雷码次序：000起点→001→011→010→110→111→101→100终点

四位格雷码

AB╲CD	00	01	11	10
00	0→	1→	3→	2↓
01	↓4	←5	←7	←6
11	12→	13→	15→	14↓
10	8	←9	←11	←10

格雷码次序：0000起点→0001→0011→0010→0110→0111→0101→0100→1100→1101→

1111→1110→1010→1011→1001→1000终点

使用异或乘除

用异或代替加减进行二进制竖式乘除，称为异或乘除，它的特点是无进退位。

如：10101除以11将变成1100余1。

二进制转格雷码：

只要异或乘以二分之三，即二进制的1.1，然后忽略小数部分；也可以理解成异或乘以三（即11），再右移一位。

格雷码转二进制：

异或除以三分之二，即除以1.1，忽略余数；或者左移一位，再异或除以三，忽略余数。

虽然最近一直在三藩活动，但是还是有一些一直想去而没机会去的地方。比如今天终于找到时间和理由来一家历史悠久的书店——City Lights Booksellers and Publishers。听名字就是一个很有感觉的地方，城市光线城市之光。以前留意过一些有趣的活动在这里，就大概估计这会是一个类似于诚品那样的书店，可以坐下来慢慢喝茶看书的。结果，我又以东方思维来推测西方生活了。这就是一个书店，一个纯粹的书店，一个不卖任何饮料的书店。

这个书店一进去就是密密麻麻的书架，然后随处可见的标示提醒着浮躁如我的都市人：坐下来，看本书(have a seat + read a book)。于是果断收起手机，随手翻起一本书，寻了一个安静的位置坐下。大概是也没想读一本很厚的诗集，于是随手选了一本薄薄的册子。本也没想读完，却没想到一打开就停不下来，然后竟然就坐在那里翻完了，伴随着心中感慨万千。严格地说，虽然这本书被放在poetry那边，但更像是文化记录的缩影，只是用诗歌这种形式来记录罢了。

这本书的名字是，I am the Beggar of the World, Landays from contemporary Afghanistan. landay大概相当于歌谣，并不一定是用文字（暂不论语言）记录在纸上，而是以声音的形式进行传播和电子记录（这里还要多谢智能手机的发明，很多语音资料得以被很方便的记录下来）。我不是很了解这个书店的定位，粗粗扫了一眼大致以文学书籍为多，还有很多他们自己出版的。直接的感慨就是美国的文化还是很开放的，至少允许人们出版自己的所见所闻。一句题外话，原来只是感慨有墙的存在所以大家翻墙查阅电子资料不方便，上周一个朋友说其实在国内做人文研究也很苦，因为很多国外出版的书籍都买不到、甚至连图书馆里面也很难借到。这样想想，貌似像我这样随便走进书店然后偶遇一本自己可能感兴趣的书，可能也是一件比较奢侈的事情了。不知道这样的隔阂，到底会造福谁。我们甘愿沉浸于禁锢的思想和快餐娱乐吗？