代码

#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 * n - 1) C[i] = 0;
	fo(i,0,n - 1) fo(k,0,n - 1)
	{
		if(A[idx(i,k)] == 0) continue;
		fo(j,0,n - 1) C[idx(i,j)] += A[idx(i,k)] * B[idx(k,j)];
	}
}

评测结果