数学建模社区-数学中国

标题: 求教lingo整数规划问题 [打印本页]

作者: daoliyouyuan 时间: 2015-7-16 11:22
标题: 求教lingo整数规划问题
本帖最后由 daoliyouyuan 于 2015-7-16 11:24 编辑

我用lingo写了一个求解背包问题的程序，但运行时提示：
Improper use of attribute name.
23] @bin(y)
请教这个问题怎么解决? 多谢...

model:
sets:
cheque/1..2/: y, c;
order/1..2/: o;
links(cheque, order):x;
endsets

data:
c = 2500, 1500;
o = 1000, 700;
enddata

min = @sum(cheque: y);

@sum(order(j): o(j) * x(1, j)) <= c(1) * y(1);
@sum(order(j): o(j) * x(2, j)) <= c(2) * y(2);

@for(cheque(i) : @sum(order(j): o(j) * x(i, j)) <= c(i) * y(i) );
@for(order(j) : @sum(cheque(i): x(i,j)) = 1);

@for(links: @bin(x));

@bin(y);

end

作者: liwenhui 时间: 2015-7-16 13:04
数学上讲，y是一个向量，你需要用一个@for来表示对每一个y的分量都限定为0-1变量。

改为：

model:
sets:
cheque/1..2/: y, c;
order/1..2/: o;
links(cheque, order):x;
endsets

data:
c = 2500, 1500;
o = 1000, 700;
enddata

min = @sum(cheque: y);

@sum(order(j): o(j) * x(1, j)) <= c(1) * y(1);
@sum(order(j): o(j) * x(2, j)) <= c(2) * y(2);

@for(cheque(i) : @sum(order(j): o(j) * x(i, j)) <= c(i) * y(i) );
@for(order(j) : @sum(cheque(i): x(i,j)) = 1);

@for(links: @bin(x));

@for(cheque

bin(y));

end

运行结果：

  Global optimal solution found.
  Objective value:                            1.000000
  Objective bound:                            1.000000
  Infeasibilities:                            0.000000
  Extended solver steps:                            0
  Total solver iterations:                            0

  Model Class:                                     PILP

  Total variables:                   6
  Nonlinear variables:                0
  Integer variables:                   6

  Total constraints:                   7
  Nonlinear constraints:             0

  Total nonzeros:                   18
  Nonlinear nonzeros:                0

                              Variable          Value       Reduced Cost
                                 Y( 1)       1.000000          1.000000
                                 Y( 2)       0.000000          1.000000
                                 C( 1)       2500.000          0.000000
                                 C( 2)       1500.000          0.000000
                                 O( 1)       1000.000          0.000000
                                 O( 2)       700.0000          0.000000
                              X( 1, 1)       1.000000          0.000000
                              X( 1, 2)       1.000000          0.000000
                              X( 2, 1)       0.000000          0.000000
                              X( 2, 2)       0.000000          0.000000

                                 Row Slack or Surplus    Dual Price
                                    1       1.000000          -1.000000
                                    2       800.0000          0.000000
                                    3       0.000000          0.000000
                                    4       800.0000          0.000000
                                    5       0.000000          0.000000
                                    6       0.000000          0.000000
                                    7       0.000000          0.000000

作者: 数学中国YY主管 时间: 2015-7-16 15:52
稍等会有有人给你解答问题的

作者: 森之张卫东 时间: 2015-7-16 16:35
!我用lingo写了一个求解背包问题的程序，但运行时提示：
Improper use of attribute name.
23] @bin(y)
请教这个问题怎么解决? 多谢...;

model:
sets:
   cheque/1..2/: y, c;
   order/1..2/: o;
   links(cheque, order):x;
endsets

data:
   c = 2500, 1500;
   o = 1000, 700;
enddata

min = @sum(cheque: y);

@sum(order(j): o(j) * x(1, j)) <= c(1) * y(1);
@sum(order(j): o(j) * x(2, j)) <= c(2) * y(2);

@for(cheque(i) : @sum(order(j): o(j) * x(i, j)) <= c(i) * y(i) );
@for(order(j) : @sum(cheque(i): x(i,j)) = 1);

@for(links: @bin(x));

!@bin(y);
!用这种方法，限制y的取值为0或1;
@for(cheque: @bin(y));
!或者这种。
!y(1)*(y(1)-1)=0;
!y(2)*(y(2)-1)=0;

end

作者: 数学中国YY主管 时间: 2015-7-16 17:49
楼上惊险大神~~~ 楼主有福了

作者: 数学中国YY主管 时间: 2015-7-16 17:54
就是需要在结尾处加一个for 语句循环。。因为你的Y定义了2个值明白了么？

作者: wujianjack2 时间: 2015-7-16 19:13
二楼给的解答是值得推荐的，言简意赅。
四楼给的注释中的另一种方案不好，指定决策变量为binary不必整成复杂的complementary constraints，能简化模型一般经验以简化模型为先，而不是复杂化。

作者: 1345389794 时间: 2015-7-16 23:54
真的有福·····

作者: 1345389794 时间: 2015-7-19 23:28
大神················

作者: 我爱我的晓~~ 时间: 2018-1-21 18:33
来刷个体力

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)