矩阵乘法c++

我想编写一个能够将2x2矩阵上升到k次方的程序。我知道如何使它平方，但一旦我想达到更高的权力，我努力保存我的结果，并在下一个方程中使用它。a b c d是矩阵中的数字，n是我想在一次内完成的矩阵的数量，k是我想对矩阵取的幂m是我想对这些数字取的模。我知道有一种方法可以使它相当简单，但我想不出一个好方法来在下一个方程中使用我之前做的方程的结果。

#include <iostream>
using namespace std;
int mnoz(int a, int b, int c,int d,int m){
    int a2 = (a*a + c*b) % m;
    int b2 = (a*b + b*d) % m;
    int c2 = (c*a + d*c) % m;
    int d2 = (c*b + d*d) % m;
    return a2, b2, c2, d2;
}
int main()
{
    int a, b, c, d, k, m, n;
    int e, f, g, h;
    int e2, f2, g2, h2;
    cin >> n;

    // a^k = (a^2)^k/2
    for (int i = 0; i<n; i++){
        cin >> a >> b >> c >> d >> k >> m;
        if (k == 1){
            cout << a << " " << b << " " << c << " " << d << endl;

        }

        else if (k == 2){
            e = (a*a + c*b) % m;
            f = (a*b + b*d) % m;
            g = (c*a + d*c) % m;
            h = (c*b + d*d) % m;
            cout << e << " " << f << " " << g << " " << h << endl;
        }
        else{
            if (k % 2 == 0){
                e = (a*a + c*b) % m;
                f = (a*b + b*d) % m;
                g = (c*a + d*c) % m;
                h = (c*b + d*d) % m;
                int z = (k/2)-1;
                for (int j = 0; j < z; j++){
                    int e2 = e;
                    int f2 = f;
                    int g2 = g;
                    int h2 = h;
                    mnoz(e2, f2, g2, h2, m);
                }
                cout << e << " " << f << " " << g << " " << h << endl;
            }
        }
    }
    system("pause");
    return 0;
}