代码

#include <bits/stdc++.h>
using namespace std;
typedef unsigned int u32;
constexpr int K=32;
void matrix_multiply(int n,const int *A,const int *B,int *C){
	for(int ii=0;ii<n;ii+=8){
		for(int kk=0;kk<n;kk+=8){
			for(int jj=0;jj<n;jj+=8){
				int pa=ii*n+kk;
				int pb=kk*n+jj;
				int pc=ii*n+jj;
				for(int i=0;i<K;++i){
					for(int k=0;k<K;++k){
						for(int j=0;j<K;++j){
							C[pc+i*n+j]+=(u32)A[pa+i*n+k]*B[pb+k*n+j];
						}
					}
				}
			}
		}
	}
}

评测结果