дано число p и квадратные матрицы A,B и C 4-го порядка. Получить (ABC)^p

#include <conio.h>
#include <iostream>
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
 
using namespace std;
 
int main()
{
  float **A,**B,**C,**D,**K;
   double p,n;
 
 
   cout<<"Please enter dimension matrix"<<endl;
   cin>>n;
 
   A=(float**)calloc(n,sizeof(float));
        for (int i=0; i<n; i++)
        {
        A[i]=(float*)calloc(n,sizeof(float));
      }
   B=(float**)calloc(n,sizeof(float*));
   for (int i=0; i<n; i++)
        {
        B[i]=(float*)calloc(n,sizeof(float));
      }
 
     C=(float**)calloc(n, sizeof(float*));
     for (int i=0; i<n; i++)
     {C[i]=(float*)calloc(n,sizeof(float));}
 
 
     D=(float**)calloc(n, sizeof(float*));
     for (int i=0; i<n; i++)
     {D[i]=(float*)calloc(n,sizeof(float));}
 
 
     K=(float**)calloc(n, sizeof(float*));
     for (int i=0; i<n; i++)
     {K[i]=(float*)calloc(n,sizeof(float));}
 
        for(int i=0;i<n;i++)
        {
            for(int j=0; j<n;j++)
            {
                cout<<"PLease enter A["<<i<<"]["<<j<<"] ="<<endl;
                cin>>A[i][j];
            }
 
        }
 
        for(int i=0;i<n;i++)
        {
            for(int j=0; j<n;j++)
            {
                cout<<"PLease enter B["<<i<<"]["<<j<<"] ="<<endl;
                cin>>B[i][j];
 
            }
 
        }
        for(int i=0;i<n;i++)
        {
            for(int j=0; j<n;j++)
            {
                cout<<"PLease enter C["<<i<<"]["<<j<<"] ="<<endl;
                cin>>C[i][j];
 
            }
 
        }
 
 
 
        for(int i=0;i<n;i++)
        {
            for(int j=0; j<n;j++)
            {
                D[i][j]=0;
                for(int k=0;k<n;k++)
                {
                    D[i][j]+=A[i][k]*B[k][j];
                }
 
 
             }
        }
 
        for(int i=0;i<n;i++)
        {
            for(int j=0; j<n;j++)
            {
               cout<<"\n"<<D[i][j];
 
            }
        }
 
        for(int i=0;i<n;i++)
        {
            for(int j=0; j<n;j++)
            {
                K[i][j]=0;
                for(int k=0;k<n;k++)
                {
                    K[i][j]+=D[i][k]*C[k][j];
                }
 
 
             }
        }
 
        for(int i=0;i<n;i++)
        {
            for(int j=0; j<n;j++)
            {
               cout<<"\n"<<K[i][j];
 
            }
        }
 
 
}