(暴躁)

昨晚一晚上在跑，早上起来发现系统崩溃了。
Maximize[{(Subscript[p,11] + Subscript[p,12] + Subscript[p,13])*3100 + (Subscript[p,21] + Subscript[p,22] + Subscript[p,23])*3800 + (Subscript[p,31] + Subscript[p,32] + Subscript[p,33])*3500 + (Subscript[p,41] + Subscript[p,42] + Subscript[p, 43])*2850,Subscript[p,11] + Subscript[p,21] + Subscript[p,31] + Subscript[p,41] <= 10 && Subscript[p,12] + Subscript[p,22] + Subscript[p,32] + Subscript[p,42] <= 16 && Subscript[p,13] + Subscript[p,23] + Subscript[p,33] + Subscript[p,43] <= 8 && Subscript[p,11]*480 + Subscript[p,21]*650 + Subscript[p,31]*580 + Subscript[p, 41]*390 <= 6800 && Subscript[p,12]*480 + Subscript[p,22]*650 + Subscript[p,32]*580 + Subscript[p,42]*390 <= 8700 && Subscript[p,13]*480 + Subscript[p,23]*650 + Subscript[p,33]*580 + Subscript[p,43]*390 <= 5300 && Subscript[p,11] + Subscript[p,12] + Subscript[p,13] <= 18 && Subscript[p,21] + Subscript[p,22] + Subscript[p,23] <= 15 && Subscript[p,31] + Subscript[p,32] + Subscript[p,33] <= 23 && Subscript[p,41] + Subscript[p,42] + Subscript[p,43] <= 12 && (Subscript[p,11] + Subscript[p,21] + Subscript[p,31] + Subscript[p,41])/(Subscript[p,12] + Subscript[p,22] + Subscript[p,32] + Subscript[p,42]) == 10/16 && (Subscript[p,12] + Subscript[p,22] + Subscript[p,32] + Subscript[p,42])/(Subscript[p,13] + Subscript[p, 23] + Subscript[p,33] + Subscript[p,43]) == 16/8},{Subscript[p,11], Subscript[p,12], Subscript[p,13], Subscript[p,21], Subscript[p,22], Subscript[p,23],Subscript[p,31], Subscript[p,32], Subscript[p,33], Subscript[p,41], Subscript[p,42], Subscript[p,43]}]

很奇怪，为什么用这么多Subscript，直接p11不久行了。。。

对于这种情况，个人建议优化模型……

这貌似不是mma的专长

我用Lingo做了


mathematica结果：{2308800/19, {Subscript[p, 11] -> 0, Subscript[p, 12] -> 0, Subscript[p, 13] -> 0, Subscript[p, 21] -> 7, Subscript[p, 22] -> 0, Subscript[p, 23] -> 8, Subscript[p, 31] -> 3, Subscript[p, 32] -> 246/19, Subscript[p, 33] -> 0, Subscript[p, 41] -> 0, Subscript[p, 42] -> 58/19, Subscript[p, 43] -> 0}}

(暴躁)