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

#define time_point 17

main()
{
 float exp_normX[time_point+1], exp_normY[time_point+1];
 int num_gen, i, j, k, m, n, len, found, startx, starty, pos_or_neg;
 float ave, std, tt, max_score, mat_pos[time_point+1][time_point+1], mat_neg[time_point+1][time_point+1];
 char dum[500], gene1[15], gene2[15], file_name[20];
 FILE *fp;


 printf("Which data file:\n");
 scanf("%s", file_name);

  if((fp = fopen(file_name, "r"))==NULL){
     printf("can't open file\n");
     exit(0);
     }

 num_gen = 0;
 while(!feof(fp)){
    fgets(dum, 500, fp);
    num_gen++;
    }
  num_gen--;

  rewind(fp);

  for(m=1; m<num_gen; m++){

    rewind(fp);
    for(k=0; k<=m; k++){
      fscanf(fp, "%s", gene1);
      for(i=1; i<=time_point; i++)
         fscanf(fp, "%f", &exp_normX[i]);
      }


   rewind(fp);
   for(n=0; n<m; n++){

      fscanf(fp, "%s", gene2);
      for(i=1; i<=time_point; i++)
         fscanf(fp, "%f", &exp_normY[i]);

   for(i=0; i<time_point; i++){
      mat_pos[0][i]=0.;
      mat_pos[i][0]=0.;
      mat_neg[0][i]=0.;
      mat_neg[i][0]=0.;
      }
   
   max_score = 0.;

   for(i=1; i<=time_point; i++)
      for(j=1; j<=time_point; j++){
           
           mat_pos[i][j] = mat_pos[i-1][j-1] + (exp_normX[i]*exp_normY[j]);
           if(mat_pos[i][j]<0.)
               mat_pos[i][j] = 0.;
 
           mat_neg[i][j] = mat_neg[i-1][j-1] - (exp_normX[i]*exp_normY[j]);
           if(mat_neg[i][j]<0.)
               mat_neg[i][j] = 0.;

        
           if(mat_pos[i][j]>max_score){
             max_score = mat_pos[i][j];
             startx = i;
             starty = j;
             pos_or_neg = 1;
             }

           if(mat_neg[i][j]>max_score){
             max_score = mat_neg[i][j];
             startx = i;
             starty = j;
             pos_or_neg = 0;
             }

           }


    if(pos_or_neg == 1){
      for(i=1; i<=time_point; i++)
        if(mat_pos[startx-i][starty-i]<=0.000001){
            break;
            }
      }
    else{
      for(i=1; i<=time_point; i++)
        if(mat_neg[startx-i][starty-i]<=0.000001){
            break;
            }
      }

    len = i;

printf("%s %s %f %d %d %d %d\n", gene1, gene2, max_score, startx, starty, len, pos_or_neg);

  }  //n

}   //m


}

fclose(fp);