#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<errno.h>
#include<sys/stat.h>
 
off_t fsize(const char *filename);
unsigned char max(unsigned char item,unsigned char mx);
unsigned char min(unsigned char item,unsigned char mn);
 
int main()	{
	int **sums;
	int **counts;
	float **diff;
	float **promedios;
	float **desviacion;
	float *sumdesviaciones;
	float *promedioDesviaciones;
	float *diffdesviaciones;
	float *desviacion2;
	FILE *file;
	char *buffer;
	float temp;
	int l,i,j,k;
	file = fopen("salida.txt","r");
	if(file)	{
		l = fsize("salida.txt");
		buffer = calloc(l+1,sizeof(char));
		fread(buffer,sizeof(char),l,file);
		fclose(file);
		sums = calloc(l,sizeof(int*));
		counts = calloc(l,sizeof(int*));
		diff = calloc(l,sizeof(float*));
		promedios = calloc(l,sizeof(float*));
		desviacion = calloc(l,sizeof(float*));
		sumdesviaciones = calloc(l,sizeof(float));
		promedioDesviaciones = calloc(l,sizeof(float));
		diffdesviaciones = calloc(l,sizeof(float));
		desviacion2 = calloc(l,sizeof(float));
		i = 0;
 
		while(i < l)	{
			sums[i] = calloc(i+1,sizeof(int));
			counts[i] = calloc(i+1,sizeof(int));
			diff[i] = calloc(i+1,sizeof(float));
			promedios[i] = calloc(i+1,sizeof(float));
			desviacion[i] = calloc(l,sizeof(float));
			i++;
		}
 
		//rellenar array de sumas
		j = 0;
		while(j < l)	{
			k =0;
			while(k < l)	{
				//printf("Procesando sums[%i][%i]\n",j,(k % (j+1)));
				sums[j][k % (j+1)] += buffer[k];
				counts[j][k % (j+1)] +=1;
				k++;
			}
			j++;
		}
		//mostramos las sumas y los contadores
		i=1;
		j=0;
		while(j < l)	{
			k =0;
			while(k < i)	{
				//printf("sums[%i][%i] = %i\n",j,k,sums[j][k]);
				//printf("counts[%i][%i] = %i\n",k,i,counts[j][k]);
				k++;
			}
			i++;
			j++;
		}
		//calcular promedios
		i=1;
		j=0;
 
		while(j < l)	{
			k =0;
			while(k < i)	{
				promedios[j][k] = ((float)sums[j][k])/((float)counts[j][k]);				
				//printf("promedios[%i][%i] = %f\n",j,k,promedios[j][k]);
				k++;
			}
			i++;
			j++;
		}
		//calcular Diferencias
		j = 0;
		while(j < l)	{
			k =0;
			while(k < l)	{
				temp = promedios[j][k % (j+1)] - buffer[k];
				if(temp < 0.0)	{
					temp = -1.0 *temp;
				}
				diff[j][k % (j+1)] += temp;
				k++;
			}
			j++;
		}
		//calcular desviacion
		i=1;
		j=0;
		while(j < l)	{
			k =0;
			while(k < i)	{
				desviacion[j][k] = diff[j][k]/counts[j][k];				
				//printf("desviacion[%i][%i] = %f\n",j,k,desviacion[j][k]);
				k++;
			}
			i++;
			j++;
		}
 
		//calcular suma desviaciones
		i=1;
		j=0;
		while(j < l)	{
			k =0;
			while(k < i)	{
				sumdesviaciones[j]+= desviacion[j][k];
				k++;
			}
			promedioDesviaciones[j] =  sumdesviaciones[j]/i;
			printf("promedio desviacion[%i] = %f\n",j,promedioDesviaciones[j]);
			i++;
			j++;
		}
		//calcular differencia de desviaciones
		i=1;
		j=0;
		while(j < l)	{
			k =0;
			while(k < i)	{
				diffdesviaciones[j]+= promedioDesviaciones[j] - desviacion[j][k];
				k++;
			}
			desviacion2[j] = diffdesviaciones[j]/(float)k;
			printf("desviacion2[%i] = %f\n",j,desviacion2[j]);
			i++;
			j++;
		}
 
		printf("l : %i\n",l);
	}
	return 0;
}
 
unsigned char max(unsigned char item,unsigned char mx)	{
	return (item >= mx)? item:mx;
}
 
unsigned char min(unsigned char item,unsigned char mn)	{
	return (item <= mn)? item:mn;
}
 
off_t fsize(const char *filename) {
	struct stat st;
	if (stat(filename, &st) == 0)
		return st.st_size;
	fprintf(stderr, "Cannot determine size of %s: %s\n",filename, strerror(errno));
	return -1;
}