2

私は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];

しかし、私はこの結果から遠く離れています。

誰かアイデアがありますか？

ありがとう

4

2 に答える 2

2

おそらく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を2次元配列として宣言する必要があります。

IloArray x(env, 4);
for (int k = 0; k < 4; ++k)
      x[k] = IloIntVarArray(env, 5, 0, 1);

また、C ++では、配列インデックスは1からサイズではなく0からサイズ1までです。あなたの目的は書かれるべきです

model.add(IloMaximize(env, 1000 * p[0] + 2000 * p[1] + 500 * p[2] + 1500 * p[3]));
于 2013-03-19T14:20:17.427 に答える
2

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();
}
于 2015-02-13T18:15:31.073 に答える