Ejemplo n.º 1
0
        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);
        }