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

#define Pi 3.14159
#define FITNESSMINIMUM -1000.0
	double maxfitness=-DBL_MAX;


int	geneindex=0;
int globalshowflag=1;
int flag=0;

#define V 10
#define M 1
#define RANDMAX 32768.0 
double Xn[V];


double Qual(int g)
{

	return -fabs(fabs(Xn[0])+fabs(1.0-Xn[1])+fabs(2.0-Xn[2]));

}
void evolution(void)
{
	int round;
	double Qn,Dn;
	double xn,bestdn;
	round=0;
	Dn=0.3;
	double A1=1.3; int L=18, G=3000;
	double Xe[V][M], Xb[V][M], De[M], Db[M], Qb[M];
	int l,g,m,r,j,v;
	double z,h;
	double average;
	average=0.0;
	srand(5678765);
	for (m=0;m<M;m++){
		De[m]=Dn;
		Qb[m]=-DBL_MAX;
		
		for (v=0;v<V;v++){
			Xe[v][m]=sqrt(-2*log(1-(rand()/(RANDMAX+1.0)))/V)*sin(6.2832*(rand()/(RANDMAX+1.0)));
			
		}
		
	}
	maxfitness=-DBL_MAX;
	for (g=1;g<=G;g++){
		average=0.0;
		for (m=0;m<M;m++){
			Qb[m]=-DBL_MAX;
		}
		for (l=0;l<L;l++){
			for (Dn=0,m=0;m<M;m++) 
				Dn+=De[m]/(double)M;
			
			Dn*=pow(A1, 2*floor((rand()/(RANDMAX+1.0))+0.5)-1.0);
			
			g=g;
			for (v=0;v<V;v++){
				r=(int)((double)M*((double)rand()/(double)(RANDMAX+1.0)));
				z=(double)sqrt(-2.0*log(1.0-(rand()/(RANDMAX+1.0)))/V)*sin(6.2832*(rand()/(RANDMAX+1.0)));
				
				Xn[v]=Xe[v][r]+Dn*z;
			}	
			
			round++;
			Qn=Qual(g);
			average+=Qn;
			for (j=0,h=Qb[0],m=1;m<M;m++){
				if (Qb[m] < h){
					h=Qb[m];
					j=m;
				}
			}
			if (Qn> Qb[j]){
				Qb[j]=Qn;
				Db[j]=Dn;
				bestdn=Dn;
				for (v=0;v<V;v++){
					Xb[v][j]=Xn[v];
				}
			}
		}
		for (m=0;m<M;m++){
			De[m]=Db[m];
			for (v=0;v<V;v++){
				Xe[v][m]=Xb[v][m];
			}
		}
		maxfitness=-DBL_MAX;
		
		
	}
	for (v=0;v<V;v++){
		for (Xn[v]=0,m=0;m<M;m++){
			Xn[v]+=Xe[v][m]/M;
		}
	}
	for (Dn=0,m=0;m<M;m++){
		Dn+=De[m]/(double)M;
	}
	printf("%i, %e %e \n",G,Qual(G),Dn);
}





void main(void)
{
	evolution();
}