Вычислить сумму матриц E = A + A*A + А*А*А+. до степени к

# include <bits/stdc++.h>

using namespace std;

const int N = 100;

int a [N][N], b [N][N], n;

void Fill (int [][N], int);
void Print (int [][N], int, const char * = "");
void Mul (int [][N], int [][N], int [][N], int); // a=b*c
void Sum (int [][N], int [][N], int [][N], int); // a=b+c
void Copy (int [][N], int [][N], int); // a=b
void Task03 (int [][N], int [][N], int, int);

int main ()
{
setlocale (LC_ALL, "rus");

cin >> n;

Fill (a, n); Print (a, n, "A");

Task03 (b, a, k, n);

Print (b, n, "A+A*A+A*A*A");


return 0;
}

void Task03 (int b [][N], int a [][N], int n, int k)
{
;
}

void Mul (int a [][N], int b [][N], int c [][N], int n) // a=b*c
{
for ( int i=0; i < n; i++ )
for ( int j=0; j < n; j++ )
{
a [i][j] = 0;
for ( int k=0; k < n; k++ )
a [i][j] += b [i][k] * c [k][j];
}
}

void Sum (int a [][N], int b [][N], int c [][N], int n)// a=b+c
{
for ( int i=0; i < n; i++ )
for ( int j=0; j < n; j++ )
a [i][j] = b [i][j] + c [i][j];
}

void Copy (int a [][N], int b [][N], int n) // a=b
{
for ( int i=0; i < n; i++ )
for ( int j=0; j < n; j++ )
a [i][j] = b [i][j];
}

void Fill (int a [][N], int n)
{
for ( int i=0; i < n; i++ )
for ( int j=0; j < n; j++ )
a [i][j] = rand () % 2;
}

void Print (int a [][N], int n, const char * title)
{
cout << title << endl;
for ( int i=0; i < n; i++ )
{
for ( int j=0; j < n; j++ )
cout << setw (3) << a [i][j];
cout << endl;
}
cout << endl;
}

0