(MIT6.006) lecture11 数值-I

UTF8gbsn

本讲其实主要是在讲数值部分.所以内容偏理论也比较少.

卡塔兰数(Catalan Numbers)

关于这个Catalan numbers的直接定义我们先给出来.
然后再来看看它在实际生活中的一些例子.

$$\begin{array}{rl}{C}_0& =1\\ C_{n+1}& =\sum _{k=0}^n{C}_k{C}_{n-k},n⩾0\end{array}$$

这个数列的枚举如下.

1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 208012, 742900,
2674440, 9694845, 35357670, 129644790, 477638700, 1767263190,
6564120420, 24466267020, 91482563640, 343059613650, 1289904147324,
4861946401452, 18367353072152, 69533550916004, 263747951750360,
1002242216651368

我们来看几个例子.

• 出栈次序问题:一个栈(无穷大)的进栈序列为1,2,3,…n,有多少个不同的出栈序列?

• 找零钱（找一半）:
  有2n个人排成一行进入剧场。入场费5元。其中只有n个人有一张5元钞票，另外n人只有10元钞票，剧院无其它钞票，问有多少中方法使得只要有10元的人买票，售票处就有5元的钞票找零？

• 三角网格 :
  形如这样的直角三角形网格，从左上角开始，只能向右走和向下走，问总共有多少种走法？

• 括号排列: 矩阵连乘：
  $P=a_0\cdot a_1\cdot a_2\cdots a_n$，共有（n+1）项，依据乘法结合律，不改变其顺序，只用括号表示成对的乘积，试问有几种括号化的方案？或者说：有n对括号，可以并列或嵌套排列，共有多少种情况？

• 球盒问题:
  球分两种颜色，黑色和白色分别各有n只，盒子数量和球的个数相同，每个盒子里面只能放一只球，并且必须满足如下限制,即每一个白球必须和一只黑球配对（白球在黑球前，允许嵌套）。

  例.用0表示白猫，1表示黑猫，则：

  0011，010101，001011合法。

  1100，101010，010101不合法。

• 最适合理解的模型:
  n个0和n个1组成一个2n位的2进制数，要求从左到右扫描时，1的累计数始终都小于等于0的累计数，求满足条件的数有多少？

• 类似题目

  • (1).将多边行划分为三角形问题。将一个凸多边形区域分成三角形区域的方法数?

  • (2).有N个节点的二叉树共有多少种情形？

  • (3).一位大城市的律师在她住所以北n个街区和以东n个街区处工作。每天她走2n个街区去上班。如果他从不穿越（但可以碰到）从家到办公室的对角线，那么有多少条可能的道路？

  • (4).在圆上选择2n个点,将这些点成对连接起来使得所得到的n条线段不相交的方法数?

上面的例子来自于百度百科.

牛顿法

牛顿法是一个求根的方法. 我们先来看看迭代公式.

$$x_{i+1}=x_i-\frac{f\left(x_i\right)}{f^{\prime }\left(x_i\right)}$$

乘法的复杂度

朴素的两个n位乘法的复杂度为$\theta (n^2)$, 即便采用下面的方法.

$$\begin{array}{rl}x& =x_1\cdot r^{n/2}+x_0{x}_1=\text{ high half }\\ y& =y_1\cdot r^{n/2}+y_0{x}_0=\text{ low half }\\ 0& \le x_0,x_1<r^{n/2}\\ 0& \le y_0,y_1<r^{n/2}\end{array}$$

$$z=x\cdot y=x_1{y}_1\cdot r^n+\left(x_0\cdot y_1+x_1\cdot y_0\right)r^{n/2}+x_0\cdot y_0$$

时间复杂度也是$\theta (n^2)$

更快的乘法

Karatsuba’s Method:

$$\begin{array}{rl}{z}_0& =x_0\cdot y_0\\ z_2& =x_1\cdot y_1\\ z_1& =\left(x_0+x_1\right)\cdot \left(y_0+y_1\right)-z_0-z_2\\ & =x_0{y}_1+x_1{y}_0\\ z& =z_2\cdot r^n+z_1\cdot r^{n/2}+z_0\end{array}$$

这个时间复杂度为:

$$T(n)=3T(n/2)+\theta (n)=\theta (n^{lo{g}_{2}3})=\theta (n^{1.5849625\cdots })$$

还有更快的算法么? 请自行搜索 Integer multiplication in time $O(nlogn)$