第四章 4.2矩阵乘法的Strassen算法（分治）

注意：虽然联系4.2-3、4.2-4和4.2-5是关于Strassen算法的变形的，但你应该先阅读4.5节，然后再尝试求解这几个问题。
####4.2-1 使用Strassen算法计算如下矩阵乘法：$$\left[\begin{array}{cc}1& 3\\ 7& 5\end{array}\right]\left[\begin{array}{cc}6& 8\\ 4& 2\end{array}\right]$$给出计算过程。

$$\begin{gathered} S_1=B_{12}-B_{22}=8-2=6 \\ {S}_2=A_{11}+A_{12}=1+3=4 \\ {S}_3=A_{21}+A_{22}=7+5=12 \\ {S}_4=B_{21}-B_{11}=4-6=-2 \\ {S}_5=A_{11}+A_{22}=1+5=6 \\ {S}_6=B_{11}+B_{22}=6+2=8 \\ {S}_7=A_{12}-A_{22}=3-5=-2 \\ {S}_8=B_{21}+B_{22}=4+2=6 \\ {S}_9=A_{11}-A_{21}=1-7=-6 \\ {S}_{10}=B_{11}+B_{12}=6+8=14 \\ {P}_1=A_{11}\cdot S_1=6 \\ {P}_2=S_2\cdot B_{22}=8 \\ {P}_3=S_3\cdot B_{11}=72 \\ {P}_4=A_{22}\cdot S_4=-10 \\ {P}_5=S_5\cdot S_6=48 \\ {P}_6=S_7\cdot S_8=-12 \\ {P}_7=S_9\cdot S_{1}0=-84 \\ {C}_{11}=P_5+P_4-P_2+P_6=48-10-8-12=18 \\ {C}_{12}=P_1+P_2=6+8=14 \\ {C}_{21}=P_3+P_4=72-10=62 \\ {C}_{22}=P_5+P_1-P_3-P_7=48+6-72+84=66 \\ 综上所述，计算出的矩阵值为\left[\begin{array}{cc}18& 14\\ 62& 66\end{array}\right] \end{gathered}$$

####4.2-2 为Strassen算法编写伪代码。

Strassen(A,B)
n=A.rows
let C be a new nxn matrix
if n==1
	c=a*b
else partition A,B,and C as in equations(4,9)
	S1=B12-B22
	S2=A11-A12
	S3=A21+A22
	S4=B21-B11
	S5=A11+A22
	S6=B11+B22
	S7=A12-A22
	S8=B21+B22
	S9=A11-A21
	S10=B11+B12
	P1=Strassen(A11,S1)
	P2=Strassen(A11,B22)
	P3=Strassen(S3,B11)
	P4=Strassen(A22,S4)
	P5=Strassen(S5,S6)
	P6=Strassen(S7,S8)
	P7=Strassen(S9,S10)
	C11=P5+P4-P2+P6
	C12=P1+P2
	C21=P3+P4
	C22=P5+P1-P3-P7
	return C

####4.2-3 如何修改Strassen算法，使之适应矩阵规模n不是2的幂的情况？证明：算法的运行时间为$\theta (n^{lg7})$。
将矩阵先补成2的幂情况。
证明：http://www.mamicode.com/info-detail-673908.html
####4.2-4 如果可以用k次乘法操作（假定乘法的交换律不成立）完成两个3x3矩阵相乘，那么你可以在O（$n^{lg7}$）时间内完成nxn的矩阵相乘，满足这一条件的最大的k是多少？此算法的运行时间是怎样的？
将矩阵补成44的矩阵
7+47=35=k
O（$4^{lg7}$）
####4.2-5 V.Pan发现一种方法，可以用132464次乘法操作完成68x68的矩阵相乘，发现另一种方法，可以用143 640次乘法操作完成70x70的矩阵相乘，还发现一种方法，可以用155 424 次乘法操作完成72x72的矩阵相乘。当用于矩阵相乘的分治算法时，上述哪种方法会得到最佳的渐进运行时间？与Strassen算法相比，性能如何？
Strassen运行次数为：135215，146647，158682（就比例来说第三种方法好一些，但算法的的控制因素有两个，不能直接因此做判断）
目前次数的运算来说，好于Strassen，但算法好坏好要靠其他因素综合来判断。
####4.2-6 用Strassen算法作为子进程来进行一个kn X n矩阵和一个nXkn的矩阵相乘，最快需要花费多长时间？对两个输入矩阵规模互换的情况，回答相同的问题。
1）先补成kn*kn，再补成齐二次幂。
2）相同
####4.2-7 设计算法，仅使用三次实数乘法即可完成复数a+bi和c+di相乘。算法需接受a、b、c和d为输入，分别生成实部ac-bd和虚部ad+bc。
(a+b)(c-d)、ad、bc