C++嵌套用于循环-内积-marix乘法

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

#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;
}