代码

#define idx(i,j) i * n + j
#define fo(i,x,y) for(int i = x;i < y;++i)

void matrix_multiply(int n, const double *A, const double *B, double *C)
{
	fo(i,0,n) fo(k,0,n)
	{
		double s = A[idx(i,k)];
		fo(j,0,n) C[idx(i,j)] += s * B[idx(k,j)];
	}
}

评测结果