Gambar yang lebih besar di belakang pilihan matriks dalam algoritma Strassen

Dalam algoritma Strassen, untuk menghitung produk dari dua matriks dan , matriks dan dibagi menjadi matriks blok dan algoritma menghasilkan perhitungan secara rekursif blok produk matriks-matriks sebagai lawan dari blok matriks- naif produk matriks, yaitu, jika kita ingin , di mana $\mathbf{A}$ $\mathbf{B}$ $\mathbf{A}$ $\mathbf{B}$ $2 \times 2$ $7$ $8$ $\mathbf{C}=\mathbf{A} \mathbf{B}$ maka kita memiliki

A = [\begin{matrix} A_{1, 1} & A_{1, 2} \\ A_{2, 1} & A_{2, 2} \end{matrix}], B = [\begin{matrix} B_{1, 1} & B_{1, 2} \\ B_{2, 1} & B_{2, 2} \end{matrix}], C = [\begin{matrix} C_{1, 1} & C_{1, 2} \\ C_{2, 1} & C_{2, 2} \end{matrix}]

$\mathbf{A} =\begin{bmatrix} \mathbf{A}_{1,1} & \mathbf{A}_{1,2} \\ \mathbf{A}_{2,1} & \mathbf{A}_{2,2} \end{bmatrix} \mbox { , } \mathbf{B} = \begin{bmatrix} \mathbf{B}_{1,1} & \mathbf{B}_{1,2} \\ \mathbf{B}_{2,1} & \mathbf{B}_{2,2} \end{bmatrix} \mbox { , } \mathbf{C} = \begin{bmatrix} \mathbf{C}_{1,1} & \mathbf{C}_{1,2} \\ \mathbf{C}_{2,1} & \mathbf{C}_{2,2} \end{bmatrix}$

yang membutuhkan

perkalian. Sebagai gantinya di Strassen, kami menghitung

C_{1, 1} = {SEBUAH}_{1, 1} B_{1, 1} + {SEBUAH}_{1, 2} B_{2, 1} C_{1, 2} = {SEBUAH}_{1, 1} B_{1, 2} + {SEBUAH}_{1, 2} B_{2, 2} C_{2, 1} = {SEBUAH}_{2, 1} B_{1, 1} + {SEBUAH}_{2, 2} B_{2, 1} C_{2, 2} = {SEBUAH}_{2, 1} B_{1, 2} + {SEBUAH}_{2, 2} B_{2, 2}

$\mathbf{C}_{1,1} = \mathbf{A}_{1,1} \mathbf{B}_{1,1} + \mathbf{A}_{1,2} \mathbf{B}_{2,1}\\ \mathbf{C}_{1,2} = \mathbf{A}_{1,1} \mathbf{B}_{1,2} + \mathbf{A}_{1,2} \mathbf{B}_{2,2}\\ \mathbf{C}_{2,1} = \mathbf{A}_{2,1} \mathbf{B}_{1,1} + \mathbf{A}_{2,2} \mathbf{B}_{2,1}\\ \mathbf{C}_{2,2} = \mathbf{A}_{2,1} \mathbf{B}_{1,2} + \mathbf{A}_{2,2} \mathbf{B}_{2,2}$

8

$8$

dan dapatkan

menggunakan

sebagai

{M.}_{1} : = ({SEBUAH}_{1, 1} + {SEBUAH}_{2, 2}) (B_{1, 1} + B_{2, 2}) {M.}_{2} : = ({SEBUAH}_{2, 1} + {SEBUAH}_{2, 2}) B_{1, 1} {M.}_{3} : = {SEBUAH}_{1, 1} (B_{1, 2} - B_{2, 2}) {M.}_{4} : = {SEBUAH}_{2, 2} (B_{2, 1} - B_{1, 1}) {M.}_{5} : = ({SEBUAH}_{1, 1} + {SEBUAH}_{1, 2}) B_{2, 2} {M.}_{6} : = ({SEBUAH}_{2, 1} - {SEBUAH}_{1, 1}) (B_{1, 1} + B_{1, 2}) {M.}_{7} : = ({SEBUAH}_{1, 2} - {SEBUAH}_{2, 2}) (B_{2, 1} + B_{2, 2})

$\mathbf{M}_{1} := (\mathbf{A}_{1,1} + \mathbf{A}_{2,2}) (\mathbf{B}_{1,1} + \mathbf{B}_{2,2})\\ \mathbf{M}_{2} := (\mathbf{A}_{2,1} + \mathbf{A}_{2,2}) \mathbf{B}_{1,1}\\ \mathbf{M}_{3} := \mathbf{A}_{1,1} (\mathbf{B}_{1,2} - \mathbf{B}_{2,2})\\ \mathbf{M}_{4} := \mathbf{A}_{2,2} (\mathbf{B}_{2,1} - \mathbf{B}_{1,1})\\ \mathbf{M}_{5} := (\mathbf{A}_{1,1} + \mathbf{A}_{1,2}) \mathbf{B}_{2,2}\\ \mathbf{M}_{6} := (\mathbf{A}_{2,1} - \mathbf{A}_{1,1}) (\mathbf{B}_{1,1} + \mathbf{B}_{1,2})\\ \mathbf{M}_{7} := (\mathbf{A}_{1,2} - \mathbf{A}_{2,2}) (\mathbf{B}_{2,1} + \mathbf{B}_{2,2})$

C_{i, j}

$\mathbf{C}_{i,j}$

M_{k}

$\mathbf{M}_{k}$

Namun, pilihan matriks

tampak sewenang-wenang bagi saya. Apakah ada gambaran yang lebih besar mengapa kami memilih produk-produk khusus dari sub-matriks

dan

? Juga, saya harapkan

's untuk melibatkan

' s dan

's secara simetris, yang tampaknya tidak menjadi kasus di sini. Misalnya, kami memiliki

C_{1, 1} = {M.}_{1} + {M.}_{4} - {M.}_{5} + {M.}_{7} C_{1, 2} = {M.}_{3} + {M.}_{5} C_{2, 1} = {M.}_{2} + {M.}_{4} C_{2, 2} = {M.}_{1} - {M.}_{2} + {M.}_{3} + {M.}_{6}

$\mathbf{C}_{1,1} = \mathbf{M}_{1} + \mathbf{M}_{4} - \mathbf{M}_{5} + \mathbf{M}_{7}\\ \mathbf{C}_{1,2} = \mathbf{M}_{3} + \mathbf{M}_{5}\\ \mathbf{C}_{2,1} = \mathbf{M}_{2} + \mathbf{M}_{4}\\ \mathbf{C}_{2,2} = \mathbf{M}_{1} - \mathbf{M}_{2} + \mathbf{M}_{3} + \mathbf{M}_{6}$

M_{k}

$\mathbf{M}_k$

A

$\mathbf{A}$

B

$\mathbf{B}$

M_{k}

$\mathbf{M}_k$

A_{i, j}

$\mathbf{A}_{i,j}$

B_{i, j}

$\mathbf{B}_{i,j}$

M_{2} := (A_{2, 1} + A_{2, 2}) B_{1, 1}

$\mathbf{M}_2: = (\mathbf{A}_{2,1}+\mathbf{A}_{2,2})\mathbf{B}_{1,1}$

A_{1, 1} (B_{1, 2} + B_{2, 2})

$\mathbf{A}_{1,1} (\mathbf{B}_{1,2} + \mathbf{B}_{2,2})$

M_{k}

$\mathbf{M}_k$

Saya akan sangat menghargai jika seseorang dapat menjelaskan hal ini.

ds.algorithms linear-algebra matrices matrix-product Komunitas
sumber

Jawaban:

De Groote (Pada Varietas Algoritma Optimal untuk Penghitungan Pemetaan Bilinear. II. Algoritma Optimal untuk Penggandaan Matriks 2x2. Teori. Komputasi. Sci. 7: 127-148, 1978) membuktikan bahwa hanya ada satu algoritma yang dapat dikalikan $2 \times 2$ -matrices dengan 7 perkalian hingga kesetaraan. Ini mungkin fitur unik dari $2 \times 2$ perkalian -matrix. (Catatan: Anda akan menemukan varian berbeda dari algoritma Strassen dalam literatur. Semuanya setara dengan gagasan persamaan yang tepat.)

Jika Anda sekarang mulai membuktikan batas bawah untuk $2 \times 2$ perkalian -matrix - lihat buku karya Bürgisser, Clausen, dan Shokrollahi bagaimana melakukan itu - kemudian algoritma Strassen atau beberapa varian muncul secara alami. Anda akan menemukan banyak identitas yang menentukan bagaimana produk terlihat. Kemudian Anda bisa menyelesaikannya dengan menebak-nebak. (Bukti De Groote menunjukkan bahwa menebak sekalipun tidak perlu.)

Schönhage pernah mengatakan kepada saya bahwa Strassen pernah mengatakan kepadanya bahwa ia menemukan algoritme dengan cara ini, dengan mencoba membuktikan batas bawah.

Markus Bläser
sumber

Ada semacam penjelasan dalam buku Teori Kompleksitas Aljabar oleh Bürgisser, Clausen dan Shokrollahi (hlm. 11-12). Idenya adalah memulai dengan dua basis $A_0,A_1,A_2,A_3$ dan $B_0,B_1,B_2,B_3$ dari ruang $2\times 2$ matriks nyata yang memenuhi properti berikut: $A_iB_j \in \{0,A_0,A_1,A_2,A_3,B_0,B_1,B_2,B_3\}$ . Selanjutnya, $A_0 = B_0$ . Untuk melipatgandakan dua matriks $A$ dan $B$ , mewakili masing-masing dari mereka dalam dasar yang sesuai, dan mengevaluasi produk. Karena hanya tujuh matriks non-nol yang berbeda muncul dalam hasil ( $A_0=B_0,A_1,A_2,A_3,B_1,B_2,B_3$ ), hanya tujuh produk yang dibutuhkan. Itu $M$ matriks hanyalah pangkalan-pangkalan ini.

Saya tidak tahu apakah Strassen datang dengan cara memandangnya seperti ini. Mengingat identitas lain yang mendasari algoritma perkalian matriks cepat, tidak jelas apakah ada sesuatu yang mendalam, di bawah daripada beberapa rumus yang bekerja. Kami telah mengalaminya sebelumnya - Lagrange menggunakan identitas empat persegi (yang telah dikenal sebelumnya) untuk membuktikan teorema empat persegi. Pada awalnya itu pasti hanya identitas aljabar yang aneh, tetapi sekarang kita tahu bahwa itu menyatakan properti multiplikatif dari norma angka empat. Mengingat keadaan pengetahuan saat ini, sulit untuk mengatakan apakah interpretasi di atas sama produktifnya.

Yuval Filmus
sumber

Basis semacam itu disebut pasangan-M, lihat bab tentang aljabar peringkat minimal dalam buku karya Bürgisser, Clausen, dan Shokrollahi. Saya pikir cukup sulit untuk mengemukakan gagasan bahwa pasangan-M ada tanpa mengetahui teorema Alder-Strassen (lihat lagi buku di atas). Khususnya,

2 \times 2

$2\times 2$ -matrices adalah satu-satunya aljabar matriks yang memiliki pasangan-M.

Markus Bläser