数学建模社区-数学中国

标题: 求lingo高手帮我看下我的lingo程序到底哪里出错啦 [打印本页]

作者: 廖文强 时间: 2012-7-14 19:54
标题: 求lingo高手帮我看下我的lingo程序到底哪里出错啦
是2002年的国赛彩票问题，用lingo8计算了5个小时都没有出结果，用lingo9以上的版本就总是报错。。。跪求高手解答！！跪谢

Model:
Sets:
Ways/1..7/:p,x,ux;
ab/1..6/:a,b;
Price/1..3/:r;
Endsets
Max=@sum(ways:p*ux);
data:
a,b=10,233,4,54,3,17,4,20,2,10,2,10;
enddata

p(1)=@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/@exp(@lgm(n+1));
p(2)=m*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/@exp(@lgm(n+1));
p(3)=m*(n-m-1)*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/@exp(@lgm(n+1));
p(4)=0.5*m*(m-1)*(n-m-1)*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/@exp(@lgm(n+1));
p(5)=0.25*m*(m-1)*(n-m-1)*(n-m-2)*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/@exp(@lgm(n+1));
p(6)=m*(m-1)*(m-2)*(n-m-1)*(n-m-2)*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/(12*@exp(@lgm(n+1)));
p(7)=m*(m-1)*(m-2)*(n-m-1)*(n-m-2)*(n-m-3)*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/(36*@exp(@lgm(n+1)));

@for(ways(i)|i#le#3:x(i)=(1-p(4)*x(4)+p(5)*x(5)+p(6)*x(6)+p(7)*x(7))*r(i)/p(i));
@for(ways(i):ux(i)=1-@exp(-(x(i)/ 630589)^2));
r(1)+r(2)+r(3)=1;
@bnd(0.5,r(1),0.8);
@bnd(600000,x(1),5000000);
@for(ways(i)|i#le#6:a(i)<(x(i)/x(i+1)));
@for(ways(i)|i#le#6: b(i)>(x(i)/x(i+1)));
@for(ways(i)|i#le#6:p(i)<p(i+1));
@bnd(5,m,7);
@bnd(29,n,60);
@for(ways(i)

bnd(0,p(i),1));
@for(ways(i):x(i)>0);
@for(price(j):r(j)>0);
@gin(m);@gin(n);
end

作者: 廖文强 时间: 2012-7-14 19:55
这个是题目

捕获.JPG (144.22 KB, 下载次数: 199)

捕获.JPG

作者: 廖文强 时间: 2012-7-14 19:56
这个是概率P

捕获2.JPG (69.39 KB, 下载次数: 182)

捕获2.JPG

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