#include <stdio.h>
#include <iostream>
#include <stdlib.h>
#include <math.h>
#include <time.h>
using namespace std;
int Perevod(int m[10])
{
    int rez=0;
    for(int i=9;i>=0;i--)
        rez+=m[i]*pow(2,9-i);
    return rez;
};
int main(int argc, char *argv[])
{
    FILE *f;
    int mas[40][20],flag=0,kol=0,par[10],par2[10],q,x,y,child[40][20],p1,p2;
    float mas3[40],mas2[40][2],q1;
    f=fopen("d:\\siit.txt","w");
    srand(time(NULL));
    for(int i=0;i<20;i++)
    {
        mas[0][i]=rand()%2;
    }
    for(int i=1;i<40;i++)
    {
        //flag=0;
        do
        {
            for(int j=0;j<20;j++)
                mas[i][j]=rand()%2;
            for(int j=0;j<i;j++)
            {
                kol=0;
                for(int k=0;k<20;k++)
                    if(mas[i][k]==mas[j][k])kol++;
                if(kol==20)flag=1;
            }
        }
        while(flag);
    }
    for(int i=0;i<40;i++)
    {
        for(int j=0;j<10;j++)
        {
            par[j]=mas[i][j];
            par2[j]=mas[i][j+10];
        }
        x=Perevod(par);
        y=Perevod(par2);
        mas2[i][0]=x*2/1024;
            mas2[i][1]=y*2/1024;
    }
    for(int i=0;i<40;i++)
mas3[i]=mas2[i][1]*mas2[i][1]*exp(-(mas2[i][1]* mas2[i][1]))* mas2[i][1]*mas2[i][1]*exp(-(mas2[i][1]* mas2[i][1]));
    for(int i=0;i<39;i++)
        for(int j=0;j<39;j++)
            if(mas3[j]<mas3[j+1])
            {
                for(int k=0;k<20;k++)
                {
                    q=mas[j][k];
                    mas[j][k]=mas[j+1][k];
                    mas[j+1][k]=q;
                }
                q1=mas3[j];
                mas3[j]=mas3[j+1];
                mas3[j+1]=q1;
                q1=mas2[j][0];
                mas2[j][0]=mas2[j+1][0];
                mas2[j+1][0]=q1;
                q1=mas2[j][1];
                mas2[j][1]=mas2[j+1][1];
                mas2[j+1][1]=q1;
            }
    for(int i=0;i<40;i++)
        fprintf(f,"%f    %f %f\n",mas3[i],mas2[i][0],mas2[i][1]);
    for(int k=0;k<20;k++)
    {
        for(int i=0;i<40;i++)
        {
            p1=rand()%20;
            do
            {
                p2=rand()%20;
            }
            while(p1==p2);
            for(int j=0;j<20;j++)
            {
                q=rand()%2;
                if(q==0)child[i][j]=mas[p1][j];
                else child[i][j]=mas[p2][j];
            }
        }
        for(int i=0;i<40;i++)
            for(int j=0;j<20;j++)
                mas[i][j]=child[i][j];
        for(int i=0;i<40;i++)
        {
            for(int j=0;j<10;j++)
            {
                par[j]=mas[i][j];
                par2[j]=mas[i][j+10];
            }
            x=Perevod(par);
            y=Perevod(par2);
            mas2[i][0]=x*2/1024;
            mas2[i][1]=y*2/1024;
        }
        for(int i=0;i<40;i++)
            mas3[i]=mas2[i][1]*mas2[i][1]*exp(-(mas2[i][1]* mas2[i][1]))* mas2[i][1]*mas2[i][1]*exp(-(mas2[i][1]* mas2[i][1]));
        for(int i=0;i<39;i++)
            for(int j=0;j<39;j++)
                if(mas3[j]<mas3[j+1])
                {
                    for(int k=0;k<20;k++)
                    {
                        q=mas[j][k];
                        mas[j][k]=mas[j+1][k];
                        mas[j+1][k]=q;
                    }
                    q1=mas3[j];
                    mas3[j]=mas3[j+1];
                    mas3[j+1]=q1;
                    q1=mas2[j][0];
                    mas2[j][0]=mas2[j+1][0];
                    mas2[j+1][0]=q1;
                    q1=mas2[j][1];
                    mas2[j][1]=mas2[j+1][1];
                    mas2[j+1][1]=q1;
                }
        fprintf(f,"%f\n",mas3[0]);
    }
    for(int i=0;i<40;i++)
        fprintf(f,"%f    %f %f\n",mas3[i],mas2[i][0],mas2[i][1]);
    fclose(f);
    return 0;
}