Задача о ранце

//Сам класс CDiophantine:
#include <stdlib.h>
#include <time.h>
#include <iostream>
#include <fstream>
#include <string>
 
#define MAXPOP  25
using namespace std;
class predmet
{
private:
    int w, p;
    string name;
public:
    predmet(){}
    void In(ifstream &ifst) 
    {
        ifst>>name;
        ifst>>w;
        ifst>>p;
    }
    void prin()
    {
        cout<<name<<" (Объем:"<<w<<" , Цена: "<<p<<" )\n";
    }
    int get_w(){return w;}
    int get_p(){return p;}
};
 
struct gene 
{
    int alleles[5];
    int fitness;
    float likelihood;
 
    // Проверка на равенство.
    int operator==(gene gn) 
    {
        for (int i=0;i<5;i++) 
        {
            if (gn.alleles[i] != alleles[i]) return false;
        }
 
        return true;
    }
};
 
class CDiophantine {
    public:
        CDiophantine(int, int, int, int, int, int);// Конструктор с коэффициентами A, B, C, D.
        int Solve();// Решить уравнение.
        
        //Возвращает данного гена.
        gene GetGene(int i) { return population[i];}
 
    protected:
        int ca,cb,cc,cd,ce;// The coefficients.
        int result;
        gene population[MAXPOP];// Population.
 
        int Fitness(gene &);// Fitness function.
        void GenerateLikelihoods(); // Generate likelihoods.
        float MultInv();// Creates the multiplicative inverse.
        int CreateFitnesses();
        void CreateNewPopulation();
        int GetIndex(float val);
 
        gene Breed(int p1, int p2);
 
};
void print(int* b)
{
    predmet A[100];
    ifstream ifst("input.txt");
    int n=0;
    for(;!ifst.eof();)
    {
        A[n].In(ifst);
        n++;
    }
    for(int i=0; i<n; i++) if(b[i]==1)
    A[i].prin();
}
CDiophantine::CDiophantine(int a, int b, int c, int d, int e, int res) : ca(a), cb(b), cc(c), cd(d), ce(e), result(res) {}
 
int CDiophantine::Solve() 
{
    int fitness = -1;
    // Генерация начальной популяции.
    srand((unsigned)time(NULL));
    for(int i=0;i<MAXPOP;i++)// Заполните населения с номерами между
    {
        for (int j=0;j<5;j++)// 0 and the result.
        {
            population[i].alleles[j] = rand() % 2;
            cout<<population[i].alleles[j];
        }
        cout<<endl;
    }
    if (fitness = CreateFitnesses()) 
    {
        return fitness;
    }
    int iterations = 0;// Вести учет итераций.
    while (fitness != 0 || iterations < 50)// Повторяйте, пока решение не найдено, или более 50 итераций.
    {
        GenerateLikelihoods();// Create the likelihoods.
        CreateNewPopulation();
        if (fitness = CreateFitnesses()) 
        {
            return fitness;
        }
        iterations++;
    }
    return -1;
}
 
int CDiophantine::Fitness(gene &gn)
{
    int total = ca * gn.alleles[0] + cb * gn.alleles[1] + cc * gn.alleles[2] + cd * gn.alleles[3]+ ce * gn.alleles[4];
 
    return gn.fitness = abs(total - result);
}
 
int CDiophantine::CreateFitnesses() 
{
    float avgfit = 0;
    int fitness = 0;
    for(int i=0;i<MAXPOP;i++) {
        fitness = Fitness(population[i]);
        avgfit += fitness;
        if (fitness == 0) {
            return i;
        }
    }
    return 0;
}
float CDiophantine::MultInv()
{
    float sum = 0;
 
    for(int i=0;i<MAXPOP;i++) {
        sum += 1/((float)population[i].fitness);
    }
 
    return sum;
}
 
void CDiophantine::GenerateLikelihoods() 
{
    float multinv = MultInv();
 
    float last = 0;
    for(int i=0;i<MAXPOP;i++) {
        population[i].likelihood = last = last + ((1/((float)population[i].fitness) / multinv) * 100);
    }
}
int CDiophantine::GetIndex(float val) {
    float last = 0;
    for(int i=0;i<MAXPOP;i++) {
        if (last <= val && val <= population[i].likelihood) return i;
        else last = population[i].likelihood;
    }
    
    return 5;
}
gene CDiophantine::Breed(int p1, int p2)
{
    int crossover = rand() % 4+1;// Создать точку пересечения (не первый).
    int first = rand() % 100;// Кто из родителей на первом месте?
 
    gene child = population[p1];// Child is all first parent initially.Ребенок все первый родитель на начальном этапе.
 
    int initial = 0, final = 4;// Кроссовер границ.
    if (first < 50) initial = crossover;    // If first parent first. start from crossover.  Если первый родитель в первую очередь. начать с кроссовером.
    else final = crossover+1;// Else end at crossover.
 
    for(int i=initial;i<final;i++)
    {// Crossover!
        child.alleles[i] = population[p2].alleles[i];
        if (rand() % 101 < 6) child.alleles[i] = rand() % (result + 1);
    }
 
    return child;// Return the kid...
}
 
void CDiophantine::CreateNewPopulation()
{
    gene temppop[MAXPOP];
 
    for(int i=0;i<MAXPOP;i++) 
    {
        int parent1 = 0, parent2 = 0, iterations = 0;
        while(parent1 == parent2 || population[parent1] == population[parent2])
        {
            parent1 = GetIndex((float)(rand() % 101));
            parent2 = GetIndex((float)(rand() % 101));
            if (++iterations > 25) break;
        }
    temppop[i] = Breed(parent1, parent2);// Create a child.
    }
 
    for(int i=0;i<MAXPOP;i++) population[i] = temppop[i];
}
void main() 
{
    setlocale(LC_ALL,"russian");
    ifstream ifst("input.txt");
    predmet A[100];
    int g,ans,P=0,n=0,m=0,W=0;
    int *a = new int[m]; //массив содержащий веса
    int *b = new int[n];
    for(;!ifst.eof();){A[m].In(ifst);m++;}
    for (int i=0;i<m;i++) a[i]=A[i].get_w();
    
    cout<<" Введите вес рюкзака:"<<endl;cin>>g;
    cout<<"\t** Полученные хромосомы:"<<endl;
   CDiophantine dp(a[0],a[1],a[2],a[3],a[4],g);
   ans = dp.Solve();
   if (ans == -1) 
   {
      cout << "No solution found." << endl;
   }
   else 
   {
      cout << "\t** Вес рюкзака равен "<<g<<" **"<<endl;
      cout << "\n Все предметы: \n"<<endl;
      ifstream ifst("input.txt");
      for(;!ifst.eof();)
      {
            A[n].In(ifst);
            A[n].prin();
            n++;
      }
      gene gn = dp.GetGene(ans);
      cout << "\n\t*** Наилучшая хромосома: ";
      cout << "( "<<gn.alleles[0] << " . ";
      cout << gn.alleles[1] << " . ";
      cout << gn.alleles[2] << " . ";
      cout << gn.alleles[3] << " . ";
      cout << gn.alleles[4] << " ) ***" << endl;
      for(int i=0;i<5;i++) b[i]=gn.alleles[i];
      cout<<"\n В рюкзаке присутствуют:\n\n";
      print(b);
      for (int i=0;i<n;i++) P+=A[i].get_p()*gn.alleles[i];
      cout<<"\n\t** Стоимость рюкзака равна "<<P<<" **"<<endl;
       cin.get();
       cin.get();
 
   }
}