マトリックス乗算テンプレートC/C++
1858 ワード
マトリクス乗算は多くのアルゴリズムの基礎であり,高速べき乗に乗算するように熟練する必要があるため,以下のテンプレートは個人的には使いやすいと思うが,時間の移行に伴って精進する必要があるに違いない.
コードの例:
完全なコードが必要かもしれません.
コードの例:
struct Matrix{
int n,m;
int v[maxn][maxn];
Matrix(int n,int m):n(n),m(m){}
void init(){ //
memset(v,0,sizeof v);
}
Matrix operator* (const Matrix B) const {
Matrix C(n,B.m); //
C.init();
for(int i = 0;i < n;i++)
for(int j = 0;j < B.m;j++)
for(int k = 0;k < m;k++)
C.v[i][j] += v[i][k]*B.v[k][j];
return C;
}
void print(){// ,
for(int i = 0;i < n;i++){
for(int j = 0;j < m;j++)
cout << v[i][j] << " ";
cout << endl;
}
}
};
完全なコードが必要かもしれません.
#include
#include
#include
using namespace std;
const int maxn = 105;
struct Matrix{
int n,m;
int v[maxn][maxn];
Matrix(int n,int m):n(n),m(m){}
void init(){ //
memset(v,0,sizeof v);
}
Matrix operator* (const Matrix B) const {
Matrix C(n,B.m); //
C.init();
for(int i = 0;i < n;i++)
for(int j = 0;j < B.m;j++)
for(int k = 0;k < m;k++)
C.v[i][j] += v[i][k]*B.v[k][j];
return C;
}
void print(){// ,
for(int i = 0;i < n;i++){
for(int j = 0;j < m;j++)
cout << v[i][j] << " ";
cout << endl;
}
}
};
int main(){
FILE *fp = fopen("Matrix.txt","r");
int n1,m1,n2,m2;
fscanf(fp,"%d%d",&n1,&m1);
Matrix A(n1,m1);
for(int i = 0;i < n1;i++)
for(int j = 0;j < m1;j++){
fscanf(fp,"%d",&A.v[i][j]);
}
fscanf(fp,"%d%d",&n2,&m2);
Matrix B(n2,m2);
for(int i = 0;i < n2;i++)
for(int j = 0;j < m2;j++){
fscanf(fp,"%d",&B.v[i][j]);
}
Matrix C = A*B;
C.print();
fclose(fp);
return 0;
}