矩阵乘法 C++实现

#include <iostream>
#include <iomanip>
using namespace std;
int main()
{
    int M=0;
    int a;
    int i,j,k;
    cout<<"请输入矩阵维数"<<endl;
    cin>>a;

    int **A;
    A=new int*[a];  //指针数组
    for (i=0;i<a;i++)
    {
        A[i]=new int[a];
    }

    int **B;
    B=new int*[a];
    for (i=0;i<a;i++)
    {
        B[i]=new int[a];
    }

    int **C;
    C=new int*[a];
    for (i=0;i<a;i++)
    {
        C[i]=new int[a];
    }
    
    cout<<"请输入矩阵A"<<endl;
    for(i=0;i<a;i++)   //给A矩阵赋值
    {
        for(j=0;j<a;j++)
        {
            cin>>A[i][j];
        }
    }
    cout<<"您输入的第一个矩阵A为："<<endl;
    for(i=0;i<a;i++)                                               
    {
        for(j=0;j<a;j++)
        {
            cout<<'\t';
            if((j+1)%a==0)
                cout<<A[i][j]<<endl;
            else
                cout<<A[i][j];
        }
    }
    cout<<"请输入矩阵B"<<endl;
    for(i=0;i<a;i++)    //给B矩阵赋值
    {
        for(j=0;j<a;j++)
        {
            cin>>B[i][j];
        }
    }
    cout<<"您输入的第二个矩阵B为："<<endl;
    for(i=0;i<a;i++)                                               
    {
        for(j=0;j<a;j++)
        {
            cout<<'\t';
            if((j+1)%a==0)
                cout<<B[i][j]<<endl;
            else
                cout<<B[i][j];
        }

    }

    for(i=0;i<a;i++)     //矩阵乘法
    {
        for(j=0;j<a;j++)
        {
            for(k=0;k<a;k++)
            { 
                M+=A[i][k]*B[k][j];
            }
            C[i][j]=M;
            M=0;
        }
    }

    cout<<"两个矩阵相乘得 "<<endl;
    for(i=0;i<a;i++)
    {
        for(j=0;j<a;j++)
        {
            cout<<'\t';
            if((j+1)%a==0)
                cout<<C[i][j]<<endl;
            else
                cout<<C[i][j];
        }
    }
    for(i = 0;i<a;i++)
    {
        delete []A[i];
    }
    delete []A;

    for(i = 0;i<a;i++)
    {
        delete []B[i];
    }
    delete []B;

    for(i = 0;i<a;i++)
    {
        delete []C[i];
    }
    delete []C;
    return 0;
}