西南交通大学数据结构第五次实验报告-

实验目的：掌握半三角矩阵的顺序存储结构。

实验代码：

#include "iostream"

using namespace std;

int posi(int i, int j) {
    int length = 0;
    for (int count = 0; count < i; count++) {
        length += count;
    }
    length += j - 1;
    return length;
}

void matrix_multiplication(const int *matrix1, const int* matrix2, int* matrix_sum, int n) {
    for (int i = 0; i <= n; i++) {
        for (int j = 0; j < n; j++) {
            for (int k = j; k <= i; k++) {
                matrix_sum[posi(i + 1, j + 1)] += matrix1[posi(k + 1, j + 1)] * matrix2[posi(i + 1, k + 1)];
            }
        }
    }
}

void readMatrix(int matrix[], int n) {
    for (int i = 0; i < n; i++)
        for (int j = 0; j <= i; j++) {
            cin >> matrix[i * (i + 1) / 2 + j];
        }
}

int main() {
    int n;
    cin >> n;
    int *matrix1, *matrix2, *matrix_sum;
    int length = 0;
    for (int i = 1; i <= n; i++) {
        length += i;
    }
    matrix1 = new int[length];
    matrix2 = new int[length];
    matrix_sum = new int[length];
    readMatrix(matrix1, n);
    readMatrix(matrix2, n);
    for (int i = 0; i < length; i++) {
        matrix_sum[i] = 0;
    }
    matrix_multiplication(matrix1, matrix2, matrix_sum, n);
    for (int i = 0; i < n; i++) {
        for (int j = 0; j <= i; j++) {
            cout << matrix_sum[posi(i + 1, j + 1)] << " ";
        }
        cout << endl;
    }
    return 0;
}

程序有一点小瑕疵，在输入用例过后你就知道是什么了，结果是正确的，但是不知道怎么解决掉这个小瑕疵，有知道的佬告诉一声，谢谢。