Qurak

Numerical nonlinear global optimization examples queens on a chessboard

Підручник 6 кроків, що виконуються по порядку в одній сесії. Кожен крок було виконано на рушії, і вивід нижче — те, що він дав.

Крок 1
Clear[f, sol];

Немає виводу — цей крок готує щось для наступного.

Крок 2
attackQ[{i_Integer, xi_Integer}, {j_Integer, xj_Integer}] := Or[i - j == 0, xi - xj == 0, i - j == xi - xj, i - j == xj - xi]

Немає виводу — цей крок готує щось для наступного.

Крок 3
countAttacks[q : {_Real..}] := Block[{qns = Ordering[q], n = Length[q]}, Count[Flatten[Table[attackQ[{i, qns[[i]]}, {j, qns[[j]]}], {i, 1, n - 1}, {j, i + 1, n}]], True]]

Немає виводу — цей крок готує щось для наступного.

Крок 4
queenPlot[vec_]  := Module[{veclength = Length[vec]}, Show[Graphics[Join[ Table[Line[{{1, i}, {veclength + 1, i}}], {i, veclength + 1}], Table[Line[{{i, 1}, {i, veclength + 1}}], {i, veclength + 1}]]], Graphics[Table[Disk[{i + .5, vec[[i]]  + .5}, 0.25], {i, veclength}]], AspectRatio -> 1]]

Немає виводу — цей крок готує щось для наступного.

Крок 5
vars  = Array[Subscript[q, #]&, 8];
{f, sol} = NMinimize[countAttacks[vars], vars, Method -> {"DifferentialEvolution", "PostProcess" -> False}];
{f, Thread[vars -> (sol = Ordering[sol[[All, 2]]])]}
Вивід
{f, {Subscript[q, 1] -> Ordering[sol[[2]]], Subscript[q, 2] -> Ordering[sol[[2]]], Subscript[q, 3] -> Ordering[sol[[2]]], Subscript[q, 4] -> Ordering[sol[[2]]], Subscript[q, 5] -> Ordering[sol[[2]]], Subscript[q, 6] -> Ordering[sol[[2]]], Subscript[q, 7] -> Ordering[sol[[2]]], Subscript[q, 8] -> Ordering[sol[[2]]]}}
Крок 6
queenPlot[sol]
Вивід
-Graphics-

Використані функції

Пов'язані рецепти

Усі рецепти · Довідник функцій · Використати з MCP-клієнта