#include <stdio.h>

#define MAXDIM 50

typedef struct {
   unsigned int r;          /* row dimension */
   unsigned int c;          /* column dimension */
   int elt[MAXDIM][MAXDIM]; /* 2-d array of the elements of the matrix */
} matrix;

void printMat ( matrix M )
{
}

int printMatZ ( matrix M, int neltp )
{
   matrix N;
   int i,j,n;

   /* No recursive calls are made for 1 x 1 matrices */
   if (M.r == 1) {
      printf("%5d", M.elt[0][0]);
      if (++neltp == 15) { printf("\n"); neltp = 0; }
      return neltp;
   }

   /* Set the dimensions of the submatrices */
   n = (M.r) / 2;
   N.r = N.c = n;

   /* Recursively print the top left submatrix */
   for (i=0; i<n; ++i) for (j=0; j<n; ++j) N.elt[i][j] = M.elt[i][j];
   neltp = printMatZ(N,neltp);

   /* Recursively print the top right submatrix */
   for (i=0; i<n; ++i) for (j=0; j<n; ++j) N.elt[i][j] = M.elt[i][n+j];
   neltp = printMatZ(N,neltp);

   /* Recursively print the bottom left submatrix */
   for (i=0; i<n; ++i) for (j=0; j<n; ++j) N.elt[i][j] = M.elt[n+i][j];
   neltp = printMatZ(N,neltp);

   /* Recursively print the bottom right submatrix */
   for (i=0; i<n; ++i) for (j=0; j<n; ++j) N.elt[i][j] = M.elt[n+i][n+j];
   neltp = printMatZ(N,neltp);

   return neltp;
}

matrix initMatRec ( unsigned int k, unsigned int *cnt )
{
}

matrix initMatItr ( unsigned int k )
{
}

int main ()
{
   matrix M, N;
   unsigned int k;
   unsigned int cnt = 0;

   printf("k = "); scanf("%u", &k); printf("%u\n", k);

   M = initMatRec(k,&cnt);
   printf("\nRow-major listing of M:\n");
   printMat(M);
   printf("\nZ-listing of M:\n");
   printMatZ(M,0);
   printf("\n");

   N = initMatItr(k);
   printf("\nRow-major listing of N:\n");
   printMat(N);
   printf("\nZ-listing of N:\n");
   printMatZ(N,0);
   printf("\n");
}