克罗内克积的c++简单实现

// 函数实现

 vector<vector<int>> Kroneck(const vector<vector<int>> &a, const vector<vector<int>> &b, int i, int j, int m, int n)
{

    vector<vector<int>> c(i*m,vector<int>(j*n,0));
    //std::cout << c.size() << std::endl;

    int index_x = c.size();
    int index_y = c[0].size();

    for (int x = 0;x<index_x;x++)
    { 
        for (int y = 0; y < index_y; y++)
        {
            int a_i = x / m;
            int a_j = y / n;
            int b_i = x % m;
            int b_j = y % n;
            c[x][y] = a[a_i][a_j] * b[b_i][b_j];
        }
    }

    return c;
}

// 简单调用

    vector<vector<int>> a = { {1,2,3},{4,5,6} };
    vector<vector<int>> b = { {1,2},{1,2} };

    vector<vector<int>> c = Kroneck(a, b,2,3,2,2);