C++实现矩阵乘法

#include<iostream>

using namespace std;

void print(int *p,int a,int b){
for(int i=0;i<a;i++){
for (int j = 0; j < b; j++)
{
cout<<p[i*b+j]<<"\t";
}
cout<<endl;
}
}

void mul(int *p,int *q,int a,int b,int c,int *r){
for (int i = 0; i < a; i++){
for (int j = 0; j < c; j++){
int temp=0;
for(int k=0;k<b;k++){
temp=temp+p[i*b+k]*q[k*c+j];//注意p和q两个矩阵的列数是不一样的，在这里debug了好久。T_T
}
r[i*c+j]=temp;
}
}
}

int main(){
int a,b,c,d;
int *p,*q;
cin>>a>>b;
while(a!=0){
p=new int[a*b];
for (int i = 0; i < a*b; i++){
cin>>p[i];
}
print(p,a,b);

cin>>c>>d;
q=new int[c*d];
for (int i = 0; i < c*d; i++){
cin>>q[i];
}
print(q,c,d);
int *r=new int[a*d];
mul(p,q,a,b,d,r);

print(r,a,d);
cin>>a>>b;
}
}