#include <stdio.h>
#include <stdlib.h>

int * makeVector(int n)
{  return (int *)malloc(n*sizeof(int));  }

int ** makeMatrix(int m,int n)
{  int ** newm = (int **)malloc(m*sizeof(int *));
   int i;
   for(i=0;i<m;i++)
    newm[i] = makeVector(n);
   return newm;
}

readMatrix(FILE * fin,int * * M,int m,int n)
{  int i,j;

   for(i=0;i<m;i++)
    for(j=0;j<n;j++)
     fscanf(fin,"%d",&(M[i][j]));
}

writeMatrix(FILE * fout,int * * M,int m,int n)
{  int i,j;

   for(i=0;i<m;i++)
   {  for(j=0;j<n;j++)
       fprintf(fout,"%3d ",M[i][j]);
      putc('\n',fout);
   }
}

int max(int a[],int n)
{  int m,i;
   m = 0;
   for(i=1;i<n;i++)
   if(a[i]>a[m])
    m = i;
   return m;
}

writeMax(char * s,int a[],int n)
{  int i,m;
   m = max(a,n);
   for(i=0;i<n;i++)
    printf("%d ",a[i]);
   printf("- max = %s[%d] = %3d\n",s,m,a[m]);
}

main(int argc,char ** argv)
{  FILE * fin;
   int ** M;
   int i,j,m,n;
   int * subjectav,* trialav;

   if((fin = fopen(argv[1],"r"))==NULL)
   {  printf("can't open %s\n",argv[1]);
      exit(0);
   }
   fscanf(fin,"%d %d",&m,&n);
   M = makeMatrix(m,n);
   readMatrix(fin,M,m,n);
   fclose(fin);
   writeMatrix(stdout,M,m,n);
   subjectav = makeVector(m);
   trialav = makeVector(n);
   for(i=0;i<m;i++)
   {  subjectav[i] = 0;
      for(j=0;j<n;j++)
       subjectav[i] = subjectav[i]+M[i][j];
      subjectav[i] = subjectav[i]/n;
   }
   for(i=0;i<n;i++)
   {  trialav[i] = 0;
      for(j=0;j<m;j++)
       trialav[i] = trialav[i]+M[j][i];
      trialav[i] = trialav[i]/m;
   }
   writeMax("Subject",subjectav,m);
   writeMax("Trial",trialav,n);  
}