csintegerprogramming.cs 4.65 KB
Newer Older
Valentin Platzgummer's avatar
Valentin Platzgummer committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108
// Copyright 2010-2018 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.

using System;
using Google.OrTools.LinearSolver;

public class CsIntegerProgramming
{
  private static void RunIntegerProgrammingExample(String solverType)
  {
    Solver solver = Solver.CreateSolver("IntegerProgramming", solverType);
    if (solver == null)
    {
      Console.WriteLine("Could not create solver " + solverType);
      return;
    }
    // x1 and x2 are integer non-negative variables.
    Variable x1 = solver.MakeIntVar(0.0, double.PositiveInfinity, "x1");
    Variable x2 = solver.MakeIntVar(0.0, double.PositiveInfinity, "x2");

    // Minimize x1 + 2 * x2.
    Objective objective = solver.Objective();
    objective.SetMinimization();
    objective.SetCoefficient(x1, 1);
    objective.SetCoefficient(x2, 2);

    // 2 * x2 + 3 * x1 >= 17.
    Constraint ct = solver.MakeConstraint(17, double.PositiveInfinity);
    ct.SetCoefficient(x1, 3);
    ct.SetCoefficient(x2, 2);

    Solver.ResultStatus resultStatus = solver.Solve();

    // Check that the problem has an optimal solution.
    if (resultStatus != Solver.ResultStatus.OPTIMAL)
    {
      Console.WriteLine("The problem does not have an optimal solution!");
      return;
    }

    Console.WriteLine("Problem solved in " + solver.WallTime() +
                      " milliseconds");

    // The objective value of the solution.
    Console.WriteLine("Optimal objective value = " + objective.Value());

    // The value of each variable in the solution.
    Console.WriteLine("x1 = " + x1.SolutionValue());
    Console.WriteLine("x2 = " + x2.SolutionValue());

    Console.WriteLine("Advanced usage:");
    Console.WriteLine("Problem solved in " + solver.Nodes() +
                      " branch-and-bound nodes");
  }

  private static void RunIntegerProgrammingExampleNaturalApi(String solverType)
  {
    Solver solver = Solver.CreateSolver("IntegerProgramming", solverType);
    if (solver == null)
    {
      Console.WriteLine("Could not create solver " + solverType);
      return;
    }
    // x1 and x2 are integer non-negative variables.
    Variable x1 = solver.MakeIntVar(0.0, double.PositiveInfinity, "x1");
    Variable x2 = solver.MakeIntVar(0.0, double.PositiveInfinity, "x2");

    solver.Minimize(x1 + 2 * x2);
    solver.Add(2 * x2 + 3 * x1 >= 17);

    Solver.ResultStatus resultStatus = solver.Solve();

    // Check that the problem has an optimal solution.
    if (resultStatus != Solver.ResultStatus.OPTIMAL)
    {
      Console.WriteLine("The problem does not have an optimal solution!");
      return;
    }

    Console.WriteLine("Problem solved in " + solver.WallTime() +
                      " milliseconds");

    // The objective value of the solution.
    Console.WriteLine("Optimal objective value = " +
                      solver.Objective().Value());

    // The value of each variable in the solution.
    Console.WriteLine("x1 = " + x1.SolutionValue());
    Console.WriteLine("x2 = " + x2.SolutionValue());

    Console.WriteLine("Advanced usage:");
    Console.WriteLine("Problem solved in " + solver.Nodes() +
                      " branch-and-bound nodes");
  }

  static void Main()
  {
    Console.WriteLine("---- Integer programming example with GLPK ----");
109
    RunIntegerProgrammingExample("GLPK");
Valentin Platzgummer's avatar
Valentin Platzgummer committed
110
    Console.WriteLine("---- Linear programming example with CBC ----");
111
    RunIntegerProgrammingExample("CBC");
Valentin Platzgummer's avatar
Valentin Platzgummer committed
112
    Console.WriteLine("---- Linear programming example with SCIP ----");
113 114 115
    RunIntegerProgrammingExample("SCIP");
    Console.WriteLine("---- Linear programming example with SAT ----");
    RunIntegerProgrammingExample("SAT");
Valentin Platzgummer's avatar
Valentin Platzgummer committed
116 117
    Console.WriteLine(
        "---- Integer programming example (Natural API) with GLPK ----");
118
    RunIntegerProgrammingExampleNaturalApi("GLPK");
Valentin Platzgummer's avatar
Valentin Platzgummer committed
119 120
    Console.WriteLine(
        "---- Linear programming example (Natural API) with CBC ----");
121
    RunIntegerProgrammingExampleNaturalApi("CBC");
Valentin Platzgummer's avatar
Valentin Platzgummer committed
122 123
    Console.WriteLine(
        "---- Linear programming example (Natural API) with SCIP ----");
124 125 126 127
    RunIntegerProgrammingExampleNaturalApi("SCIP");
    Console.WriteLine(
        "---- Linear programming example (Natural API) with SAT ----");
    RunIntegerProgrammingExampleNaturalApi("SAT");
Valentin Platzgummer's avatar
Valentin Platzgummer committed
128 129
  }
}