Pi Hybrids问题
清华大学肖秀波梁湧老师翻译的Rardin教授的《运筹学》[1]已于今年年中出版，感谢机械工业出版社张有利老师的推荐和赠书，让我能看到如此完美的千页级宏篇译著。该书的翻译质量非常高，书中内容深入浅出，附有大量的应用案例（Application）和练习题库。尤其让人欣喜的是该著作能与计算实践密切结合，凡有计算机图案标记的练习都是与计算和软件应用相关的，彰显了运筹学的应用数学本质。有鉴于此，笔者计划对该书中的问题进行大量实践，争取达到双位数的规模。
Pi Hybrids问题
Pi Hybrids问题，见原书第二章，案例2-2(Pi Hybrids, Page 33,37,56-58)。
Pi Hybrids问题是一个叠加了生产问题的运输问题，是一个包含了三脚标标引变量的模型，该著作在一开始就介入如此有一定深度的模型展现了作者应用导向的理念。
+Leapms生成的模型摘录
对任何问题，+Leapms都建议直接用+Leapms建模语言直接写出模型，并进行模型调试。
当模型调试完毕，对模型的正确性有充分的信心后，+Leapms系统可生成模型摘录，包括数学概念模型（即使用标引符号表示的数学模型）和+Leapms源码供交流使用。
我们的经验是，模型调试能够发现模型的不足、促进建模的完美性，甚至可以促进对问题的更加深入的了解。
以下是问题的+Leapms模型摘录（pdf屏幕截图）：
+Leapms模型求解结果
+Leapms>load
Current directory is "ROOT".
.........
RLX2-.leap
.........
please input the filename:RLX2-
================================================================
: //Rardin R. L 著，肖勇波、梁湧译. 运筹学. 北京:机械工业出版社,2018
: //Rardin-Liang&Xiao:案例2-2(Pi Hybrids, Page 33,37,56-58)
:
: min sum{f=,..,l;h=,..,m}p[f][h]x[f][h]+ -->
: sum{f=,..,l;h=,..,m;r=,..,n}s[f][h][r]y[f][h][r]
: s.t.
: sum{h=,..,m}a[h]x[f][h]<=u[f] | f=,..,l
: sum{f=,..,l}y[f][h][r]=d[h][r]| h=,..,m;r=,..,n
: sum{r=,..,n}y[f][h][r]=x[f][h]| f=,..,l;h=,..,m
: where
: l,m,n are integers
: a[h] is a number|h=,..,m
: u[f] is a number|f=,..,l
: p[f][h] is a number|f=,..,l;h=,..,m
: d[h][r] is a number|h=,...,m;r=,...,n
: s[f][h][r] is a number|f=,..,l;h=,..,m;r=,..,n
: x[f][h] is a variable of nonnegative number-->
: |f=,..,l;h=,..,m
: y[f][h][r] is a variable of nonnegative number-->
: |f=,..,l;h=,..,m;r=,..,n
: data
: l=
: m=
: n=
: u={ } // corn processing capacity of facility f
: a={ } // number of corn to obtain a hybrid h
: p={
: 1.10 0.89 2.05 1.45
: 1.55 1.13 2.15 1.56
: 0.95 0.83 1.80 1.22
: }
: d={ //number of hybrid h demanded in region r
:
:
:
:
: }
: s={ //cost of shipping hybrid h from facility f to region r
: 0.89 0.91 0.77
: 1.00 0.84 0.89
: 0.77 0.76 0.78
: 0.99 1.03 0.85
:
: 0.92 0.89 0.92
: 0.87 0.95 0.90
: 0.91 0.83 0.77
: 0.89 0.79 0.86
: }
:
================================================================
>>end of the file.
Parsing model:
1D
2R
3V
4O
5C
6S
7End.
..................................
number of variables=
number of constraints=
..................................
+Leapms>solve
The LP is solved to optimal.
找到线性规划最优解.非零变量值和最优目标值如下：
.........
x1_1*=
x1_2*=
x1_3*=
x1_4*=93.7222
x2_4*=736.278
y1_1_1*=
y1_1_2*=
y1_1_3*=
y1_2_1*=
y1_2_2*=
y1_2_3*=
y1_3_1*=
y1_3_2*=
y1_3_3*=
y1_4_3*=93.7222
y2_4_1*=
y2_4_2*=
y2_4_3*=346.278
.........
Objective*=7400.61
.........
+Leapms>
讨论
在原书中，Pi Hybrids问题是讲述标引模型（即对大规模模型中的变量和常量符号等实用脚标标引形式给出）的举例。+Leapms是一种以标引模型为主要表达方式的语言。
另外，在英文原版书中，参数数据可能有误，导致生产能力不足以满足市场需求，此处把生产能力参数 u 从英文原书中的{2200 2555}更改为{22000 25550}。
参考文献
[1] Rardin R. L 著，肖勇波、梁湧译. 运筹学. 北京:机械工业出版社,2018