defmodule CPSolver.Examples.Zebra do @moduledoc """ https://en.wikipedia.org/wiki/Zebra_Puzzle There are five houses. The Englishman lives in the red house. The Spaniard owns the dog. Coffee is drunk in the green house. The Ukrainian drinks tea. The green house is immediately to the right of the ivory house. The Old Gold smoker owns snails. Kools are smoked in the yellow house. Milk is drunk in the middle house. The Norwegian lives in the first house. The man who smokes Chesterfields lives in the house next to the man with the fox. Kools are smoked in the house next to the house where the horse is kept. The Lucky Strike smoker drinks orange juice. The Japanese smokes Parliaments. The Norwegian lives next to the blue house. Now, who drinks water? Who owns the zebra? """ @doc """ The model is a clone of MiniZinc model (https://github.com/hakank/hakank/blob/39a2b9e868011df38bd7c90f530f8c5b3e0740cb/minizinc/zebra_inverse.mzn) by HÃ¥kan Kjellerstrand """ alias CPSolver.IntVariable, as: Variable alias CPSolver.Model import CPSolver.Constraint.Factory import CPSolver.Variable.View.Factory def model() do ### Variables domain = 0..4 color_vars = [red, green, ivory, yellow, blue] = create_vars([:red, :green, :ivory, :yellow, :blue], domain) nationality_vars = [englishman, spaniard, ukrainian, norwegian, japanese] = create_vars([:englishman, :spaniard, :ukrainian, :norwegian, :japanese], domain) animal_vars = [dog, snail, fox, horse, _zebra] = create_vars([:dog, :snail, :fox, :horse, :zebra], domain) drink_vars = [coffee, tea, milk, orange_juice, _water] = create_vars([:coffee, :tea, :milk, :orange_juice, :water], domain) brand_vars = [old_gold, kool, chesterfield, lucky_strike, parliament] = create_vars([:old_gold, :kool, :chesterfield, :lucky_strike, :parliament], domain) ### Constraints inverse_constraints = Enum.map( [color_vars, nationality_vars, animal_vars, drink_vars, brand_vars], fn vars -> inverse(create_tmp_vars(domain), vars) end ) constraints = [ equal(englishman, red), equal(spaniard, dog), equal(coffee, green), equal(ukrainian, tea), equal(green, inc(ivory, 1)), equal(old_gold, snail), equal(kool, yellow), equal(milk, house(3)), equal(norwegian, house(1)), next_to(chesterfield, fox), next_to(kool, horse), equal(lucky_strike, orange_juice), equal(japanese, parliament), next_to(norwegian, blue), inverse_constraints ] Model.new( color_vars ++ nationality_vars ++ animal_vars ++ drink_vars ++ brand_vars, constraints ) end def solve(opts \\ []) do CPSolver.solve(model(), opts) end defp create_tmp_vars(domain) do Enum.map(1..5, fn _ -> Variable.new(domain) end) end defp create_vars(names, domain) do Enum.map(names, fn name -> Variable.new(domain, name: name) end) end defp next_to(var1, var2) do {tmp, c} = subtract(var1, var2) [c, absolute(tmp, 1)] end defp find_match(block1, block2, position_to_match) do value_to_match = Enum.at(block2, position_to_match) Enum.find_index(block1, fn x -> x == value_to_match end) end def puzzle_solution(res) do solution = hd(res.solutions) nationality_names = Enum.slice(res.variables, 5, 5) [_colors, nationalities, animals, drinks, _brands] = Enum.chunk_every(solution, 5) |> Enum.take(5) zebra_owner = find_match(nationalities, animals, 4) water_drinker = find_match(nationalities, drinks, 4) %{ zebra_owner: Enum.at(nationality_names, zebra_owner), water_drinker: Enum.at(nationality_names, water_drinker) } end ## Shift to 0-based defp house(number) do number - 1 end end