#include<bits/stdc++.h>
using namespace std;
#define N 8
int i,j;
int A[N][N],B[N][N],C[N][N];
void input(int n,int x[][N]){
for(i=0;i<n;i++)
for(j=0;j<n;j++){
cin>>x[i][j];
}
}
void output(int n,int x[][N]){
for(i=0;i<n;i++){
for(j=0;j<n;j++){
cout<<x[i][j]<<' ';
}
cout<<endl;
}
}
void Matrix_M(int n,int A[][N],int B[][N],int C[][N]){
for(i=0;i<n;i++)
for(j=0;j<n;j++){
C[i][j]=0;
for(int k=0;k<n;k++)
C[i][j]+=A[i][k]*B[k][j];
}
}
void Matrix_A(int n,int A[][N],int B[][N],int C[][N]){
for(i=0;i<n;i++)
for(j=0;j<n;j++){
C[i][j]=A[i][j]+B[i][j];
}
}
void Matrix_S(int n,int A[][N],int B[][N],int C[][N]){
for(i=0;i<n;i++)
for(j=0;j<n;j++){
C[i][j]=A[i][j]-B[i][j];
}
}
void Strassen(int n,int A[][N],int B[][N],int C[][N]){
int A11[N][N],A12[N][N],A21[N][N],A22[N][N];
int B11[N][N],B12[N][N],B21[N][N],B22[N][N];
int C11[N][N],C12[N][N],C21[N][N],C22[N][N];
int M1[N][N],M2[N][N],M3[N][N],M4[N][N],M5[N][N],M6[N][N],M7[N][N];
int T1[N][N],T2[N][N];
if(n==2){
Matrix_M(2,A,B,C);
return;
}
for(i=0;i<n/2;i++){
for(j=0;j<n/2;j++){
A11[i][j]=A[i][j];
A12[i][j]=A[i][j+n/2];
A21[i][j]=A[i+n/2][j];
A22[i][j]=A[i+n/2][j+n/2];
B11[i][j]=B[i][j];
B12[i][j]=B[i][j+n/2];
B21[i][j]=B[i+n/2][j];
B22[i][j]=B[i+n/2][j+n/2];
}
}
//M1=A11(B12-B22)
Matrix_S(n/2,B12,B22,T1);
Strassen(n/2,A11,T1,M1);
//M2=(A11+A12)B22
Matrix_A(n/2,A11,A12,T1);
Strassen(n/2,T1,B22,M2);
//M3=(A21+A22)B11
Matrix_A(n/2,A21,A22,T1);
Strassen(n/2,T1,B11,M3);
//M4=A22(B21-B11)
Matrix_S(n/2,B21,B11,T1);
Strassen(n/2,A22,T1,M4);
//M5=(A11+A22)(B11+B22)
Matrix_A(n/2,A11,A22,T1);
Matrix_A(n/2,B11,B22,T2);
Strassen(n/2,T1,T2,M5);
//M6=(A12-A22)(B21+B22)
Matrix_S(n/2,A12,A22,T1);
Matrix_A(n/2,B21,B22,T2);
Strassen(n/2,T1,T2,M6);
//M7=(A11-A21)(B11+B12)
Matrix_S(n/2,A11,A21,T1);
Matrix_A(n/2,B11,B12,T2);
Strassen(n/2,T1,T2,M7);
//C11=M5+M4-M2+M6
Matrix_A(n/2,M5,M4,T1);
Matrix_S(n/2,T1,M2,T2);
Matrix_A(n/2,T2,M6,C11);
//C12=M1+M2
Matrix_A(n/2,M1,M2,C12);
//C21=M3+M4
Matrix_A(n/2,M3,M4,C21);
//C22=M5+M1-M3-M7
Matrix_A(n/2,M5,M1,T1);
Matrix_S(n/2,T1,M3,T2);
Matrix_S(n/2,T2,M7,C22);
for(i=0;i<n/2;i++){
for(j=0;j<n/2;j++){
C[i][j]=C11[i][j];
C[i+n/2][j]=C21[i][j];
C[i][j+n/2]=C12[i][j];
C[i+n/2][j+n/2]=C22[i][j];
}
}
}
int main(){
int i,j,n;
cout<<"输入矩阵阶数："<<endl;
cin>>n;
cout<<"输入第一个矩阵："<<endl;
input(n,A);
cout<<"输入第二个矩阵："<<endl;
input(n,B);
Strassen(n,A,B,C);
cout<<"结果为："<<endl;
output(n,C);
return 0;
}
OJ.rar （29个子文件） 
OJ 
蓄栏保留.cpp 966B 求解迷宫问题.cpp 665B n皇后.cpp 900B Strassen矩阵乘法.cpp 3KB 加工生产调度.cpp 780B 分形宇宙.cpp 654B 智力大冲浪.cpp 955B 1268.cpp 390B Cow Party.cpp 1KB 田忌赛马.cpp 591B 真—回溯01背包.cpp 2KB 夺取宝藏.cpp 453B 食物链.cpp 821B 流水作业调度.cpp 798B n皇后摆放.cpp 898B 马走棋问题.cpp 387B 01背包.cpp 392B 1204.cpp 411B 饼干.cpp 355B 1719.cpp 747B 守望者.cpp 430B 最长公共子序列.cpp 460B 回溯01背包.cpp 993B 逆序对z.cpp 921B 单调不降.cpp 445B a+b+c+d==0.cpp 692B 前中转后.cpp 1KB wzy跑步.cpp 720B 部分背包.cpp 963B资源评论
