C#调用CPLEX求解简单线性规划问题样例

这次用的例子是《运筹学 第四版》清华大学出版社 上的例子

using System;
using System.Drawing;
using System.Collections;
using System.ComponentModel;
using System.Windows.Forms;
using System.Data;
using System.Data.SqlClient;
using System.IO;
using System.Text;
using System.Drawing.Drawing2D; 
using System.Drawing.Text;
using System.Runtime.InteropServices;
using System.Collections.Generic;
using ILOG.Concert;
using ILOG.CPLEX;


public class CTransportProblem//节点
{
    public int m_nSupply;//行数
    public int m_nDemand;//列数
    public double[,] Cost=null;
    public double[] x = null;
    public double[] Supply;

    public CTransportProblem()
    {        
        
    }
    public void SetParameters()
    {
        m_nSupply = 3;
        m_nDemand = 3;
        Supply = new double[m_nSupply];
        Cost = Matrix.LoadData("TPCost.txt", '\t');
        //Cost = {{ 1,-2,1 },
        //{ -4,1,2 },
        //{ -2,0,1 }
        //};//决策变量系数
        Supply[0] = 11;Supply[1] = 3;Supply[2] = 1;
        x = new double[m_nDemand];
    }
    public double SolveTransportProblembyGurobi()
    {
        int i, j;
        double[] Coef = {-3,1,1 };//改动
        double MaxKeyNumber = 0;
        Cplex cplexModel = new Cplex();
        //定义变量，逐维度new
        INumVar[] varX = new INumVar[m_nSupply];
        NumVarType[] TypeVar = new NumVarType[m_nDemand];

        //变量下界与上界
        double LbVar = 0.0;
        double UbVar = double.MaxValue;
        for (j = 0; j < m_nDemand; j++)
        {
            TypeVar[j] = NumVarType.Float;
            varX[j] = cplexModel.NumVar( LbVar, UbVar, TypeVar[j]);
        }
        //目标函数
        ILinearNumExpr expr = cplexModel.LinearNumExpr();
       
        for (j = 0; j < m_nDemand; j++)
        {
            expr.AddTerm(varX[j], Coef[j]);
        }
        
        cplexModel.AddMinimize(expr);
        //约束条件       
        expr = cplexModel.LinearNumExpr();
        for (j = 0; j < m_nDemand; j++)
        {
            expr.AddTerm(varX[j], Cost[0,j]);
        }
        cplexModel.AddLe(expr, Supply[0]);

        expr = cplexModel.LinearNumExpr();
        for (j = 0; j < m_nDemand; j++)
        {
            expr.AddTerm(varX[j], Cost[1, j]);
        }
        cplexModel.AddGe(expr, Supply[1]);

        expr = cplexModel.LinearNumExpr();
        for (j = 0; j < m_nDemand; j++)
        {
            expr.AddTerm(varX[j], Cost[2, j]);
        }
        cplexModel.AddEq(expr, Supply[2]);

        //求解
        double objvalue = double.MaxValue;
        if (cplexModel.Solve())
        {
            //最后求的目标函数的值
            objvalue = cplexModel.ObjValue;
            for (j = 0; j < m_nDemand; j++)
            {
                x[j] = cplexModel.GetValue(varX[j]);
            }

            //Matrix.Output(x, '\t', "x.txt");
        }
        cplexModel.End();
        return objvalue;
    }
}

写的比较匆忙，添加约束条件呢里直接复制粘贴的...有些粗糙 见谅...