public static CombinatorialSudoku New(CombinatorialSudoku sudoku)
{
return(new CombinatorialSudoku(MatrixHelper.DuplicateMatrix(sudoku.CellValues)));
}
private CombinatorialSudoku SolveInternal(CombinatorialSudoku sudoku, int numOrganisms, int maxEpochs)
{
var numberOfWorkers = (int)(numOrganisms * 0.90);
var hive = new Organism[numOrganisms];
var bestError = int.MaxValue;
CombinatorialSudoku bestSolution = null;
for (var i = 0; i < numOrganisms; ++i)
{
var organismType = i < numberOfWorkers
? OrganismType.Worker
: OrganismType.Explorer;
var randomSudoku = CombinatorialSudoku.New(MatrixHelper.RandomMatrix(_rnd, sudoku.CellValues));
var err = randomSudoku.Error;
hive[i] = new Organism(organismType, randomSudoku.CellValues, err, 0);
if (err >= bestError)
{
continue;
}
bestError = err;
bestSolution = CombinatorialSudoku.New(randomSudoku);
}
var epoch = 0;
while (epoch < maxEpochs)
{
if (epoch % 1000 == 0)
{
Console.WriteLine($"Epoch: {epoch}, Best error: {bestError}");
}
if (bestError == 0)
{
break;
}
for (var i = 0; i < numOrganisms; ++i)
{
if (hive[i].Type == OrganismType.Worker)
{
var neighbor = MatrixHelper.NeighborMatrix(_rnd, sudoku.CellValues, hive[i].Matrix);
var neighborSudoku = CombinatorialSudoku.New(neighbor);
var neighborError = neighborSudoku.Error;
var p = _rnd.NextDouble();
if (neighborError < hive[i].Error || p < 0.001)
{
hive[i].Matrix = MatrixHelper.DuplicateMatrix(neighbor);
hive[i].Error = neighborError;
if (neighborError < hive[i].Error)
{
hive[i].Age = 0;
}
if (neighborError >= bestError)
{
continue;
}
bestError = neighborError;
bestSolution = neighborSudoku;
}
else
{
hive[i].Age++;
if (hive[i].Age <= 1000)
{
continue;
}
var randomSudoku = CombinatorialSudoku.New(MatrixHelper.RandomMatrix(_rnd, sudoku.CellValues));
hive[i] = new Organism(0, randomSudoku.CellValues, randomSudoku.Error, 0);
}
}
else
{
var randomSudoku = CombinatorialSudoku.New(MatrixHelper.RandomMatrix(_rnd, sudoku.CellValues));
hive[i].Matrix = MatrixHelper.DuplicateMatrix(randomSudoku.CellValues);
hive[i].Error = randomSudoku.Error;
if (hive[i].Error >= bestError)
{
continue;
}
bestError = hive[i].Error;
bestSolution = randomSudoku;
}
}
// merge best worker with best explorer into worst worker
var bestWorkerIndex = 0;
var smallestWorkerError = hive[0].Error;
for (var i = 0; i < numberOfWorkers; ++i)
{
if (hive[i].Error >= smallestWorkerError)
{
continue;
}
smallestWorkerError = hive[i].Error;
bestWorkerIndex = i;
}
var bestExplorerIndex = numberOfWorkers;
var smallestExplorerError = hive[numberOfWorkers].Error;
for (var i = numberOfWorkers; i < numOrganisms; ++i)
{
if (hive[i].Error >= smallestExplorerError)
{
continue;
}
smallestExplorerError = hive[i].Error;
bestExplorerIndex = i;
}
var worstWorkerIndex = 0;
var largestWorkerError = hive[0].Error;
for (var i = 0; i < numberOfWorkers; ++i)
{
if (hive[i].Error <= largestWorkerError)
{
continue;
}
largestWorkerError = hive[i].Error;
worstWorkerIndex = i;
}
var merged = MatrixHelper.MergeMatrices(_rnd, hive[bestWorkerIndex].Matrix, hive[bestExplorerIndex].Matrix);
var mergedSudoku = CombinatorialSudoku.New(merged);
hive[worstWorkerIndex] = new Organism(0, merged, mergedSudoku.Error, 0);
if (hive[worstWorkerIndex].Error < bestError)
{
bestError = hive[worstWorkerIndex].Error;
bestSolution = mergedSudoku;
}
++epoch;
}
return(bestSolution);
}
public static CombinatorialSudoku New(int[,] cellValues)
{
return(new CombinatorialSudoku(MatrixHelper.DuplicateMatrix(cellValues)));
}