-3

我有一个 4 行 3 列的矩阵 A 和一个 8 行 4 列的矩阵 B。B 的第一列中的系数应该表示我想用 B 的行作为内积的 A 的那些行。

#include <iostream>
#include <vector>
#include <iomanip>
using namespace std;

int main()
{
int ik, il, im;
//A Matrix
vector<vector<int>> A(4, vector<int>(3));
A={ {0, 0, 0}, {1, 1, 0}, {0, 0, 1}, {1, 1, 1} };

//B Matrix
vector<vector<int>> B(8, vector<int>(4));
B={ {1, 0, 0, 0}, {2, 1, 0, 0}, {2, 0, 1, 0}, {2, 1, 1, 0} ,{3, 0, 0, 1}, {4, 1, 0, 1}, {4, 0, 1, 1}, {4, 1, 1, 1} };

vector<int> BZI; 
BZI={{1},{2},{2},{2},{3},{4},{4},{4}};

//calculate inner products
vector<vector<int>> SKP_AB(4, vector<int>(8));
int skp_temp;
for(ik=0; ik<4; ik++)
{
  for(im=0; im<8; im++)
  {
    for(il=0; il<3; il++)
    {
    skp_temp=A[BZI[ik]][il]*B[im][il+1];
    SKP_AB[ik][im]+=skp_temp;
    }
  }
}
//Ausgabe von SKP_AB
cout << "\n" << "#SKP_AB" << "\n";
for(ik=0; ik<4; ++ik)
{
 for(il=0; il<8; ++il)
 {
 cout << setw(2) << SKP_AB[ik][il];
 }
cout << "\n";
}

return 0;
}
4

1 回答 1

0

在 BZI 和 B 的第一列中,您使用基于行的行索引,但 C++ 使用基于零的行索引。

所以你需要补偿,你有BZI[ik]你需要的地方BZI[ik]-1

另请注意,在不被告知您期望的输出(只有您的程序不工作)的情况下达到该答案的难度与曾经告知您期望的输出有多么明显。

于 2015-11-16T20:05:10.197 回答