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

int main(){

int **a,**b,**c;

int a_r,a_c,b_r,b_c,inc;
int i,j,k;
printf("\n insert number of rows and columns for a matrix:");
scanf("%d%d",&a_r,&a_c);
printf("\n insert number of columns for b matrix:");
scanf("%d",&b_c);
b_r=a_c;

/* alloca memoria per la matrice a */
a=(int **) malloc(sizeof(int*)*a_r);
for( i=0;i<a_r; i++) {
     a[i]=(int *) malloc(sizeof(int)*a_c);
     }

/* alloca memoria per la matrice b */
b=(int **) malloc(sizeof(int*)*b_r);
for( i=0;i<b_r; i++){
     b[i]=(int *) malloc(sizeof(int)*b_c);
     }
/* alloca memoria per la matrice c */
c=(int **) malloc(sizeof(int*)*a_r);
for( i=0;i<a_r; i++) {
     c[i]=(int *) malloc(sizeof(int)*b_c);
     }

for(i=0;i<a_r; i++) {
    for(j=0;j<a_c; j++) {
       a[i][j]= 2;
       }
}
for(i=0;i<b_r; i++) {
    for(j=0;j<b_c; j++) {
       b[i][j]= 3;
       }
}


/*inizializza c */
for(i=0;i<a_r; i++) {
    for(j=0;j<b_c; j++) {
       c[i][j]=0;
       }
}

 clock_t begin= clock();
/* matrix multiplication */
for(i=0;i<a_r; i++) {
    for(j=0;j<b_c; j++) {
       for(k=0;k<a_c; k++){
          if (i%2==0) {
            inc=4*pow((double)i,(double)i);
            c[i][j]=c[i][j]+a[i][k]*b[k][j] +inc;
            }
          else {
            inc=9*pow((double)i,(double)i);
            c[i][j]=c[i][j]+a[i][k]*b[k][j] +inc;
            }
          }
   }
}
 clock_t end = clock();
 double elapsed;
 elapsed = ((double) (end - begin)) / CLOCKS_PER_SEC;
 printf (" Time = %f \n",elapsed);


/*printf("\n Il prodotto matriciale di a e b vale\n");
for(i=0;i<a_r; i++) {
    for(j=0;j<b_c; j++) {
        printf("%d\t",c[i][j]);
        }
        printf("\n");
} */

/*libera memoria */
while (--a_r > -1) { free(a[a_r]); }
free(a);
while (--b_r > -1) { free(b[b_r]); }
free(b);
while (--a_r > -1) { free(c[a_r]); }
free(c);

return (0);
}