Data Structures实习四 稀疏矩阵.pptVIP

Data Structures實習四 稀疏矩陣 Instructor: Ching-Chi Lin林清池 助理教授 chingchi.lin@ Department of Computer Science and EngineeringNational Taiwan Ocean University Outline * 稀疏矩陣介紹 稀疏矩陣的儲存 稀疏矩陣的乘法 稀疏矩陣介紹 * 稀疏矩陣： 很多元素都是零。 只存矩陣的大小和非零的元素以節省空間。 請寫一個程式，儲存並處理稀疏矩陣的乘法。 將稀疏矩陣中非零的元素存在陣列中。 每個元素需要儲存的資料為：所在的行，列，以及值。 矩陣的大小與元素量也需要被儲存。 稀疏矩陣的儲存 * #define MAX_TERMS 101 每個矩陣最多可以有的非零變數+1 typedef struct { 　　unsigned int row; 　　unsigned int col; 　　int value; 　}term; term matrix a[MAX_TERMS]; 稀疏矩陣的儲存 * row col value a[0] 3 4 5 a[1] 0 0 3 a[2] 0 2 1 a[3] 1 0 1 a[4] 1 1 3 a[5] 2 0 -2 a[6] 0 0 0 … … … … 為一3(列)*4(行)的矩陣，並且有5個非零元素。 0 1 2 3 0 3 0 1 0 1 1 3 0 0 2 -2 0 0 0 矩陣的乘法 * 0 1 2 0 3 0 1 1 1 0 0 2 2 0 3 3 0 0 4 0 1 2 3 0 3 0 1 0 1 1 3 0 0 2 -2 0 0 0 * 0 1 2 0 11 0 6 1 6 0 1 2 -6 0 -1 = 稀疏矩陣的乘法 * 在乘法（a*b）時： 需要找出b矩陣中同column的元素。 儲存資料的排序方法為row - column。 為了加快速度，先對b陣列做轉置（transpose）。 對矩陣中的所有元素，Ai.j=ATj,i。 轉置過後的矩陣資料，依然依照row - column的順序排序。 稀疏矩陣的轉置 * 0 1 2 3 0 0 3 0 1 1 0 1 3 4 2 0 0 -2 0 row col value a[0] 3 4 6 a[1] 0 1 3 a[2] 0 3 1 a[3] 1 1 1 a[4] 1 2 3 a[5] 1 3 4 a[6] 2 2 -2 0 1 2 3 0 3 -1 5 0 1 1 0 0 2 0 1 0 3 1 1 -1 0 2 5 0 3 0 2 row col value b[0] 4 2 5 b[1] 0 0 3 b[2] 0 1 1 b[3] 1 0 -1 b[4] 2 0 5 b[5] 3 1 2 row col value b_T[0] 2 4 5 b_T[1] 0 0 3 b_T{2] 0 1 -1 b_T[3] 0 2 5 b_T[4] 1 0 1 b_T[5] 1 3 2 轉置 稀疏矩陣轉置後相乘 * 0 1 2 3 0 0 3 0 1 1 0 1 3 4 2 0 0 -2 0 row col value a[0] 3 4 6 a[1] 0 1 3 a[2] 0 3 1 a[3] 1 1 1 a[4] 1 2 3 a[5] 1 3 4 a[6] 2 2 -2 0 1 2 3 0 3 -1 5 0 1 1 0 0 2 row col value b_T[0] 2 4 5 b_T[1] 0 0 3 b_T{2] 0 1 -1 b_T[3] 0 2 5 b_T[4] 1 0 1 b_T[5] 1 3 2 row col Sum= 0 -3 (0,0)=-3 col 2 (0,1)=2 row 稀疏矩陣乘法的PseudoCode * 請參考課程投影片Ch2： The Sparse Matrix Abstract Data Type之部分。 稀疏矩陣乘法的練習 * 一般矩陣相乘：A為m*n矩陣、B為n*p矩陣，D（m*p）＝A*B，其元素i,j為： 　 　 稀疏矩陣乘法的練習 * 輸入：matrix_a.txt / matrix_b.txt 為兩個矩陣A/B。 請至教學網站下載範例檔。 每一元素之間以換行符號(\n)相隔。 第一行記錄矩陣的行列數以及非0元素個數。 每一元素的格式為： 行(非負整數)\t列(非負整數)\t值(整數)。 檔案為排序的。 總長度不超過100個元素。 值的大小介於 0 ~ 16,777,215 (FFFFFF)。 輸出：分別印出A與B矩陣， B的轉置矩陣，