查看: 4392|回复: 9

求教lingo整数规划问题

字体大小: 正常放大

1 主题	7 听众	2 积分

升级 40%

该用户从未签到

自我介绍: 软件工程师

电梯直达

1^#

发表于 2015-7-16 11:22 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

本帖最后由 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

zan

lingo

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

相关帖子

使用道具举报

liwenhui

70 主题	65 听众	5197 积分

独孤求败

TA的每日心情

	擦汗 2018-4-26 23:29

签到天数: 1502 天

[LV.Master]伴坛终老

自我介绍: 紫薇软剑，三十岁前所用，误伤义士不祥，乃弃之深谷。重剑无锋，大巧不工。四十岁前恃之横行天下。四十岁后，不滞于物，草木竹石均可为剑。自此精修，渐进至无剑胜有剑之境。

群组: 计量经济学之性

群组: LINGO

2^#

发表于 2015-7-16 13:04 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

数学上讲，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主管

715 主题	213 听众	8601 积分

TA的每日心情

	开心 2017-4-28 17:18

签到天数: 415 天

[LV.9]以坛为家II

群组: 乐考无忧考研公益讲座

群组: 2017美赛两天强训

群组: 模友会交流视频

群组:

群组: 国赛讨论

3^#

发表于 2015-7-16 15:52 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

稍等会有有人给你解答问题的

我也说一句

发表

使用道具举报

森之张卫东

413 主题	36 听众	1854 积分

升级 85.4%

TA的每日心情

	开心 2019-9-18 21:55

签到天数: 258 天

[LV.8]以坛为家I

群组: 2015国赛冲刺

群组: 2016美赛公益课程

群组: 国赛讨论

群组: 第三届数模基础实训

群组: Matlab讨论组

4^#

发表于 2015-7-16 16:35 |只看该作者 |招呼Ta 关注Ta

!我用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