#!/usr/bin/env python3 # Copyright 2010-2025 Google LLC # Licensed under the Apache License, Version 2.0 (the "License"); # you may not use this file except in compliance with the License. # You may obtain a copy of the License at # # http://www.apache.org/licenses/LICENSE-2.0 # # Unless required by applicable law or agreed to in writing, software # distributed under the License is distributed on an "AS IS" BASIS, # WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. # See the License for the specific language governing permissions and # limitations under the License. """Test linear sum assignment on a 4x4 matrix. Example taken from: http://www.ee.oulu.fi/~mpa/matreng/eem1_2-1.htm with kCost[0][1] modified so the optimum solution is unique. """ from typing import Sequence from absl import app from ortools.graph.python import linear_sum_assignment def run_assignment_on_4x4_matrix(): """Test linear sum assignment on a 4x4 matrix.""" num_sources = 4 num_targets = 4 cost = [ [90, 76, 75, 80], [35, 85, 55, 65], [125, 95, 90, 105], [45, 110, 95, 115], ] expected_cost = cost[0][3] + cost[1][2] + cost[2][1] + cost[3][0] assignment = linear_sum_assignment.SimpleLinearSumAssignment() for source in range(0, num_sources): for target in range(0, num_targets): assignment.add_arc_with_cost(source, target, cost[source][target]) solve_status = assignment.solve() if solve_status == assignment.OPTIMAL: print("Successful solve.") print("Total cost", assignment.optimal_cost(), "/", expected_cost) for i in range(0, assignment.num_nodes()): print( "Left node %d assigned to right node %d with cost %d." % (i, assignment.right_mate(i), assignment.assignment_cost(i)) ) elif solve_status == assignment.INFEASIBLE: print("No perfect matching exists.") elif solve_status == assignment.POSSIBLE_OVERFLOW: print("Some input costs are too large and may cause an integer overflow.") def main(argv: Sequence[str]) -> None: if len(argv) > 1: raise app.UsageError("Too many command-line arguments.") run_assignment_on_4x4_matrix() if __name__ == "__main__": app.run(main)