用更复杂的模型做简单的实验

现在让我们考虑一个更困难的模式，glass4.mps．我们再次读取模型并开始优化:

Gurobi > m = 'c: / gurobi811 / win64/ /数据/例子glass4”) 从文件中读取MPS格式模型c: / gurobi811 / win64/ /数据/ glass4.mps例子

Reading time = 0.00 seconds glass4: 396 Rows, 322 Columns, 1815 NonZeros gurobi> m.optimize() Optimize a model with 396 Rows, 322 Columns and 1815 NonZeros矩阵范围[1e+00, 8e+06]目标范围[1e+00, 1e+06]边界范围[1e+00, 8e+02] RHS范围[1e+00, 8e+06]发现启发式解:目标2.40002e+09预解删除4行5列预解时间:0.01s预解:392行，317列，1815非零变量类型:19连续，298整数(298二进制)根松弛:目标8.000024e+08, 72次迭代，0.00秒节点当前节点| |客观界限|工作Expl Unexpl | Obj深度IntInf | |现任BestBd差距/节点时间0 0 8.0000 2.4000 e + e + 08年0 72 09年8.0000 e + 08年66.7% - 2.283353 0 0 0 H e + 09年8.0000 e + 08年65.0% - 0 H 0 0 2.283353 e + 09年8.0000 e + 08年65.0% - 0 0 0 8.0000 2.2834 e + e + 08年0 72 09年8.0000 e + 08年65.0% - 0 0 0 8.0000 2.2834 e + e + 08年0 72 09年8.0000 e + 08年65.0% - 0 0 0 8.0000 2.2834 e + e + 08年0 72 09年8.0000 e + 08年65.0% - 0 H 0 0 2.233352 e + 09年8.0000 e + 08年64.2% 2 - 0 0 8.0000 2.2334 e + e + 08年0 72 09年8.0000 e + 08年64.2% - 0 H 702 664 2.133352 e + 09年8.0000 e + 08年62.5% 5.1 0 * 779 663 81 2.000020 e + 09年8.0000 e + 08年60.0% 5.2 0 H 1396 1055 2.000019 e + 09年8.0000 e + 08年60.0% 5.4 0 H 1397 1011 2.000017 e + e + 08年09年8.000060.0% 5.4 0s * 1563 1071 46 2.000017e+09 8.0000e+08 60.0% 5.2 0s * 1566 1028 47 2.000017e+09 8.0000e+08 60.0% 5.2 0s * 1641 1034 46 2.000017e+09 8.0000e+08 60.0% 5.1 0s H 2751 1465 1.950016e+09 8.0000e+08 59.0% 4.5 1s * 3452 1815 129 1.946685e+09 8.0000e+08 58.9% 4.3 1s * 3628 1890 130 1.940020e+09 8.0000e+08 58.8% 4.3 1s * 3774 1952 103 1.912519e+09 8.0000e+08 58.2% 4.2 1s * 3778 1935 104 1.900018e+09 8.0000e+08 57.9% 4.2 1s * 5459 3246 111 1.900017e+09 8.0000e+08 57.9% 3.7 1s * 5706 3206 114 1.850017e+09 8.0000e+08 56.8% 3.7 1s * 5707 3206 114 1.850017e+09 8.0000e+08 56.8% 3.7 1s H 9808 6160 1.850017e+09 8.0000e+08 56.8% 3.3 1s * 9962 6053 75 1.800017e+09 8.0000e+08 55.6% 3.3 1s * 9963 6050 75 1.800017e+09 8.0000e+08 55.6% 3.3 1s *11813 6857 83 1.766683e+09 8.0000e+08 54.7% 3.2 2s H13611 5753 1.593348e+09 8.0000e+08 49.8% 3.2 2s H17448 7869 1.571443e+09 8.0000e+08 49.1% 3.1 2s H17455 7722 1.560015e+09 8.0000e+08 48.7% 3.1 2s H24494 11296 1.550014e+09 8.0000e+08 48.4% 3.4 3s H29764 13948 1.550014e+09 8.0001e+08 48.4% 3.4 7s 29791 13966 1.1000e+09 33 148 1.5500e+09 8.0431e+08 48.1% 3.4 10s H29794 13269 1.537514e+09 8.0448e+08 47.7% 3.4 10s 29861 13318 1.3000e+09 63 131 1.5375e+09 8.1068e+08 47.3% 3.5 15s Interrupt request received Cutting planes: Gomory: 42 Cover: 2 Implied bound: 10 Clique: 1 MIR: 6 Flow cover: 42 Explored 29862 nodes (105763 simplex iterations) in 15.38 seconds Thread count was 8 (of 8 available processors) Solution count 10: 1.53751e+09 1.55001e+09 1.55001e+09 ... 1.85002e+09 Solve interrupted Best objective 1.537514308333e+09, best bound 8.106816259965e+08, gap 47.2732%

很明显，这个模型比之前的模型要困难一些硬币模型。最优解实际上是1200000000年，但找到解决办法需要一段时间。让模型运行10秒钟后，我们中断运行(通过按CTRL-C，这将产生中断请求收到消息)，并考虑我们的选择。打字m.optimize ()会从被打断的地方重新开始。