我是新使用CPLEX，我试图找到互联网上的一些信息，但没有找到明确的东西，以帮助我在我的问题。
我有P [k]的k将等于1到4
我有一个决定变量x [I] [K]必须等于0或1（也P [k]的）
的i是介于1至5
现在我做这样的
IloEnv env;
IloModel model(env);
IloNumVarArray p(env);
p.add(IloNumVar(env, 0, 1));
p.add(IloNumVar(env, 0, 1));
p.add(IloNumVar(env, 0, 1));
IloIntVar x(env, 0, 1);
model.add(IloMaximize(env, 1000 * p[1] + 2000 * p[2] + 500 * p[3] + 1500 * p[4]));
for(int k = 1; k <= 4; k++){
for(int i = 1; i <= 5; i++){
model.add(x[i][k] + x[i][k] + x[i][k] + x[i][k] + x[i][k] => 2 * p[k]; );
}}
循环应该做这样的事情：
X [1] [1] + X [2] [1] + X [3] [1] + X [4] [1] + X [5] [1] => 2 * P [1];
X [1] [2] + X [2] [2] + X [3] [2] + X [4] [2] + X [5] [2] => 2 * P [2];
X [1] [3] + X [2] [3] + X [3] [3] + X [4] [3] + X [5] [3] => 2 * P [3];
X [1] [4] + X [2] [4] + X [3] [4] + X [4] [4] + X [5] [4] => 3 * P [4];
但我远离这个结果。
有没有人有一个想法？
谢谢
你可能想使用IloNumExpr
for(int k = 0; k < 4; k++){
IloNumExpr sum_over_i(env);
for(int i = 0; i < 5; i++){
sum_over_i += x[i][k];
}
model.add(sum_over_i >= 2 * p[k]; );
}
您还需要声明x作为一个二维数组。
IloArray x(env, 4);
for (int k = 0; k < 4; ++k)
x[k] = IloIntVarArray(env, 5, 0, 1);
另外，在C ++中，数组的下标是从0到尺寸-1，不为1〜大小。 你的目标应该写成
model.add(IloMaximize(env, 1000 * p[0] + 2000 * p[1] + 500 * p[2] + 1500 * p[3]));
Usertfwr已经给了一个很好的答案，但我想给的解决方案的另一个版本，这可能会帮助你的代码CPLEX应用程序在一个更通用的方法。 首先，我建议你用一个文本文件来保存所有这一切将被送入程序中的数据（目标函数系数）。 在你的情况，你只需要像字面上数据下面的矩阵复制到记事本，并将它命名为“coef.dat”：
[1000，2000，500，1500]
现在到了完整的代码，让我知道如果有任何的理解困难的语句：
#include <ilcplex/ilocplex.h>
#include <fstream>
#include <iostream>
ILOSTLBEGIN
int main(int argc, char **argv) {
IloEnv env;
try {
const char* inputData = "coef.dat";
ifstream inFile(inputData); // put your data in the same directory as your executable
if(!inFile) {
cerr << "Cannot open the file " << inputData << " successfully! " <<endl;
throw(-1);
}
// Define parameters (coef of objective function)
IloNumArray a(env);
// Read in data
inFile >> a;
// Define variables
IloBoolVarArray p(env, a.getSize()); // note that a.getSize() = 4
IloArray<IloBoolVarArray> X(env, 5); // note that you need a 5x4 X variables, not 4x5
for(int i = 0; i < 5; i++) {
X[i] = IloBoolVarArray(env,4);
}
// Build model
IloModel model(env);
// Add objective function
IloExpr objFun (env);
for(int i = 0; i < a.getSize(); i++){
objFun += a[i]*p[i];
}
model.add(IloMaximize(env, objFun));
objFun.end();
// Add constraints -- similar to usertfwr’s answer
for(int i = 0; i < 4; k++){
IloExpr sumConst (env);
for(int j = 0; j < 5; i++){
sumConst += x[j][i];
}
// before clearing sumConst expr, add it to model
model.add(sumConst >= 2*p[i]);
sumConst.end(); // very important to end after having been added to the model
}
// Extract the model to CPLEX
IloCplex cplex(mod);
// Export the LP model to a txt file to check correctness
//cplex.exportModel("model.lp");
// Solve model
cplex.solve();
}
catch (IloException& e) {
cerr << "Concert exception caught: " << e << endl;
}
catch (...) {
cerr << "Unknown exception caught" << endl;
}
env.end();
}
