第八周实践项目7 对称矩阵的压缩存储及基本运算
生活随笔
收集整理的這篇文章主要介紹了
第八周实践项目7 对称矩阵的压缩存储及基本运算
小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
/*
*Copyright (c) 2017,煙臺(tái)大學(xué)計(jì)算機(jī)與控制工程學(xué)院
*All rights reserved.
*文件名稱:項(xiàng)目7-
*作 者:邵雪源
*完成日期:2017年12月14日
*版 本 號(hào):v1.0
*/
#include <stdio.h>
#include <malloc.h>
#define N 4void Init(int *&b);//為N階對(duì)稱矩陣初始化存儲(chǔ)數(shù)據(jù)的一維數(shù)組b
int Value(int b[], int i, int j);//返回存儲(chǔ)在b[M]中,對(duì)應(yīng)二維數(shù)組A[i][j]的值
void Assign(int b[], int e, int i, int j);//將e賦值給對(duì)應(yīng)二維數(shù)組元素A[i][j],要存儲(chǔ)到b[M]中
void Disp(int b[]);//輸出壓縮存儲(chǔ)在b中的對(duì)稱矩陣
void Destroy(int b[]); //銷毀存儲(chǔ)空間
//為N階對(duì)稱矩陣初始化存儲(chǔ)數(shù)據(jù)的一維數(shù)組B
void Init(int *&b)
{b = (int*)malloc(sizeof(int)*(N*(N+1)/2));
}
//返回存儲(chǔ)在b[M]中,對(duì)應(yīng)二維數(shù)組A[i][j]的值
int Value(int b[], int i, int j)
{if (i>=j)return b[(i*(i+1))/2+j];elsereturn b[(j*(j+1))/2+i];
}
//將e賦值給對(duì)應(yīng)二維數(shù)組元素A[i][j],要存儲(chǔ)到B[M]中
void Assign(int b[], int e, int i, int j)
{if (i>=j)b[(i*(i+1))/2+j] = e;elseb[(j*(j+1))/2+i] = e;return;
}
//輸出壓縮存儲(chǔ)在b中的對(duì)稱矩陣
void Disp(int b[])
{int i,j;for (i=0; i<N; i++){for (j=0; j<N; j++)printf("%4d",Value(b,i,j));printf("\n");}
}
//銷毀存儲(chǔ)空間
void Destroy(int b[])
{free(b);
}
int main()
{int *b1; //指向整型的指針,待初始化int i, j;int v;Init(b1);printf("請(qǐng)輸入對(duì)稱矩陣(只需要輸入下三角部分即可)\n");for(i=0;i<N;i++){printf("輸入第%d行的%d個(gè)數(shù)據(jù)元素: ", i+1, i+1);for(j=0; j<=i; j++){scanf("%d", &v);Assign(b1, v, i, j);}}Disp(b1);Destroy(b1);return 0;
}
(2)設(shè)計(jì)算法,實(shí)現(xiàn)兩個(gè)用壓縮形式存儲(chǔ)的對(duì)稱矩陣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];elsereturn 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表示壓縮存儲(chǔ)的對(duì)稱矩陣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; }
#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];elsereturn 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表示壓縮存儲(chǔ)的對(duì)稱矩陣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; }
總結(jié)
以上是生活随笔為你收集整理的第八周实践项目7 对称矩阵的压缩存储及基本运算的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 第八周实践项目6 猴子选大王(数组版)
- 下一篇: 第八周实践项目8 稀疏矩阵的三元组表示的