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;
	}
}