CPLEXコンサート技術の双対

Question

私は制約

の双対を取得しようとしました。これはC++で実装コードです：プライマル・バリアのようなCPLEXの異なる方法で、このLP問題を解く際、

IloEnv env; 
    IloModel MasterProblem(env); 

    IloNumVarArray XX(env,Create_routes.size(),0,IloInfinity,ILOFLOAT); 
    IloNumVarArray t(env,m,0,IloInfinity,ILOFLOAT); 
    IloExpr expr(env); 

    ////defining ojective of problem 

    IloObjective masterObj(env,expr,IloObjective::Maximize); 
    expr.end(); 
    MasterProblem.add(masterObj); 

    IloRangeArray const1(env); //hala yeki yeki mahdudiyatha ro misazim  

    for (int i=0; i<n; i++){ 
     IloExpr expr(env); 
     for (int j=0; j<Create_routes.size(); j++){ 
      if (Create_routes[j]->internalnodes[i+m]==1) 
       expr+=XX[j]; 
     } 
     const1.add(1==expr); 
     MasterProblem.add(const1[i]); 
     expr.end(); 
    } 
    IloRangeArray const2(env);  
    IloRangeArray const4(env);//mahdudiate depohaye open shode 


    for (i=0; i<m; i++){ 
     IloExpr expr(env); 
     for (int j=0; j<Create_routes.size(); j++){ 
      if (Create_routes[j]->depot==i){ 
       expr+=XX[j]*Create_routes[j]->demand_collected; 
      } 
     } 

     expr-=t[i]*g[i]->QF; 
     const2.add(0>=expr); 
     MasterProblem.add(const2[i]); 
     expr.end(); 
    } 

    IloRangeArray2 const3(env,m); 

    for (i=0; i<m; i++){ 
     const3[i]=IloRangeArray(env); 
    } 

    for (int f=0; f<m; f++){ 
     for (i=0; i<n; i++){ 
      IloExpr expr(env); 
      for (int j=0; j<Create_routes.size(); j++){ 
       if ((Create_routes[j]->depot==f)&&(Create_routes[j]->internalnodes[i+m]==1)){ 
        expr+=XX[j]; 
       } 
      } 
      expr-=t[f]; 
      const3[f].add(0>=expr); 
      MasterProblem.add(const3[f][i]); 
      expr.end(); 
     }  
    } 

    IloCplex cplexM(MasterProblem); 
    cplexM.setParam(IloCplex::RootAlg, IloCplex::Barrier); 
    cplexM.setParam(IloCplex::Threads, 4); 

    if (!cplexM.solve()){ 
     env.error() << "Failed to optimize LP." << endl; 
     nodee->uperbound=0; 
     env.end(); 
     return; 
    } 
    else{ 
     if (!cplexM.isPrimalFeasible()){//agar infeasible bud bia birun 
     nodee->uperbound=0; 
     return; 
     } 
     cout<<"MasterProblem Solved"<<endl; 
     cout<<"objective="<<cplexM.getObjValue()<<endl; 
     javab=cplexM.getObjValue(); 
    } 

    IloNumArray duall(env,n); 
    IloNumArray duall1(env,m); 

    cplexM.getDuals(duall,const1); 
    cplexM.getDuals(duall1,const2); 

    IloNumArray2 duall2(env,m); 
    for (i=0; i<m; i++){ 
     duall2[i]=IloNumArray(env,n); 
     for (j=0;j<n;j++){ 
      duall2[i][j]=cplexM.getDual(const3[i][j]); 
     } 
    } 

デュアル、ネットワーク私は最後に完全に異なるデュアル値と異なるソリューションを得ました。なぜこれはこういうの？ 私の問題に平等の制約があるからですか？ 真の値がプレpleスを通過していることを確認するにはどうすればよいですか？

本当にありがとうございます。

Answer 1

複数の最適なデュアルソリューションがあります。実現可能なすべてのソリューションで、目的の値が最適な目的の値と等しい。これは、同じ制約があってもなくても、固有の最適なプライマル解が存在しても起こる可能性があります。