第九周專案2—矩陣運算
阿新 • • 發佈:2019-01-23
問題及程式碼
*Copyright(c)2016,煙臺大學計算機學院 *檔名稱:矩陣運算.cpp *作者:郗傳秀 *完成日期;2016年10月27日 *問題描述: 設計演算法,實現兩個用壓縮形式儲存的對稱矩陣A和B的加法和乘法。 實現中請使用好前面設計的基本運算。 *輸入描述:兩個矩陣 *程式輸出:兩個矩陣相加與相乘的結果 */ #include <stdio.h> #define N 4 #define M 10 int value(int a[],int i,int j) { if (i>=j) return a[(i*(i+1))/2+j]; else return a[(j*(j+1))/2+i]; } void madd(int a[],int b[],int c[][N]) { int i,j; for (i=0; i<N; i++) for (j=0; j<N; j++) c[i][j]=value(a,i,j)+value(b,i,j); } void mult(int a[],int b[],int c[][N]) { int i,j,k,s; for (i=0; i<N; i++) for (j=0; j<N; j++) { s=0; for (k=0; k<N; k++) s=s+value(a,i,k)*value(b,k,j); c[i][j]=s; } } void disp1(int a[]) { int i,j; for (i=0; i<N; i++) { for (j=0; j<N; j++) printf("%4d",value(a,i,j)); printf("\n"); } } void disp2(int c[][N]) { int i,j; for (i=0; i<N; i++) { for (j=0; j<N; j++) printf("%4d",c[i][j]); printf("\n"); } } int main() { int a[M]= {1,2,3,4,5,6,7,8,9,10}; //a表示壓縮儲存的對稱矩陣 int b[M]= {1,1,1,1,1,1,1,1,1,1}; int c1[N][N],c2[N][N]; madd(a,b,c1); mult(a,b,c2); printf("a矩陣:\n"); disp1(a); printf("b矩陣:\n"); disp1(b); printf("a+b:\n"); disp2(c1); printf("a×b:\n"); disp2(c2); printf("\n"); return 0; }
執行結果