Возведение матрицы в степень

#include <iostream>
#include <Windows.h>
#include <fstream>
#include <deque>
#include <algorithm>
 
using namespace std;
 
int** Multiply(int** a, int** b, int n);
int** MatrixPower(int** a, int n, int m);
int **SumMatrix(int**a, int**b, int n);
int main() {
    SetConsoleCP(1251);
    SetConsoleOutputCP(1251);
 
    ifstream in("Matrix.txt");
    int n;
    in >> n;
 
    int **adjacencyMatrix = new int*[n]; // матрица смежности
    for (int i = 0; i < n; i++) {
        adjacencyMatrix[i] = new int[n];
    }
 
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; j++) {
            in >> adjacencyMatrix[i][j];
        }
    }
 
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; j++) {
            cout << ' ' << adjacencyMatrix[i][j];
        }
        cout << '\n';
    }
    cout << "\n\n";
 
 
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; j++) {
            cout << ' ' << adjacencyMatrix[i][j];
        }
        cout << '\n';
    }
    cout << "\n\n";
 
    int **firstMat = new int*[n]; // единичная матрица
    for (int i = 0; i < n; i++) {
        firstMat[i] = new int[n];
    }
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; j++) {
            if (i == j) {
                firstMat[i][j] = 1;
            }
            else firstMat[i][j] = 0;
        }
    }
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; j++) {
            cout << ' ' << firstMat[i][j];
        }
        cout << '\n';
    }
    cout << "\n\n";
 
        for (int h = 1; h <= n; h++) {
        cout << h << endl;
        adjacencyMatrix = MatrixPower(adjacencyMatrix, n, h);
        cout << n<<endl;
        for (int i = 0; i < n; i++) {
            for (int j = 0; j < n; j++) {
                cout << ' ' << adjacencyMatrix[i][j];
            }
            cout << '\n';
        }
        cout << "\n\n";
    }
 
    system("pause");
    return 0;
}
int** Multiply(int** a, int** b, int n)
{
    int sum;
    int** c = new int*[n];
    for (int i = 0; i < n; i++)
    {
        c[i] = new int[n];
    }
    for (int j = 0; j < n; j++)
    {
        for (int i = 0; i < n; i++)
        {
            sum = 0;
            for (int r = 0; r < n; r++)
            {
                sum += a[j][r] * b[r][i];
            }
            c[j][i] = sum;
        }
    }
    return c;
    for (int i = 0; i < n; i++)
    {
        delete[] c[i];
    }
    delete[] c;
}
 
int** MatrixPower(int** a, int n, int m)
{
    if (m == 1)
        return a;
    if (m % 2 == 1)
        return Multiply(MatrixPower(a, n, m - 1), a, n);
    else {
        a = MatrixPower(a, n, m / 2);
        return Multiply(a, a, n);
    }
}