02年彩票的LINGO程序，请教下

model: sets: col/1..10/:n; !n的十种选择1～10对应29～38; row/1..3/:m; !m的3种选择1,2,3对应5，6，7; rank/1..7/:r; !奖的等级1～7; pos(col,row,rank):p; !p(n,m,j)对应n选m的第j等奖概率; dimension/1..2/:nm; !nm(1)=n,nm(2)=m; endsets data: p=2.1051745189676224E-6 0.000012631047113805734 0.0002778830365037262 0.0006947075912593154 0.007294429708222812 0.009725906277630416 0.06483937518420277 6.407052883814503E-7 4.484937018670152E-6 0.0000941836773920732 0.00028255103217621956 0.0028255103217621957 0.0047091838696036595 0.02982483117415651 2.329837412296183E-7 1.8638699298369464E-6 0.000037277398596738924 0.00013047089508858625 0.0012394735033415692 0.0024789470066831384 0.014873682040098831 1.684139615174098E-6 0.000010104837691044588 0.0002324112668940255 0.0005810281672350638 0.006391309839585702 0.008521746452780936 0.05965222516946655 4.912073877591119E-7 3.438451714313783E-6 0.00007564593771490323 0.0002269378131447097 0.0023828470380194517 0.003971411730032419 0.026476078200216132 1.7085474356838674E-7 1.366837948547094E-6 0.000028703596919488972 0.0001004625892182114 0.001004625892182114 0.002009251784364228 0.012725261300973446 1.3581771090113694E-6 8.149062654068216E-6 0.00019557750369763718 0.000488943759244093 0.005622853231307069 0.007497137641742759 0.05497900937278023 3.802895905231834E-7 2.6620271336622838E-6 0.00006122662407423253 0.0001836798722226976 0.0020204785944496735 0.003367464324082789 0.023572250268579523 1.2676319684106113E-7 1.014105574728489E-6 0.00002231032264402676 0.00007808612925409366 0.0008199043571679834 0.0016398087143359669 0.010932058095573112 1.1035189010717375E-6 6.6211134064304255E-6 0.00016552783516076065 0.0004138195879019016 0.004965835054822819 0.006621113406430425 0.050761869449299925 2.97101242596237E-7 2.0797086981736593E-6 0.00004991300875616782 0.00014973902626850346 0.0017219988020877898 0.0028699980034796496 0.02104665202551743 9.507239763079586E-8 7.605791810463668E-7 0.000017493321164066438 0.00006122662407423253 0.0006734928648165578 0.0013469857296331156 0.009428900107431809 9.028791008768762E-7 5.417274605261257E-6 0.00014084913973679269 0.0003521228493419817 0.004401535616774772 0.005868714155699695 0.04694971324559756 2.3407976689400493E-7 1.6385583682580346E-6 0.00004096395920645086 0.0001228918776193526 0.0014747025314322312 0.002457837552387052 0.0188434212349674 7.202454365969383E-8 5.761963492775506E-7 0.000013828712382661215 0.00004840049333931425 0.0005566056734021139 0.0011132113468042277 0.008163549876564338 7.435474948397804E-7 4.461284969038682E-6 0.00012045469416404442 0.00030113673541011103 0.003914777560331444 0.005219703413775258 0.043497528448127154 1.858868737099451E-7 1.3012081159696156E-6 0.00003383141101521001 0.00010149423304563002 0.0012686779130703753 0.0021144631884506253 0.016915705507605003 5.50775922103541E-8 4.406207376828328E-7 0.000011015518442070821 0.00003855431454724787 0.0004626517745669745 0.000925303549133949 0.007093993876693608 6.160822100101038E-7 3.6964932600606224E-6 0.00010350181128169743 0.00025875452820424356 0.0034931861307572884 0.0046575815076763845 0.040365706399862 1.4870949896795608E-7 1.0409664927756925E-6 0.000028106095304943697 0.00008431828591483109 0.0010961377168928043 0.0018268961948213404 0.015224134956844503 4.248842827655888E-8 3.39907426212471E-7 8.837593081524247E-6 0.00003093157578533486 0.0003866446973166858 0.0007732893946333716 0.006186315157066973 5.134018416750865E-7 3.0804110500505186E-6 0.00008933192045146505 0.0002233298011286626 0.0031266172158012765 0.004168822954401702 0.03751940658961532 1.1979376305752016E-7 8.385563414026412E-7 0.000023479577559273953 0.00007043873267782186 0.0009509228911505951 0.0015848714852509918 0.013735552872175262 3.3046555326212464E-8 2.643724426096997E-7 7.138055950461891E-6 0.000024983195826616622 0.0003247815457460161 0.0006495630914920322 0.005413025762433602 4.3014748897101836E-7 2.5808849338261104E-6 0.0000774265480147833 0.00019356637003695827 0.002806712365535895 0.00374228315404786 0.034927976104446694 9.713007815474609E-8 6.799105470832226E-7 0.000019717405865413456 0.00005915221759624037 0.0008281310463473651 0.0013802184105789418 0.012421965695210478 2.5901354174598956E-8 2.0721083339679165E-7 5.801903335110166E-6 0.000020306661672885582 0.00027413993258395533 0.0005482798651679107 0.004751758831455226 3.6222946439664705E-7 2.1733767863798825E-6 0.00006737468037777636 0.0001684367009444409 0.002526550514166613 0.0033687340188888176 0.0325644288492585 7.923769533676654E-8 5.546638673573658E-7 0.000016639916020720974 0.00004991974806216292 0.0007238363469013624 0.0012063939115022706 0.011259676507354526 2.0448437506262334E-8 1.6358750005009867E-7 4.744037501452861E-6 0.000016604131255085016 0.0002324578375711902 0.0004649156751423804 0.0064915675142380; enddata init: x=5E6 6E4 5E3 4E2 30 10 5; r=0.8 0.1 0.1 0 0 0 0; nm=1; !初始化n=1,m=1即29选5; endinit @gin(nm(1)); max=@sum(rank(i):p(nm(1),nm(2),i)*(1-@exp(-1*(x(i)/630589)^2))); r(1)+r(2)+r(3)=1; r(1)<0.8; r(1)>0.5; x(1)<5E6; x(1)>6E5; x(4)>=100; @for(rank(j):x(j)>=0); !@for(rank(j)|j#le#6:x(j+1)/x(j)<=23); x(2)/x(3)<20; x(2)/x(3)>5; @for(rank(i)|i#le#6:x(i+1)/x(i)<10); @for(rank(i)|i#le#6:x(i)>=1); @for(rank(j):r(j)>0); @for(rank(j)|j#le#3:x(j)=(1-@sum(rank(j)|j#ge#4:p(nm(1),nm(2),j)*x(j)))*r(j)/p(nm(1),nm(2),j)); end

2002高教社杯全国大学生数学建模竞赛 B题 彩票中的数学 参考答案 注意：以下答案是命题人给出的，仅供参考。各评阅组应根据对题目的理解及学生的解答，自主地进行评阅。 评价一个方案的优劣，或合理性如何，主要取决于彩票公司和广大彩民两方面的利益。事实上，公司和彩民各得销售总额的50%是确定的，双方的利益主要就取决于销售总额的大小，即双方的利益都与销售额成正比。因此，问题是怎样才能有利于销售额的增加？即公司采用什么样的方案才能吸引广大的彩民积极踊跃购买彩票？具体地讲，问题涉及到一个方案的设置使彩民获奖的可能性有多大、奖金额有多少、对彩民的吸引力有多大、广大彩民如何看待各奖项的设置，即彩民的心理曲线怎样？另外，一个方案对彩民的影响程度可能与区域有关，即与彩民所在地区的经济状况以及收入和消费水平有关。为此，我们要考查一个方案的合理性问题，需要考虑以上这些因素的影响，这是我们建立模型的关键所在。 1. 模型假设与符号说明 彩票摇奖是公平公正的，各号码的出现是随机的； 彩民购买彩票是随机的独立事件； 对同一方案中高级别奖项的奖金比例或奖金额不应低于相对低级别的奖金比例或奖金额； 根据我国的现行制度，假设我国居民的平均工作年限为T =35年。 ---第 等（高项）奖占高项奖总额的比例， ； ----第 等奖奖金额均值， ； ----彩民中第 等奖 的概率， ； ----彩民对某个方案第 等奖的满意度，即第 等奖对彩民的吸引力， ； ----某地区的平均收入和消费水平的相关因子，称为“实力因子”，一般为常数； ----彩票方案的合理性指标，即方案设置对彩民吸引力的综合指标； 2. 模型的准备 （1）彩民获各项奖的概率 从已给的29种方案可知，可将其分为四类, :10选6+1（6+1/10）型、 : 选 型、 : 选 型和 : 选 无特别号型，分别给出各种类型方案的彩民获各奖项的概率公式： &#61548; ：10选6+1（6+1/10）型 ， ， , &#61548; ： 选 型 ， ， ， ， ， ， 。 &#61548; ： 选 型 ， ， ， ， ， ， 。 &#61548; ： 选 无特别号型 ， ， ， ， 。 各种方案的各个奖项获奖概率及获奖总概率 计算如表一。 表一： 方案 6+1/10 2×10-7 8×10-7 1.8×10-5 2.61×10-4 3.42×10-3 4.1995×10-2 ----- 0.045695 7/29 6.40705×10-7 4.48494×10-6 9.4184×10-5 2.8255×10-4 2.8255×10-3 4.7092×10-3 0.029825 0.037742 6+1/29 6.40705×10-7 1.4096×10-5 8.4573×10-5 8.8880×10-4 2.2200×10-3 1.4800×10-2 0.019734 0.037742 7/30 4.91207×10-7 3.43845×10-6 7.5646×10-5 2.2694×10-4 2.3828×10-3 3.9714×10-3 0.026476 0.033137 7/31 3.80290×10-7 2.66203×10-6 6.1227×10-5 1.8368×10-4 2.0205×10-3 3.3675×10-3 0.023572 0.029208 7/32 2.97101×10-7 2.07971×10-6 4.09913×10-5 1.4974×10-4 1.722×10-3 2.8700×10-3 0.021047 0.025832 7/33 2.34080×10-7 1.63856×10-6 4.0964×10-5 1.2289×10-4 1.4747×10-3 2.4578×10-3 0.018843 0.022941 7/34 1.85887×10-7 1.30121×10-6 3.3831×10-5 1.0149×10-4 1.2687×10-3 2.1145×10-3 0.016916 0.020436 7/35 1.48709×10-7 1.04097×10-6 2.8106×10-5 8.4318×10-5 1.0961×10-3 1.8269×10-3 0.015224 0.018261 7/36 1.19794×10-7 8.38556×10-7 2.3480×10-5 7.0439×10-5 9.5092×10-4 1.5849×10-3 0.013736 0.016367 6+1/36 1.19794×10-7 3.47402×10-6 2.0844×10-5 2.9182×10-4 7.2954×10-4 6.5659×10-3 0.008755 0.016367 7/37 9.71301×10-8 6.79911×10-7 1.9717×10-5 5.9152×10-5 8.2813×10-4 1.3802×10-3 0.012422 0.014710 6/40 2.6053×10-7 1.5632×10-6 5.1584×10-5 1.2896×10-4 2.0634×10-3 2.7512×10-3 0.028428 0.033425 5/60 1.831×10-7 9.155×10-7 4.9437×10-5 9.8874×10-5 2.6202×10-3 2.6202×10-3 0.045416 0.050806 （2）确定彩民的心理曲线 一般说来，人们的心理变化是一个模糊的概念。在此，彩民对一个方案的各个奖项及奖金额的看法（即对彩民的吸引力）的变化就是一个典型的模糊概念。由模糊数学隶属度的概念和心理学的相关知识，根据人们通常对一件事物的心理变化一般遵循的规律，不妨定义彩民的心理曲线为 其中 表示彩民平均收入的相关因子，称为实力因子，一般为常数。 （3）计算实力因子 实力因子是反应一个地区的彩民的平均收入和消费水平的指标，确定一个地区的彩票方案应该考虑所在地区的实力因子，在我国不同地区的收入和消费水平是不同的，因此，不同地区的实力因子应有一定的差异，目前各地区现行的方案不尽相同，要统一来评估这些方案的合理性，就应该对同一个实力因子进行研究。为此，我们以中等地区的收入水平（或全国平均水平）为例进行研究。根据相关网站的统计数据，不妨取人均年收入1.5万元，按我国的现行制度，平均工作年限T=35年，则人均总收入为52.5万元，于是，当 万元时，取 （即吸引力的中位数），则有 。 同理，可以算出年收入1万元、2万元、2.5万元、3万元、4万元、5万元、10万元的实力因子如表二。 表二： 年收入 指标 1 万元 1.5 万元 2 万元 2.5 万元 3 万元 4 万元 5 万元 10 万元 420393 630589 840786 1050982 1261179 1681571 2101964 4203928 3．模型的建立与求解 问题（一） 要综合评价这些方案的合理性，应该建立一个能够充分反应各种因素的合理性指标函数。因为彩民购买彩票是一种风险投资行为，为此，我们根据决策分析的理论，考虑到彩民的心理因素的影响，可取 为风险决策的益损函数，于是作出如下的指标函数： （1） 即表示在考虑彩民的心理因素的条件下，一个方案的奖项和奖金设置对彩民的吸引力。 另一方面，由题意知，单注所有可能的低项奖金总额为 ，根据高项奖的计算公式得单注可能的第 项（高项奖）奖金额为 ， 故平均值为 （2） 于是由（1），（2）式得 （3） 利用Matlab可计算出29种方案的合理性指标值 及高项奖的期望值，排在前三位的如下表三。 表三： 指 标 方 案 排 序 9 7/30 4.009×10-7 1.086×106 20679 1410 1 11 7/31 3.784×10-7 1.704×106 32448 2116 2 5 7/29 3.637×10-7 7.557×105 35984 1714 3 问题(二) 根据问题（一）的讨论，现在的问题是取什么样的方案 / （ 和 取何值）、设置哪些奖项、 高项奖的比例 为多少和低项奖的奖金额 为多少时，使目标函数 有最大值。 设以 ， ， ， 为决策变量，以它们之间所满足的关系为约束条件，则可得到非线性规划模型： 关于约束条件的说明： 1）条件（1）（2）同问题（一）； 2）条件（3）（4）是对高项奖的比例约束， 的值不能太大或太小，（4）是根据已知的方案确定的； 3）条件（5）是根据题意中一等奖的保底额和封顶额确定的； 4）条件（6）中的 分别为 等奖的奖金额 比 等奖的奖金额 高的倍数，可由问题（一）的计算结果和已知各方案的奖金数额统计得： 5）条件（7）是根据实际问题确定的，实际中高等奖的概率 应小于低等奖的概率 ，它的值主要有 确定。 6）条件（8）（9）是对方案中 取值范围的约束，是由已知的方案确定的； 这是一个较复杂的非线性（整数）规划，其中概率 的取值分为四种不同的情况 ，者有整数变量 确定，一般的求解是困难的。为此，利用Matlab可求解得最优解为 ，最优值为 。故对应的最优方案为：32选6（6/32），一、二、三等奖的比例分别为80%、9%、11%，四、五、六、七等奖的金额分别为200、10、1、0元。 前面是针对中等收入水平的彩民情况考虑的，对于经济发达地区和欠发达地区应有所不同。这里分别对年收入 万元、2万元、2.5万元、3万元、4万元、5万元、10万元,工作年限均35年的情况进行了讨论，给出适用于相应各种情况的最优方案，如下面的表四。 表四： 年收入 指标 1 万元 2 万元 2.5 万元 3 万元 4 万元 5 万元 10 万元 420393 840786 1050982 1261179 1681571 2101964 4203928 最优方案 5+1/33 6/32 7/30 6/37 6+1/32 7/33 7/35 8.255×10-7 4.623×10-7 4.103×10-7 3.223×10-7 2.475×10-7 2.075×10-7 1.828×10-7 0.80 0.80 0.73 0.70 0.73 0.73 0.80 0.10 0.9 0.17 0.15 0.19 0.18 0.13 0.10 0.11 0.10 0.15 0.07 0.09 0.07 6.5×105 6．18×105 1.38×106 1．46×106 2.23×106 2.99×106 3.91×106 3037 120004 47506 52172 22721 1.07×105 94252 607 600 1235 1739 1507 1974 1746 138 200 100 200 100 200 103 7 10 10 20 20 10 20 1 1 5 2 2 2 5 0 0 0 0 0 0 3 问题（三）（略） 说明： （1）研究此问题必须要考虑心理曲线，但心理曲线的可能会有不同的形式，主要是看对问题解释是否合理，实力因子 在不同地区可以取不同的值，对方案的评判结果也会有差别。 （2）问题的合理性指标函数的一定与心理曲线有关，但应该在风险决策的意义下确定出益损函数，益损函数的确定不是唯一的。 （3）问题中的概率公式的形式应该是唯一的。 参考文献 中国彩票网：http://www.cpiao.com

不能上传文档~~

恐怖 【 在 zhaoweideyx (DemHt) 的大作中提到: 】 : 这是我们组写的02年彩票的LINGO程序， : 将根据n,m,i(n选m，中i等奖的概率)做成表格（二维的），n一维，m+i一维,换算了一下。 : 然后就去做求解，根据m,n,查P表，得一个概率值这样子的，然后接着算。 : ...................

还得等几天才能申请附件。。 【 在 zhaoweideyx (DemHt) 的大作中提到: 】 : 不能上传文档~~

待会改下程序~~

数据还是放文本文件里吧。。。 【 在 zhaoweideyx (DemHt) 的大作中提到: 】 : 待会改下程序~~

偷了一下懒~~

上次解一题数据量巨大，，贴进去lingo直接罢工了 【 在 zhaoweideyx (DemHt) 的大作中提到: 】 : 偷了一下懒~~