Java自学者论坛

 找回密码
 立即注册

手机号码,快捷登录

恭喜Java自学者论坛(https://www.javazxz.com)已经为数万Java学习者服务超过8年了!积累会员资料超过10000G+
成为本站VIP会员,下载本站10000G+会员资源,会员资料板块,购买链接:点击进入购买VIP会员

JAVA高级面试进阶训练营视频教程

Java架构师系统进阶VIP课程

分布式高可用全栈开发微服务教程Go语言视频零基础入门到精通Java架构师3期(课件+源码)
Java开发全终端实战租房项目视频教程SpringBoot2.X入门到高级使用教程大数据培训第六期全套视频教程深度学习(CNN RNN GAN)算法原理Java亿级流量电商系统视频教程
互联网架构师视频教程年薪50万Spark2.0从入门到精通年薪50万!人工智能学习路线教程年薪50万大数据入门到精通学习路线年薪50万机器学习入门到精通教程
仿小米商城类app和小程序视频教程深度学习数据分析基础到实战最新黑马javaEE2.1就业课程从 0到JVM实战高手教程MySQL入门到精通教程
查看: 858|回复: 0

1.1python解决数学建模之席位分配问题

[复制链接]
  • TA的每日心情
    奋斗
    2024-11-24 15:47
  • 签到天数: 804 天

    [LV.10]以坛为家III

    2053

    主题

    2111

    帖子

    72万

    积分

    管理员

    Rank: 9Rank: 9Rank: 9

    积分
    726782
    发表于 2021-9-2 13:52:10 | 显示全部楼层 |阅读模式
    一:上代码
    #比例法
    def rate_method(p,n):
        lst =[] #保存各组席位数
        sum_ =sum(p)    #人数和
        k =0#临时变量
        for i in p:
            lst.append(i/sum_*n)
            k += int(i/sum_*n)
        while k!=n:
            max_ =0
            for i in lst:
                max_ =i if (i -int(i)) > max_- int(max_)else max_  #小鼠的比较
            lst[lst.index(max_)] =int(max_) +1
            k +=1
        for i in lst:
            lst[lst.index(i)] =int(i)  
        print(lst)
               
               
    #Q值法:
    def q_value(p,n):       #p:保存各组人数的列表,n:席位数
        lst =[] #保存各组席位数
        for i in p:
            lst.append(1) #初始席位数都为1
        lst_ =lst[:]     #临时列表,保存各组Q值
        for i in range(n-len(p)):
            for i in range(len(p)):
                lst_ =p*p/(lst+1)/lst
            max_index=lst_.index(max(lst_))
            lst[max_index] +=1
        print(lst)
    #d'Hondt 方法
    def d_Hondt(p,n):       #p必须以大到小顺序排序
        lst =[] #保存各组席位数
        a,b=0,0 #保存比值,下一项比值
        for i in p:
            lst.append(0) #初始席位数都为0
        for i in range(n):
            for j in range(len(p)):
                a =p[j]/(i+1)
                if b>a:
                    break
                b =p[j]/(i+2)
                lst[j] += 1
            if sum(lst) ==n:    #放在第一个循环内,是考虑到不太可能只除了1就分配好了
                break
        print(lst)
    if __name__ =='__main__':
        p =[432,333,235] #从大到小比较好,为了一次性使用所有方法
        n =10
        rate_method(p,n)
        q_value(p,n)
        d_Hondt(p,n)

     

    二:席位分配常用的三种方法:

    1.比例分配 :公平而又简单

    2.Q值法:(摘自:数学模型-姜启源-)

            设两方人数分别 p1 和 p2 ,占有席位分别是 n1 和 n2 ,则两方每个席位代表的人数分别为 p1/n1 和 p2/n2 .显然仅当 p1/n1 = p2/n2 时席位的分配才是公平的.但是因为人数和席 位都是整数,所以通常 p1/n1 ≠ p2/n2 ,这时席位分配不公平,并且 pi ni (i= 1, 2)数值较大的一方吃亏,或者说对这一方不公平.

             不妨假设 p1/n1 > p2/n2 ,不公平程度可用数值 p1/n1 - p2/n2 衡量.如设 p1 = 120, p2 = 100, n1 = n2 = 10,则 p1/n1 - p2/n2 = 12 - 10 = 2,它衡量的是不 公平的绝对程度,常常无法区分两种程度明显不同的不公平情况.例如上述双方 人数增加为 p1 = 1 020 和 p2 = 1 000,而席位 n1 , n2 不变时, p1/n1 - p2/n2 = 102 - 100 = 2,即绝对不公平程度不变.但是常识告诉我们,后面这种情况的不公 平程度比起前面来已经大为改善了.

            为了改进上述绝对标准,自然想到用相对标准.仍记 p1 , p2 为 A, B 两方的 固定人数, n1 , n2 为两方分配的席位(可变),若 p1/n1 > p2/n2 ,则定义

      rA ( n1 , n2 ) =(p1 /n1 - p2 /n2)/( p2 /n2)      (1)           为对 A 的相对不公平度.

    若 p2 n2 > p1 n1 ,则定义

      rB ( n1 , n2 ) =(p2 /n2 - p1/ n1 )/(p1 n1)    (2)    为对 B 的相对不公平度

      建立了衡量分配不公平程度的数量指标 rA ,rB 后,制定席位分配方案的原 则是使它们尽可能小. 确定分配方案 假设 A, B 两方已分别占有 n1 和 n2 席,利用相对不公平 度 rA 和 rB 讨论,当总席位增加 1 席时,应该分配给 A 还是 B. 不失一般性可设 p1 /n1 > p2 /n2 ,即对 A 不公平.当再分配 1 个席位时,关 于 pi ni (i= 1,2)的不等式可能有以下 3 种情况: 1 . p1/ ( n1 + 1) > p2 /n2 ,这说明即使 A 方增加 1 席,仍然对 A 不公平,所以 这一席显然应分给 A 方. 2 . p1 /(n1 + 1) < p2 /n2 ,说明当 A 方增加 1 席时将变为对 B 不公平,参照 (2)式可计算出对 B 的相对不公平度为  rB ( n1 + 1, n2 ) =p2 /(n1 + 1) /(p1 /n2)- 1     (3)3 . p1 /n1 > p2/ ( n2 + 1),即当 B 方增加 1 席时将对 A 不公平,参照(1)式可 计算出对 A 的相对不公平度为

      rA ( n1 , n2 + 1) =p1 (/ n2 + 1) (p2/ n1) - 1     (4)

      (不可能出现 p1 /n1 < p2 /( n2 + 1)的情况.为什么 ?) 因为公平分配席位的原则是使得相对不公平度尽可能地小,所以如果 rB ( n1 + 1, n2 ) < rA ( n1 , n2 + 1) (5) 则这 1 席应分给 A 方;反之则分给 B 方.根据(3),(4)两式,(5)式等价于 (p2 ^2) /(n2 (n2 + 1) )< (p1^2 )/ n1 ( n1 + 1) (6) 还不难证明,上述第 1 种情况的 p1 /( n1 + 1) > p2 /n2 也会导致(6)式.于是我们 的结论是:当(6)式成立时增加的 1 席应分给 A 方,反之则分给 B 方.或者,若记 Qi =( p i^2)/ (ni( ni + 1)),i= 1,2,则增加的 1 席应分给 Q 值较大的一方. 上述方法可以推广到有 m 方分配席位的情况.设第 i 方人数为 pi,已占有 ni 个席位,i= 1,2,⋯, m.当总席位增加 1 席时,计算
    Qi =(pi^2)/(ni( ni + 1)), i= 1,2,⋯, m     (7)
    应将这一席分给 Q 值最大的一方.这种席位分配方法称 Q 值法.

    3.d’ Hondt 方法:

    将各组人数用正整数 n = 1, 2, 3,⋯相除,将所得商数从大到小取前 n个(n 为席位数)

     

    哎...今天够累的,签到来了1...
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 立即注册

    本版积分规则

    QQ|手机版|小黑屋|Java自学者论坛 ( 声明:本站文章及资料整理自互联网,用于Java自学者交流学习使用,对资料版权不负任何法律责任,若有侵权请及时联系客服屏蔽删除 )

    GMT+8, 2024-12-22 11:58 , Processed in 0.055533 second(s), 27 queries .

    Powered by Discuz! X3.4

    Copyright © 2001-2021, Tencent Cloud.

    快速回复 返回顶部 返回列表