Разобраться в коде. Возведение в степень

#include <iostream>
#include <stdlib.h>
#include <iostream>
#include <vector>
 
using namespace std;
 
struct M
{
    int N;
    vector <vector<int>>m;
    M(){}
    M(int N) :N(N)
    {
        m.resize(N, vector<int>(N, 0));
    }
 
    M operator *(const M & b)
    {
        M c(N);
        for (int i = 0; i < N; i++)
        {
 
            for (int j = 0; j < N; j++)
            {
                for (int s = 0; s < N; s++)
                {
                    c.m[i][j] =(c.m[i][j]+m[i][s] * b.m[s][j])%10000;
                }
            }
        }
        return c;
    }
 
};
M b_pow(M m, int n) {
    M c(m.N);
    for (int i = 0; i < m.N; i++)
        c.m[i][i] = 1;
 
    while (n){
        if (n & 1){
            n --;
            c = c*m;
        }
        m = m*m;
        n >>= 1;
    }
    return c;
}
int main()
{
    int n;
    cin >> n;
    M a(n);
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            cin >> a.m[i][j];
        }
    }
    int s;
    cin >> s;
    M b = b_pow(a, s);
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            cout<< b.m[i][j]<<" ";
        }
        cout << endl;
    }
}