代码

typedef unsigned int u32;
void matrix_multiply(int n,const int *A,const int *B,int *C){
	for(int i=0;i<n;++i){
		for(int k=0;k<n;++k){
			for(int j=0;j<n;j+=8){
				C[i*n+j]+=(u32)A[i*n+k]*B[k*n+j];
				C[i*n+j+1]+=(u32)A[i*n+k]*B[k*n+j+1];
				C[i*n+j+2]+=(u32)A[i*n+k]*B[k*n+j+2];
				C[i*n+j+3]+=(u32)A[i*n+k]*B[k*n+j+3];
				C[i*n+j+4]+=(u32)A[i*n+k]*B[k*n+j+4];
				C[i*n+j+5]+=(u32)A[i*n+k]*B[k*n+j+5];
				C[i*n+j+6]+=(u32)A[i*n+k]*B[k*n+j+6];
				C[i*n+j+7]+=(u32)A[i*n+k]*B[k*n+j+7];
			}
		}
	}
}

评测结果