12 3 / 3 页下一页

查看: 8402|回复: 20

为何公务员招聘16*4矩阵整数规划用LINGO解不出来？

[复制链接]

字体大小: 正常放大

basic_maths

61 主题	7 听众	1127 积分

数学建模

TA的每日心情

	开心 2024-2-3 09:52

签到天数: 171 天

[LV.7]常住居民III

群组: 高等数学课程学习答疑

群组: 中学生数学建模（HIMC

群组: 小波/PDE/图像分割

群组: 专套（升）本数学答疑

群组: 数学建摸协会

电梯直达

1^#

发表于 2009-8-20 10:37 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

本帖最后由 basic_maths 于 2009-8-21 18:45 编辑

为何公务员招聘的分配模型(图 1.jpg)中的16*4矩阵的0-1整数规划用LINGO的sets和link或直接列出的笨办法都解不出来？解出的结果显示的含义是“维数过高，求不出来”或版本不行，但我用的LINGO 8的版本不算低。据说理论不完善导致大规模整数规划用LINGO解很困难，但这论文中的又是怎么用LINGO解出来了呢？谢谢

1.jpg (60.39 KB, 下载次数: 239)

zan

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

使用道具举报

东方明珠-WDZYQ

16 主题	21 听众	3372 积分

无

升级 45.73%

TA的每日心情

	开心 2013-10-25 15:56

签到天数: 4 天

[LV.2]偶尔看看I

自我介绍: 活泼，爱交友

群组: 数学建模

群组: 数学趣味、游戏、IQ等

群组: 我行我数

群组: Mathematica研究小组

群组: LINGO

2^#

发表于 2009-8-20 11:27 |只看该作者 |招呼Ta 关注Ta

报错正茬啊的

我也说一句

发表

使用道具举报

lirui0081

5 主题	4 听众	711 积分

升级 27.75%

TA的每日心情

	开心 2013-8-8 20:47

签到天数: 3 天

[LV.2]偶尔看看I

3^#

发表于 2009-8-20 23:03 |只看该作者 |招呼Ta 关注Ta

就是楼主问题都没放出来

我也说一句

发表

使用道具举报

basic_maths

61 主题	7 听众	1127 积分

数学建模

TA的每日心情

	开心 2024-2-3 09:52

签到天数: 171 天

[LV.7]常住居民III

群组: 高等数学课程学习答疑

群组: 中学生数学建模（HIMC

群组: 小波/PDE/图像分割

群组: 专套（升）本数学答疑

群组: 数学建摸协会

4^#

发表于 2009-8-21 18:47 |只看该作者 |招呼Ta 关注Ta

学习还是很辛苦，论文的关键部分常常看不到，琢磨很久常常陷入困境

我也说一句

发表

使用道具举报

jinling84364336

31 主题	4 听众	1110 积分

升级 11%

TA的每日心情

	开心 2011-10-29 16:57

签到天数: 1 天

[LV.1]初来乍到

群组: 数模讨论——图论方面

群组: 南京邮电大学数模协会

群组: 数学天下

5^#

发表于 2009-8-21 19:37 |只看该作者 |招呼Ta 关注Ta

破解的变量没限制的。估计你的版本是学生版或试用版本。

我也说一句

发表

使用道具举报

basic_maths

61 主题	7 听众	1127 积分

数学建模

TA的每日心情

	开心 2024-2-3 09:52

签到天数: 171 天

[LV.7]常住居民III

群组: 高等数学课程学习答疑

群组: 中学生数学建模（HIMC

群组: 小波/PDE/图像分割

群组: 专套（升）本数学答疑

群组: 数学建摸协会

6^#

发表于 2009-8-21 19:39 |只看该作者 |招呼Ta 关注Ta

本帖最后由 basic_maths 于 2009-8-21 19:40 编辑

运行结果说维数过高的LINGO程序如下：

sets:
set1/1..4/:i;
set2/1..16/:j;
link(set1,set2):h,x;
endsets

max=@sum(link(i,j):h(i,j)*x(i,j));

@for(set1(i)

sum(link(i,j):x(i,j))=1);
@for(set2(j)

sum(link(i,j):x(i,j))<1);
@for(link(i,j)

bin(x(i,j)));

data:
h = 0.078838 0.060075 0.052251 0.065991  0.052251  0.035797 0.065911 0.052251  0.045319 0.029068  0.040378  0.077803  0.050714  0.034983
            0.065991  0.052251
0.080657 0.051387  0.058591  0.06237  0.058591  0.037773  0.06237  0.058591  0.051288  0.034763  0.058536  0.084301 0.064077  0.045746                      0.06237  0.058591
            0.091519  0.056818  0.077022  0.060233  0.077022  0.02591  0.060233  0.077022  0.04915  0.040194 0.034811 0.067051 0.040352 0.043609  0.060233 0.077022
   0.1168  0.064332  0.099849  0.075315  0.099849  0.03862  0.075315  0.099849  0.069347  0.05538  0.067568  0.097246  0.070552  0.066363  0.075315  0.099849;
enddata

结果为何算不出来，谢谢

我也说一句

发表

使用道具举报

Mar

0 主题	0 听众	25 积分

升级 21.05%

该用户从未签到

7^#

发表于 2009-8-21 20:38 |只看该作者 |招呼Ta 关注Ta

sets:
set1/1..16/;
set2/1..4/;
link(set1,set2):h,x;
endsets

max=@sum(link(i,j):h*x);

@for(set2(j)

sum(set1(i):x(i,j))>=1);
@for(set2(j)

sum(set1(i):x(i,j))<=2);
@for(set1(i)

sum(set2(j):x(i,j))>=0);
@for(set1(i)

sum(set2(j):x(i,j))<=4);
@sum(link(i,j) | j #le# @size(set2)

sum(link(i,j) | i #le# @size(set1):x(i,j)))=8;
@for(set1(i):
@for(set2(j): @bin(x(i,j))));

data:
h = 0.078838 0.060075 0.052251 0.065991  0.052251  0.035797 0.065911 0.052251  0.045319 0.029068  0.040378  0.077803  0.050714  0.034983  0.065991  0.052251
0.080657 0.051387  0.058591  0.06237  0.058591  0.037773  0.06237  0.058591  0.051288  0.034763  0.058536  0.084301 0.064077  0.045746  0.06237  0.058591
0.091519  0.056818  0.077022  0.060233  0.077022  0.02591  0.060233  0.077022  0.04915  0.040194 0.034811 0.067051 0.040352 0.043609  0.060233 0.077022
0.1168  0.064332  0.099849  0.075315  0.099849  0.03862  0.075315  0.099849  0.069347  0.05538  0.067568  0.097246  0.070552  0.066363  0.075315  0.099849;
enddata

试一下这个吧！

我也说一句

发表

已有 1 人评分	体力	收起理由
madio	+ 1	鼓励参与解决别人的问题！

总评分: 体力 + 1 查看全部评分

使用道具举报

数魔小子

10 主题	4 听众	1262 积分

升级 26.2%

TA的每日心情

	奋斗 2011-12-24 16:45

签到天数: 136 天

[LV.7]常住居民III

自我介绍: 数模爱好者 java 软件开发

群组: 数学建模

群组: 数学趣味、游戏、IQ等

群组: 聊天

群组: 数学建摸协会

群组: 数学建模培训课堂2

8^#

发表于 2009-8-22 10:06 |只看该作者 |招呼Ta 关注Ta

本帖最后由数魔小子于 2009-8-22 10:08 编辑

sets:
set1/1..4/:i;
   set2/1..16/:j;
   link(set1,set2):h,x;
endsets

max=@sum(link(i,j):h(i,j)*x(i,j));

@for(set1(i)：@sum(link(i,j):x(i,j))=1);
@for(set2(j)：@sum(link(i,j):x(i,j))<1);
@for(link(i,j)：@bin(x(i,j)));

data:
h = 0.078838 0.060075 0.052251 0.065991  0.052251  0.035797 0.065911 0.052251  0.045319 0.029068  0.040378  0.077803  0.050714  0.034983
            0.065991  0.052251
0.080657 0.051387  0.058591  0.06237  0.058591  0.037773  0.06237  0.058591  0.051288  0.034763  0.058536  0.084301 0.064077  0.045746                      0.06237  0.058591
            0.091519  0.056818  0.077022  0.060233  0.077022  0.02591  0.060233  0.077022  0.04915  0.040194 0.034811 0.067051 0.040352 0.043609  0.060233 0.077022
   0.1168  0.064332  0.099849  0.075315  0.099849  0.03862  0.075315  0.099849  0.069347  0.05538  0.067568  0.097246  0.070552  0.066363  0.075315  0.099849;
enddata
我用的是LINGO9.0进行的运算，是学生版本还是没算出，中间的错误已改正。