西电数据结构上机题——对称矩阵相乘
【摘要】
这道题拿到是懵逼的 本题最为关键的是对称矩阵相乘的算法 幸好有老哥之前探索出了 对称矩阵M的第i行和第j列的元素的数据存储在一维数组a中的位置k的计算公式: 1、当i大于或等于j时,k = (i * (i...
这道题拿到是懵逼的
本题最为关键的是对称矩阵相乘的算法
幸好有老哥之前探索出了
对称矩阵M的第i行和第j列的元素的数据存储在一维数组a中的位置k的计算公式:
1、当i大于或等于j时,k = (i * (i + 1)) / 2 + j (下三角)
2、当i小于j时,k = (j * (j + 1)) / 2 + i (上三角)
(注意这里是整除,真的是非常Amazing,有时间可以去研究一下是怎么推出来的)
链接: https://blog.csdn.net/xiezhi123456/article/details/86607261
在他的基础上顺利解决
//对称矩阵相乘的程序代码
#include<stdio.h>
#include<malloc.h>
//数组结构类型的定义.h
const int n=3;
const int size=n*(n+1)/2;
typedef int datatype;
typedef struct{
datatype A[size],B[size],C[n][n];
}array;
void input(datatype[]);
void output(datatype[][n]);
void mult(array*);
int main()
{
array*pa;
pa=(array*)malloc(sizeof(array));
printf("以行为主序输入矩阵A的下三角:\n");
input(pa->A);//以行为主序输入矩阵A的下三角
printf("以行为主序输入矩阵B的下三角:\n");
input(pa->B);//以行为主序输入矩阵B的下三角
mult(pa);
output(pa->C);//输出矩阵C
}
//对称矩阵的输入
void input(datatype x[])
{
for(int i=0;i<size;i++)
scanf("%d",&x[i]);
}
//矩阵的输出
void output(datatype x[][n])
{
for(int i=0;i<n;i++){
for(int j=0;j<n;j++)
printf("%5d",x[i][j]);
printf("\n");
}
}
datatype 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 mult(array*pa)
{
int i,j,k,sum;
for (i=0;i<n;i++)
{
for(j=0;j<n;j++)
{
sum=0;
for(k=0;k<n;k++)
{
sum=sum+value(pa->A,i,k)*value(pa->B,k,j);
}
pa->C[i][j]=sum;
}
}
}
- 1
- 2
- 3
- 4
- 5
- 6
- 7
- 8
- 9
- 10
- 11
- 12
- 13
- 14
- 15
- 16
- 17
- 18
- 19
- 20
- 21
- 22
- 23
- 24
- 25
- 26
- 27
- 28
- 29
- 30
- 31
- 32
- 33
- 34
- 35
- 36
- 37
- 38
- 39
- 40
- 41
- 42
- 43
- 44
- 45
- 46
- 47
- 48
- 49
- 50
- 51
- 52
- 53
- 54
- 55
- 56
- 57
- 58
- 59
- 60
- 61
- 62
- 63
- 64
- 65
- 66
- 67
- 68
- 69
- 70
- 71
文章来源: zstar.blog.csdn.net,作者:zstar-_,版权归原作者所有,如需转载,请联系作者。
原文链接:zstar.blog.csdn.net/article/details/109631763
【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱:
cloudbbs@huaweicloud.com
- 点赞
- 收藏
- 关注作者
评论(0)