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

int main(void)
{
	int i,j;
	int test[4321],joint;
	double corr;
	FILE *fp1,*fp2,*fp3;

	if((fp1=fopen("expcounts.txt", "r"))==NULL)
	{
		printf("Error opening file.\n");
		exit(1);
	}

	for (i=0; i<4321; i++)
		fscanf(fp1,"%*s %d",&test[i]);

	if(fclose(fp1)==EOF)
	{
		printf("Error closing file.\n");
		exit(1);
	}

	if((fp2=fopen("corr.txt", "r"))==NULL)
	{
		printf("Error opening file.\n");
		exit(1);
	}

	if((fp3=fopen("correl.txt", "w"))==NULL)
	{
		printf("Error opening file.\n");
		exit(1);
	}

	for (i=0; i<4321; i++)
	{
		for (j=0; j<4321; j++)
		{
			fscanf(fp2,"%d",&joint);
			corr=(i==j) ? 1.0: (((double)(joint)/10295.0)-((double)(test[i])/10295.0)*((double)(test[j])/10295.0))/sqrt(((double)(test[i])/10295.0)*((10295.0-(double)(test[i]))/10295.0)*((double)(test[j])/10295.0)*((10295.0-(double)(test[j]))/10295.0));
			fprintf(fp3,"%lf ",corr);
		}
		fprintf(fp3,"\n");
	}

	if(fclose(fp2)==EOF)
	{
		printf("Error closing file.\n");
		exit(1);
	}

	if(fclose(fp3)==EOF)
	{
		printf("Error closing file.\n");
		exit(1);
	}
	
	return 0;
}