Xij_Pyomo.ipynb

{
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/OptimizationExpert/Pyomo/blob/main/Xij_Pyomo.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 19,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "XUECKn64jg4v",
        "outputId": "09d28f8c-159f-4283-8a40-18f6e05c5aa7"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Requirement already satisfied: pyomo in /usr/local/lib/python3.10/dist-packages (6.7.0)\n",
            "Requirement already satisfied: ply in /usr/local/lib/python3.10/dist-packages (from pyomo) (3.11)\n",
            "Model unknown\n",
            "\n",
            "  Variables:\n",
            "    X : Size=100, Index=X_index\n",
            "        Key      : Lower : Value : Upper : Fixed : Stale : Domain\n",
            "          (1, 1) :     0 :  None :     1 : False :  True :  Reals\n",
            "          (1, 2) :     0 :   1.0 :     1 : False : False :  Reals\n",
            "          (1, 3) :     0 :   1.0 :     1 : False : False :  Reals\n",
            "          (1, 4) :     0 :   1.0 :     1 : False : False :  Reals\n",
            "          (1, 5) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (1, 6) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (1, 7) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (1, 8) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (1, 9) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (1, 10) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (2, 1) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (2, 2) :     0 :  None :     1 : False :  True :  Reals\n",
            "          (2, 3) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (2, 4) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (2, 5) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (2, 6) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (2, 7) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (2, 8) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (2, 9) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (2, 10) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (3, 1) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (3, 2) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (3, 3) :     0 :  None :     1 : False :  True :  Reals\n",
            "          (3, 4) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (3, 5) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (3, 6) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (3, 7) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (3, 8) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (3, 9) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (3, 10) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (4, 1) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (4, 2) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (4, 3) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (4, 4) :     0 :  None :     1 : False :  True :  Reals\n",
            "          (4, 5) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (4, 6) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (4, 7) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (4, 8) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (4, 9) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (4, 10) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (5, 1) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (5, 2) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (5, 3) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (5, 4) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (5, 5) :     0 :  None :     1 : False :  True :  Reals\n",
            "          (5, 6) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (5, 7) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (5, 8) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (5, 9) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (5, 10) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (6, 1) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (6, 2) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (6, 3) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (6, 4) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (6, 5) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (6, 6) :     0 :  None :     1 : False :  True :  Reals\n",
            "          (6, 7) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (6, 8) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (6, 9) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (6, 10) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (7, 1) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (7, 2) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (7, 3) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (7, 4) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (7, 5) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (7, 6) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (7, 7) :     0 :  None :     1 : False :  True :  Reals\n",
            "          (7, 8) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (7, 9) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (7, 10) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (8, 1) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (8, 2) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (8, 3) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (8, 4) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (8, 5) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (8, 6) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (8, 7) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (8, 8) :     0 :  None :     1 : False :  True :  Reals\n",
            "          (8, 9) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (8, 10) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (9, 1) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (9, 2) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (9, 3) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (9, 4) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (9, 5) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (9, 6) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (9, 7) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (9, 8) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "          (9, 9) :     0 :  None :     1 : False :  True :  Reals\n",
            "         (9, 10) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (10, 1) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (10, 2) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (10, 3) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (10, 4) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (10, 5) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (10, 6) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (10, 7) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (10, 8) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "         (10, 9) :     0 :   0.0 :     1 : False : False :  Reals\n",
            "        (10, 10) :     0 :  None :     1 : False :  True :  Reals\n",
            "\n",
            "  Objectives:\n",
            "    obj1 : Size=1, Index=None, Active=True\n",
            "        Key  : Active : Value\n",
            "        None :   True :   3.0\n",
            "\n",
            "  Constraints:\n",
            "    C1 : Size=1\n",
            "        Key  : Lower : Body : Upper\n",
            "        None :   3.0 :  3.0 :   3.0\n"
          ]
        }
      ],
      "source": [
        "!pip install pyomo\n",
        "from pyomo.environ import *\n",
        "!apt-get install -y -qq glpk-utils\n",
        "opt = SolverFactory('glpk', executable='/usr/bin/glpsol')\n",
        "\n",
        "N = 10\n",
        "model = AbstractModel()\n",
        "model.i = RangeSet(N)\n",
        "model.X = Var(model.i,model.i, bounds=(0,1), within=Reals)\n",
        "\n",
        "def rule_C1(model):\n",
        "  return sum (model.X[i,j] for i in model.i for j in model.i if i!=j) == 3\n",
        "model.C1=Constraint(rule=rule_C1)\n",
        "\n",
        "\n",
        "def rule_OF(model):\n",
        "    return sum (model.X[i,j] for i in model.i for j in model.i if i!=j)\n",
        "model.obj1 = Objective(rule=rule_OF, sense=minimize)\n",
        "instance = model.create_instance()\n",
        "results = opt.solve(instance)\n",
        "instance.display()"
      ]
    }
  ],
  "metadata": {
    "kernelspec": {
      "display_name": "Python 3",
      "language": "python",
      "name": "python3"
    },
    "language_info": {
      "codemirror_mode": {
        "name": "ipython",
        "version": 3
      },
      "file_extension": ".py",
      "mimetype": "text/x-python",
      "name": "python",
      "nbconvert_exporter": "python",
      "pygments_lexer": "ipython3",
      "version": "3.10.7"
    },
    "colab": {
      "provenance": [],
      "include_colab_link": true
    }
  },
  "nbformat": 4,
  "nbformat_minor": 0
}