#include <stdio.h>
#define NMAX 100

struct elemento
{
	float val;
	int pos[2];
	
	};

void leggimat(float m[NMAX][NMAX],int *r, int *c);
void stampamat(float m[NMAX][NMAX],int r, int c);
void vettori(float m[NMAX][NMAX],int r, int c, float v1[],float v2[], int *n);
void ordina(float v[], int n);
float prodscal(float v1[], float v2[], int n);
void creastruct(struct elemento vet[], int *n, float mat[NMAX][NMAX], int r, int c);


int main()
{
	float mat[NMAX][NMAX];
	float v1[NMAX], v2[NMAX];
	int r,c,n;
	int ps;
	struct elemento vettore[NMAX];
	
	
	leggimat(mat,&r,&c);
	stampamat(mat,r,c);
	vettori(mat,r,c,v1,v2,&n);
	ordina(v1,n);
	ordina(v2,n);
	
	
	ps = prodscal(v1,v2,n);
	int k;
	creastruct(vettore,&k,mat,r,c);
	
	 
	 FILE *fp=fopen("output.bin","w");
	  if(fp!=NULL)
	  {
		  fwrite(vettore,sizeof(struct elemento), k,fp);
		  fclose(fp);
		  }
	
	} 

void stampamat(float m[NMAX][NMAX],int r, int c)
{
	for(int i=0;i<r;i++)
	{
		printf("\n");
	 for(int j=0;j<c;j++)
	 {
		 printf("%f\t",m[i][j]);
		 }
	}
	}

void leggimat(float m[NMAX][NMAX],int *r, int *c)
{
	FILE *fp;
	fp=fopen("input.txt", "r");
	if(fp!=NULL)
	{
		fscanf(fp,"%d",r);
		fscanf(fp,"%d",c);
		for(int i=0;i<*r;i++)
		  for(int j=0;j<*c;j++)
		     fscanf(fp, "%f",&m[i][j]);
		
		fclose(fp);
		}
	
	}
	
void vettori(float m[NMAX][NMAX],int r, int c, float v1[],float v2[], int *n)
{
	int k=0;
	for(int i=0;i<r;i=i+2)
	  for(int j=0;j<c; j=j+2)
	    {
			v1[k]=m[i][j];
			k++;
		}
	k=0;
	for(int i=1;i<r;i=i+2){
	  for(int j=1;j<c; j=j+2)
	    {
			v2[k]=m[i][j];
			k++;
		}
	}
	
	  *n=k;
	}


void ordina(float v[], int n)
{
	for(int i=0;i<n;i++)
	  for(int j=0;j<n-1-i;j++)
	    if(v[j]>v[j+1])
	    {
			int temp=v[j];
			v[j]=v[j+1];
			v[j+1]=temp;
			
			}
	
	}
	
	
float prodscal(float v1[], float v2[], int n)
	{
		int ps=0;
		for(int i=0;i<n;i++)
			ps=ps+v1[i]*v2[i];
		return ps;
		}
	
void creastruct(struct elemento vet[], int *n, float mat[NMAX][NMAX], int r, int c)
{
	 
	int k=0;
	for(int i=0;i<r;i=i+2)
	  for(int j=0;j<c; j=j+2)
	    {
			vet[k].val=mat[i][j];
			vet[k].pos[0]=i;
			vet[k].pos[1]=j;
			k++;
		}

	for(int i=1;i<r;i=i+2)
	  for(int j=1;j<c; j=j+2)
	    {
			vet[k].val=mat[i][j];
			vet[k].pos[0]=i;
			vet[k].pos[1]=j;
			k++;
		}
	*n=k;
	}