#include <stdio.h>
#include <stdlib.h>
#define dim 100

void readMatrix(float mat[][dim],int row,int col);
float medianum(float mat[][dim],int fix_row,int col);
void sortRow(float mat[][dim],int fix_row,int col);
void sortMat(float mat[][dim],int fix_col,int row,int col);
void swap(float* a,float* b);

int main(int argv, char *argc[])
{
         FILE *fpin,*fpout;
         int r,c;
         int fix_col;
         float mat[dim][dim];
         float media[dim];
         
         //Apro il file input.txt
         fpin=fopen("input.txt","r");
         if(fpin!=NULL)
           {
           fscanf(fpin,"%d",&r); 
           fscanf(fpin,"%d\n\n",&c);
           
           for(int i=0;i<r;i++)
              for(int j=0;j<c;j++)
              {
              if(j==(c-1))
                fscanf(fpin,"%f\n",&mat[i][j]);
              else
                fscanf(fpin,"%f ",&mat[i][j]);
              }
           fclose(fpin);
           }
         else
           printf("Errore di apertura di input.txt.\n");
         
         readMatrix(mat,r,c);
        
         for(int i=0;i<r;i++)
            media[i]=medianum(mat,i,c);
         
         for(int i=0;i<r;i++)
            printf("\nMedia della riga %d-esima e': %6.2f\n\n",i,media[i]);
  
         //Aggiungo una colonna alla matrice con valore la corrispettiva media
         c=c+1; 
         //Ciclo che stampa la nuova matrice ed inserisce la media all'ultima colonna
         for(int i=0;i<r;i++)
            {
            for(int j=0;j<c;j++)
               { 
               if(j==(c-1))
                 {
                 mat[i][j]=media[i];
                 printf("%6.2f",mat[i][j]);
                 }
               else
                 printf("%6.2f",mat[i][j]);
               }
            printf("\n");
            }
          
         //Utilizzo il sottoprogramma che ordina la singola riga
         for(int i=0;i<r;i++)
            sortRow(mat,i,c);
         
         printf("\n\nOrdinato per riga:\n");
         for(int i=0;i<r;i++)
         {
            for(int j=0;j<c;j++)
               printf("%6.2f",mat[i][j]);
         printf("\n");
         }
         
         //Utilizzo il sottoprogramma che ordina le righe della matrice
         //per valori crescenti dell'ultima colonna
         fix_col=c-1;
         for(int i=0;i<r;i++)
            sortMat(mat,fix_col,r,c);
         
         //Stampo la matrice ordinata
         printf("\n\nMatrice ordinata:\n");
         for(int i=0;i<r;i++)
         {
           for(int j=0;j<c;j++)
             printf("%6.2f",mat[i][j]);
         printf("\n");
         }
         
         //Creo il file binario e ci salvo la matrice
         fpout=fopen("output.bin","w");
         if(fpout!=NULL)
           {
           fwrite(mat,sizeof(float),dim,fpout);
           fclose(fpout);
           printf("\n\nMatrice salvata in output.bin con successo.\n");
           }
         else
           printf("File output.bin non creato correttamente!\n");
}

//Funzione per lettura delle dimensioni e degli elementi della matrice
void readMatrix(float mat[][dim],int row,int col) {
           
           //Stampo a schermo il valore delle righe e delle col
           printf("La matrice ha:\nRighe: %d\ncol: %d\n\n",row,col);
           //Stampo la matrice
           for(int i=0;i<row;i++)
              {
              for(int j=0;j<col;j++)
                 printf("%6.2f",mat[i][j]);
              printf("\n");
              }
}

//Funzione per il calcolo della media per ogni riga
float medianum(float mat[][dim],int fix_row,int col) {
              
              float med=0;
              float sum=0;
                for(int i=0;i<col;i++)
                   sum+=mat[fix_row][i];
                med=sum/col;
              
              return med; //ritorna il valore della media della riga da posizionare all'interno del vettore
}
  
void sortRow(float mat[][dim],int fix_row,int col) {
               
  int k=col-1,scambio=1;
  
  while(k>0 && scambio==1)
  {
    scambio=0;
    for(int i=0;i<k;i++)
       if((mat[fix_row][i])>(mat[fix_row][i+1]))
       {
        swap(&mat[fix_row][i],&mat[fix_row][i+1]);
        scambio=1;
       }
    k--;
  }
}

void sortMat(float mat[][dim],int fix_col,int row,int col) {
  
  int k=row-1,scambio=1;
  while(k>0 && scambio==1)
  {
    scambio=0;
    for(int i=0;i<k;i++)
      if((mat[i][fix_col])>(mat[i+1][fix_col]))
      {
        scambio=1;
        for(int j=0;j<col;j++)
           swap(&mat[i][j],&mat[i+1][j]);
      }
  k--;
}
}
  
void swap(float* a,float* b) {
         
         float tmp;
         
         tmp=*a;
         *a=*b;
         *b=tmp;
}