席位分配问题的D’hondt模型和相对尾数模型席位分配问题的D’hondt模型和相对尾数模型摘要:讨论公平席位分配的模型已有很多。
本文首先用比例加惯例法、Q值法、D’hondt法对问题中名额进行了分配,再对D’hondt法的合理性进行了分析,并在Q值法对绝对尾数(绝对不公平度)的处理方式基础上,提出了相对尾数模型,并讨论了其满足Young公理的1,3,4条;在模型求解上,全部由MATLAB程序来实现名额分配。
关键词:相对尾数 Balinsky & Young不可能定理 MATLAB正文1 问题复述公平的席位分配问题是一个非常有趣而重要的问题,它在政治学、管理学和对策论等领域具有广泛的应用价值。
处理这个问题的最早的方法是Hamilton法,即比例加惯例法;后来出现了Q值法;1974年M.L.Balinski和H.P.Young引入了席位分配问题的公理体系研究方法,并于1982年证明了同时满足五个公理的席位分配方法是不存在的;因此,我们只能根据实际建立在一定公平准则下成立并尽量多的满足Young公理的算法。
这里,我们需要理解并运用比例加惯例法、Q值法、D’hondt法对宿舍委员会名额进行分配,继而提出更优的公平分配席位的方法。
2 模型假设2.1合理假设2.1.1 比例加惯例法、Q值法等分配模型均为已知;2.1.2 各个宿舍相互独立互不影响,人数保持不变;2.1.3 委员分配以各宿舍人数为唯一权重。
2.2 符号约定符号意义Q第个宿舍的Q值 iin第个宿舍的人数 iim第个宿舍分配的名额 iin总人数m总名额数p 第个宿舍的理想分配名额 ii,个宿舍的理想分配名额总席位增加一个时第ip iniq 第个宿舍的分配比例,即 miins 第个宿舍的绝对尾数值 iir 第个宿舍的相对尾数值 ii,个宿舍的相对尾数值总席位增加一席时第ir it按比例分配后剩余名额3 模型的建立与求解3.1按比例加惯例模型分配根据比例加惯例分配模型的原理,编写MATLAB程序实现(附录-程序1,2,3,附录-输入及运行结果1),结果如表所示:表1(比例加惯例法分配结果):10个席位的分配 15个席位的分配宿舍学生人数比例分配惯例分配比例分配惯例分配的席位的结果的席位的结果A 235 2 3 3 4B 333 3 3 4 5C 432 4 4 6 6总数 1000 9 10 13 153.2按Q值法模型分配2ni首先用比例分配法对名额进行初步分配,再根据表达式 i,A,B,CQ,im(m,1)ii对剩下的名额进行分配,编写MATLAB程序实现求解(附录-程序4,5,附录-输入及运行结果2):表2(Q值法分配结果):10个席位的分配 15个席位的分配宿舍学生人数比例分配最终分配比例分配最终分Q值 Q值名额名额名额配名额A 235 2 9204.17 2 3 4602.08 4B 333 3 9240.75 3 4 5544.45 5C 432 4 9331.2 5 6 4443.43 6 总数 1000 9 10 13 153.3 D’hondt模型3.3.1 模型建立设,分别表示宿舍总人数和总分配席位数,()表示各宿舍人数,令nnmi,1,2,3iniaa,(),则得到一个数列,将该数列按递减顺序重新排列,得ij,,1,2,3,1,2,...,,ijijj()k()k()k()kaaa到,其中表示中第大的项。
取中前项,则相应得到amk,,,,,,ijijijij()k(),m,,m即为按p,1,2,3mmaip,,(k=1,2,...,m)中的元素的个数,,,,123pijD’hondt模型分配的结果。
3.3.2 按D’hondt模型分配根据建立的D’hondt模型,编写MATLAB程序求出结果(附件-程序6,附录-输入及运行结果3):表3(D’hondt模型分配结果):宿舍人数 10个名额的分配 15个名额的分配A 235 2 3B 333 3 5C 432 5 7总数 1000 10 153.4 相对尾数模型3.4.1 模型准备讨论一般情况:个宿舍人数分别为n,,总人数为,待分配nnn,,,...kik,1,2, (i1)k的席位为个,理想化的分配结果是(),满足,记pmp,mik,1,2,...,,iii,1 ni()。
显然,若q全为整数,应有q=p(),当q不全为整数ik,1,2,...,ik,1,2,...,qm,iiiiin时,需要确定同时满足下面公理的分配方案。
pqq公理一:(),即取或之一,其中qpq,,ik,1,2,...,,,,,,,,,iiiiii,,,,qqqqq=,=q,1,表示的整数部分。
,,,,,,,,,,iiiiii,,公理二:,,即总席位增加时,各pmnnnpmnnn(,,,...,)(1,,,...,),,ik,1,2,...,ikik1212宿舍的席位数不应该减少。
公理一显然满足Balinsky & Young不可能定理 (见附录) 中的公理4(公平分摊性),公理nn,,ii二满足其的公理1(人口单调性)和公理3(名额单调性)。
令,smmqq,,,,,,iii,,,nn,,,si称其为对第个宿舍的绝对尾数值。
令,称其为对第个宿舍的相对尾数值。
ii,riq,,i,3.4.2 模型建立与求解r由于人数都是整数,为使分配趋于公平,需所有的越小越好,所以趋于公平的分配方案irr应该是最大的达到最小,即所有的达到最小。
iiqqnn,nnn,,为方便起见,首先考虑只有两个宿舍的情形,即,,且,和k,2121212 ,,pr不全是整数(实际上,他们同为整数或小数)。
记,为总席位增加一席时的分配结果和相ii对尾数。
给出定理:定理:以下分配方案满足公理一,二,nn,,,,121) 若,且,则取,,即按比例加惯例rr,ss,pm1,,pm,121212,,,,nn,,,,,,法分配;nn,,,,122) 若,则取,; rr,pm1,,pm,1212,,,,nn,,,,,,nn,,,,123) 若,则取,。
rr,pm1,pm,,1212,,,,nn,,,,,,定理证明见附录。
按照定理,对三个宿舍的情形进行讨论。
设,,全部为零(实际上,如果有一个为零,rrr123即是按两个宿舍分配),可以做以下分配:1) 当时,按比例分配取整后,剩余的席位分配给绝对尾数较大的宿舍,即按rrr,,123比例加惯例法分配;2) 当时,按比例分配后,若剩余一个席位,则分配给第一个宿舍,若剩余两rrr,,123个席位,则分配一席给第一个宿舍,另外一席分配给第二三个宿舍中绝对尾数值较大者;3) 当时,按比例分配后,若剩余一个席位分配给第一二个宿舍中绝对尾数rrr,,123值较大者,若剩余两个席位,则分配给第一二宿舍各一席;4) 当时,按比例分配后,若剩余一个席位,则分配给第一个宿舍,若剩余两rrr,,123个席位,则分配给第二个宿舍。
一般地,对个宿舍,设r,r,…,r不全为零,且,则当rr,时,rrr,,,...k12n12ktt,1kn,,i将剩余的个席位分配给第一至第个宿舍各一席,当trrrr,,,tmm,,,tttt,,,112,,n,,i,1,kn,,iss时,个席位分配给第一至第个宿舍及和较大的宿舍各一席,t,1tmm,,,tt,1,,n,,i,1,kn,,i当()时,个席位分配给第一至第rrrr,,,t,11,,,skttmm,,,tttts,,,11,,n,,i,1,个宿舍及s,s,…中较大的宿舍各一席,当(),srrr,,1,',,,sskttt,1ts,tststs,,,,1'kn,,iss 个席位分配给第一至第个宿舍及,,…s中个较大的所sts,tmm,,,tt,1ts,,,n,,i,1,对应的宿舍各一席。
最后,编写出尾数法的MATLAB程序,实现3本题中的名额分配(附录-程序7,附录-输入及运行结果4)。
表4(尾数法分配结果):宿舍人数 10个名额的分配 15个名额的分配A 235 3 4B 333 3 5C 432 4 6总数 1000 10 154 模型检验及结果分析席位分配的尾数模型满足Young公理的1、3、4条,是以严格证明了的定理形式给出。
对按上述四种分配模型分配的结果列表比较。
表5(各方法分配结果的比较1):20个席位的分配 21个席位的分配学生宿舍人数 B Q D R B Q D RA 103 10 11 11 10 11 11 11 10B 63 6 6 6 6 7 6 7 7C 34 4 3 3 4 3 4 3 4 总数 200 20 20 20 20 21 21 21 21表6(各方法分配结果的比较2):10个席位的分配 15个席位的分配学生宿舍人数 B Q D R B Q D RA 235 3 2 2 3 4 4 3 4B 333 3 3 3 3 5 5 5 5C 432 4 5 5 4 6 6 7 6 总数 1000 10 10 10 10 15 15 15 15表格中,B表示比例加惯例法,Q表示Q值法,D表示D'hondt法,R表示相对尾数法。
“比例加惯例”法用各团体人数占团体总人数的比例乘以总席位数, 取其整数位为第一次分配, 再次分配时, 则按小数位的大小分, 大的先分配, 直到席位分完。
从表4看到,当总席位数增加时,C宿舍分得的席位却减少;Q值法利用相对不公平度建立了衡量不公平程度的数量指标, 进而将席位分给最不公平的一方。
D’hondt方法将各团体的人数用正整数相除, 其商数组成一个表, 将数从大到小取, 直到取得的商数的个数等于总席位数, 统计出每个团体被取到的商数的个数, 即为该团体分得的席位数。
5 优缺点分析及改进从对模型的检验与分析可以看到,上面讨论的三个模型都有自身的不足:比例加惯例法满足公理一,却不满足公理二;Q值法满足公理二但不满足公理一;D’hondt法也不能解决对每个宿舍成员公平的大小问题;尾数法虽然满足公理一和二,但由于两个公理本身只满足Young公理体系的部分,也不尽完美。
优点:尾数模型打破Q值法的对绝对尾数的比较方法,以相对尾数来讨论,使得模型满足了Young公理体系中更多的公理,虽不尽完善,但相比之前的四种方法是很大的改进。
并且,这种对已有方法改进的思想很有启发意义。
改进:本文中只给出了尾数法对3个宿舍的名额分配程序,对不定数量宿舍的分配没能程序实现,是可以改进的。
参考文献[ 1 ] 姜启源等数学建模[M](第三版)北京高等教育出版社,2004.24—27.[ 2 ] 岳林关于Q值法的一种新定义[J]. 系统工程.1995,13(4):70—73. [ 3 ] 高尚席位分配的最大熵法[J].数学的实践与认识,1996,26(2):73—75. [ 4 ] 吴承祯,洪伟资源公平分配的遗传算法研究[J].运筹与管理,1998,7(2):23—28.2[ 5 ] 吴黎军名额分配问题中的拟合法[J].生物数学报,1995,10(3):77—81. ,[ 6 ] 严余松席位分配问题的0-1规划模型[J].系统工程,1996,14(5):51—53. [ 7 ] 林建良席位分配的最小极差法[J].华南理工大学学报,2001,29(1):21—23. [ 8 ] 杜跃鹏杜太生席位分配的最大概率法[J].数学的实践与认识,2003,33(7):15—19. [ 9 ] 王秀莲席位分配问题的相对尾数法[J].数学的实践与认识,2007,37(9):81—85.附录Balinsky & Young不可能定理公理 1 (份额单调性) 一个州人口的增加不会导致它失去席位。