基于压缩存储的半三角矩阵乘法运算的实现

这是一个简便版,供参考

已知两个n阶下半三角矩阵的乘积仍为n阶下半三角矩阵。编程输入两个n阶下半三角矩阵，输出这两个矩阵的乘积。要求n阶下半三角矩阵采用一维数组压缩存储（即只存储下半三角）。

程序先从键盘(或字符文件)输入n值，建立三个矩阵的一维数组动态存储结构，然后从键盘(或字符文件)输入两个半三角矩阵，最后输出计算结果到屏幕上(或另一个字符文件中)。

例如：键盘输入为：

3

1

2 3

4 5 6

-1

-2 -3

-4 -5 -6

则输出为：

-1

-8   -9

-38 -45 -36

#define _CRT_SECURE_NO_WARNINGS 1
#include <iostream>
#include <fstream>
#include <iomanip>
#include<string>
#define filename1 "x1.txt"
#define filename2 "x2.txt"
#define filename3 "y.txt"
using namespace std;

//读取文件,a为矩阵名,filename为文件名,n为矩阵阶数
void Read(int* a, string filename, int n)
{
    int i = 0, c;
    ifstream infile(filename);
    if (!infile)                        //打开文件失败
    {
        cout << "File open error!" << endl;
        exit(0);
    }
    while (1)
    {
        infile >> c;
        a[i] = c;
        if (infile.eof())
            break;
       
        if (i + 1 == n * (n + 1) / 2)  
            break;
        i++;
    }
    infile.close();
}
//写入文件,a为矩阵名,filename为文件名,n为矩阵阶数
void Write(int* a, string filename, int n)
{
    int i = 0, j = 0;
    ofstream outfile(filename);
    if (!outfile)                     //打开文件失败
    {
        cout << "File open error!" << endl;
        exit(0);
    }
    for (i = 0, j = 1; i < n * (n + 1) / 2; i++)    //将元素按矩阵形式存到文件中
    {
        outfile << setw(5) << a[i];
        if ((i + 1) == j * (j + 1) / 2)
        {
            j++;
            outfile << endl;           //每行元素达到最大，输出换行
        }
    }
    outfile.close();
}
//两个矩阵相乘,矩阵c=矩阵a*矩阵b,n为矩阵阶数
void Mutiply(int* a, int* b, int* c, int n)
{
    //m记录数组c的下标，sum为矩阵c每个位置元素的计算结果
    int i, j, k, m = 0, sum = 0;    
    //从 0 到 n 按行扫描
    for (i = 0; i < n; i++)        
    {
        //从 0 到 i 按列扫描
        for (j = 0; j <= i; j++)    
        {
            //从 j 到 i 将矩阵 a 与矩阵 b 相乘
            for (k = j; k <= i; k++)
            {
                sum += a[i * (i + 1) / 2 + k] * b[k * (k + 1) / 2 + j];
                c[m] = sum;
            }
            sum = 0;
            m++;
        }
    }
}
//输出函数，a为矩阵名，n为矩阵阶数
void Show(int* a, int n)          
{
    int i, j;
    for (i = 0, j = 1; i <n * (n + 1) / 2; i++)
    {
        cout << setw(5) << a[i];
        if ((i + 1) == j * (j + 1) / 2)      //将元素按矩阵形式输出
        {
            j++;
            cout << endl;
        }
    }
}
void test01() {
    int* a1, * a2, * a3, n, m = 0, j = 0;
    cout << "请输入n的值: ";             
    //输入矩阵阶数
    cin >> n;
    //分别建立三个数组，n * (n + 1) / 2为数组中的元素个数
    a1 = new int[n * (n + 1) / 2];         
    a2 = new int[n * (n + 1) / 2];
    a3 = new int[n * (n + 1) / 2];
    //初始化所有元素为0
    a1[n * (n + 1) / 2] = a2[n * (n + 1) / 2] = a3[n * (n + 1) / 2] = { 0 };
    Read(a1, filename1, n);                   
    Read(a2, filename2, n);               
    cout << "矩阵A为: " << endl;
    Show(a1, n);                           
    cout << "矩阵B为: " << endl;
    Show(a2, n);                           
    Mutiply(a1, a2, a3, n);                   
    cout << "矩阵C为: " << endl;
    Show(a3, n);                           
    Write(a3, filename3, n);                
}
int main()
{
    test01();
    return 0;
}

需要在文件x1,x2中写好矩阵!!!