public Schedule Optimize()
{
int[] taskIDs = tasks.Keys.ToArray();
Genotype[] population = new Genotype[1000];
for (int i = 0; i < 1000; i++)
{
int[] order = taskIDs.OrderBy(x => rnd.Next()).ToArray();
schedule.ChangeOrderAndUpdate(order);
population[i] = new Genotype(order, schedule.GetTotalCost());
}
population = population.OrderBy(x => x.Comp).ToArray();
double SUM = (1000 * (1 + 1000)) / 2; // sum from 1 to 1000
for (int p = 0; p < generations; p++)
{
// distrib
population[0].Beg = 0.0;
population[0].End = 1.0 / SUM;
for (int i = 1; i < 1000; i++)
{
population[i].Beg = population[i - 1].End;
population[i].End = ((double)(i + 1)) / SUM + population[i].Beg;
}
Genotype[] children = new Genotype[2000];
int insert = 0;
for (int w = 0; w < 1000; w++)
{
// get parents
Genotype parent1 = null, parent2 = null;
double x, y;
do
{
x = rnd.NextDouble();
y = rnd.NextDouble();
for (int i = 0; i < 1000; i++)
{
double beg = population[i].Beg;
double end = population[i].End;
if (beg <= x && x <= end)
{
parent1 = population[i];
}
if (beg <= y && y <= end)
{
parent2 = population[i];
}
}
} while (parent1 == parent2);
Genotype child1, child2;
if (rnd.NextDouble() > pk)
{
// make two children from conbination of parent genes
if (crossOX)
{
child1 = CrossOX(parent1, parent2);
child2 = CrossOX(parent2, parent1);
}
else
{
child1 = CrossPoint(parent1, parent2);
child2 = CrossPoint(parent2, parent1);
}
}
else
{
child1 = new Genotype(parent1.Genes, parent1.Cost);
child2 = new Genotype(parent2.Genes, parent2.Cost);
}
if (rnd.NextDouble() < pm)
{
// swap random genes
child1 = SwapRandomGenes(child1);
child2 = SwapRandomGenes(child2);
}
children[insert++] = child1;
children[insert++] = child2;
}
children = children.OrderBy(x => x.Comp).ToArray();
for (int i = 1000; i < 2000; i++)
{
population[i - 1000] = children[i];
}
}
Genotype minGen = null;
int minCost = int.MaxValue;
for (int i = 0; i < 1000; i++)
{
if (population[i].Cost < minCost)
{
minGen = population[i];
minCost = minGen.Cost;
}
}
schedule.ChangeOrderAndUpdate(minGen.Genes);
return(schedule);
}