C ++ для вложенных циклов — внутренний продукт — умножение матрицы

У меня есть матрица A с 4 строками и 3 столбцами и матрица B с 8 строками и 4 столбцами. Коэффициенты в первом столбце 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;
}