Решение Диофантова уравнения через генетический алгоритм

#include <iostream>
#include <cstdlib>
#include <time.h>
#include <math.h>
 
using namespace std;
 
int main()
{
    setlocale(LC_ALL, "Russian");
    int a,b,c,d,e;
    e = 30;
    double revcof;
    srand(time(0));
    int sol[4][5];
    int ans[5];
    int diff[5];
    int par[5][2];
    int child[4][5];
    float percentpar[5];
    float chanse[1000];
    for(int f = 0;f<5;f++) // тут задаются первоначальные массивы чисел
    {
        for(int g =0;g<4;g++)
        {
            sol[g][f] = rand() % 30;
            cout << sol[g][f] << ",";
        }
        cout << endl;
    }
    while(diff[0] != 0 || diff[1] != 0 || diff[2] != 0 || diff[3] != 0 || diff[4] != 0) // цикл должен выполнятся до тех пор, пока разница между ответами и 30 не станет равной нулю.
    {
    for(int h = 0;h<5;h++) // Здесь вычисляется разность между ответами и 30, после чего вычисляется сумма обратных коэффициентов
    {
        ans[h] = sol[0][h] + 2 * sol[1][h] + 3 * sol[2][h] + 4 * sol[3][h];
        diff[h] = abs(ans[h] - 30);
        cout << "diff." << diff[h] << endl;
        revcof += 1.000000/diff[h];
    }
    for(int h = 0;h<5;h++) // Тут вычисляется процент шанса становления родителем, и округляется до 1 знака после запятой.
    {
        percentpar[h] = (1.00/diff[h])/revcof*100;
        percentpar[h] = round(percentpar[h]*10)/10;
        cout << percentpar[h] << endl; 
    }
    cout << "FINISH" << endl;
    for(int m = 0;m<5;m++) // В этих вложенных циклах грубо говоря генерируется рандом, который в дальнейшем и будет определять вероятность становления родителем.
    {
        int f = percentpar[m] * 10;
        for(int l = 0;l<f;l++)
        {
            int k;
            k = percentpar[m-1]*10;
            if(m == 0)
            {
            chanse[l] = 0;
            }
            if(m == 1)
            {
            chanse[l+k] = 1;
            }
            if(m == 2)
            {
            k = percentpar[m-1]*10 + percentpar[m-2]*10;
            chanse[l+k] = 2;
            }
            if(m==3)
            {
            k = percentpar[m-1]*10 + percentpar[m-2]*10 + percentpar[m-3]*10;
            chanse[l+k] = 3;
            }
            if(m==4)
            {
            k = percentpar[m-1]*10 + percentpar[m-2]*10 + percentpar[m-3]*10 + percentpar[m-4]*10;
            chanse[l+k] = 4;
            }
        }
    }
    for(int a=0;a<5;a++) // цикл по сути кидает игральную кость, с 1000 сторонами и шансом который он берёт из цикла предыдущего
    {
        int f = rand()%1001;
        par[a][0] = chanse[f];
        do // этот цикл не даёт возможности быть одновременно обоими родителями одному кандидату
        {
        int e = rand()%1001;
        par[a][1] = chanse[e];
        }
        while(par[a][1] == par[a][0]);
        cout << par[a][0] << "^" << par[a][1] << endl;
    }
    for(int b = 0;b<5;b++) // Ну а этот цикл уже отвечает за перестановку "генов" в хромосомах с рандомом, сгенерированным в прошлом цикле
    {
            int rannum1 = rand()%4 + 1;
        int gene = par[b][0];
    for(int a = 0;a<rannum1;a++)
    {
        sol[a][b] = sol[a][gene];
        //cout << sol[a][b] << ",";
    }
    int rannum2 = 4 - rannum1;
    gene = par[b][1];
    for(int c = 0; c < rannum2;c++)
    {
        sol[c][b] = sol[c][gene];
        //cout << sol[c][b] << ",";
    }
    int rannum3 = rand()%30 + 1; // Здесь создаётся случайная мутация в рандомной хромосоме
    int rannum4 = rand()%4 + 1;
    sol[rannum4][b] = rannum3;
    //cout << sol[rannum4][b] << endl;
    //cout << "rannum 4 = " << rannum4 << endl;
    }
    }
    cout << "Решение: a = " << sol[0][0] << "b = " << sol[1][0] << "c = " << sol[2][0] << "d = " << sol[3][0] << endl;
    return 0;
}