有关数组与函数的矩阵运算

Question

矩阵运算
描述
实现从键盘上读取一个数n，输入n*n的矩阵A，首先判断矩阵A是否为实对称矩阵（实对称矩阵：n阶矩阵A，其矩阵的元素都为实数，且矩阵A的转置等于其本身（aij=aji），(i,j为元素的脚标），则称A为实对称矩阵。）

如果是实对称矩阵则输出 YES，再换行输出矩阵A*A^T的结果（A^T为A的转置） 。

如果不是实对称矩阵则输出NO，再输出矩阵A+A^T的结果（A^T为A的转置）。

以上矩阵输出都以n行n列的格式输出，每个元素用空格隔开，在每行最后一个元素输出换行。

（1）矩阵加法公式：每行每列对应元素相加，Cij=Aij + Bij；

（2）矩阵乘法公式：Cij = Ai0×B0j+Ai1×B1j+……+Ain×Bnj；

输入
从键盘上输入一个数n后再输入n*n的矩阵A。

输出
判断矩阵A是否为实对称矩阵（实对称矩阵：n阶矩阵A，其矩阵的元素都为实数，且矩阵A的转置等于其本身（aij=aji），(i,j为元素的脚标），则称A为实对称矩阵。）

如果是实对称矩阵则输出 YES，再换行输出矩阵A*A^T的结果（A^T为A的转置） ，

如果不是实对称矩阵则输出NO，再换行输出矩阵A+A^T的结果（A^T为A的转置）。

以上矩阵输出都以n行n列的格式输出，每个元素用空格隔开，在每行最后一个元素输出换行。（最后一个元素后面没有空格）

Answer 1

供参考：

#include <stdio.h>
#define N 20
int is_symmetricmatrix(int a[][N], int Row,int Col)//判断一个N*N矩阵，是否是对称矩阵，是返回1，否返回0
{
    int i, j, flag;
    for (i = 0, flag = 1; i < Row; i++)
    {
        for (j = 0; j < i; j++)
        {
            if (a[i][j] != a[j][i])
            {
                flag = 0;
                break;
            }
        }
        if (j < i) break;
    }
    if (flag)
        return 1;
    else
        return 0;
}
void getTransposeMatrix(int b[][N], int a[][N], int Row, int Col)//矩阵转置
{
    int i, j;
    for (i = 0; i < Row; i++)
        for (j = 0; j < Col; j++)
            b[j][i] = a[i][j];
}
void matrixMul(int a[][N], int b[][N], int c[][N], int Rowa, int Cola, int Colb)//矩阵相乘
{
    int i, j, k;
    for(i=0;i<Rowa;i++)
    {
        for (j = 0; j < Colb; j++)
        {
            for (k = 0; k < Cola; k++)
                c[i][j] += a[i][k] * b[k][j];
        }
    }
}
void matrixAdd(int a[][N], int b[][N], int c[][N], int Rowa, int Cola)//矩阵相加
{
    int i, j;
    for (i = 0; i < Rowa; i++)
        for (j = 0; j < Cola; j++)
            c[i][j] = a[i][j] + b[i][j];
}
void printmatrix(int a[][N], int Row, int Col)
{
    int i, j;
    for (i = 0; i < Row; i++){
        for (j = 0; j < Col; j++)
        {
            printf("%d ", a[i][j]);
        }
        printf("\n");
    }
}
int main()
{
    int a[N][N] = { 0 }, b[N][N] = { 0 }, c[N][N] = {0}, i, j, n = 3;
    scanf("%d", &n);
    for (i = 0; i < n; i++)
        for (j = 0; j < n; j++)
            scanf("%d", &a[i][j]);
    if (is_symmetricmatrix(a, n, n))
    {
        printf("YES\n");
        getTransposeMatrix(b, a, n, n);
        matrixMul(a, b, c, n, n, n);
        printmatrix(c, n, n);
    }
    else{
        printf("NO\n");
        getTransposeMatrix(b, a, n, n);
        matrixAdd(a, b, c, n, n);
        printmatrix(c, n, n);
    }
    return 0;
}