数学建模社区-数学中国
标题:
席位分配问题——惯例Q值法和d'hondt法的MATLAB程序
[打印本页]
作者:
佛自业障
时间:
2018-10-30 10:12
标题:
席位分配问题——惯例Q值法和d'hondt法的MATLAB程序
程序介绍了数学建模中经典问题的两种解法，即席位分配问题！
A( ^: t/ S0 p* P$ R% @
%适用于所有情况 BY Gu
1 |; z+ G2 W7 I0 L+ Y Y
clear all
. q) r5 Z! K% n. ~& ]
clc %惯例Q值法分配席位，wy为席位数，ps为人数总和，R为分配方案
5 r7 ^" h5 S: W$ I j
wy=19;
$ G) Q/ |' p, m) ?
P=[103 63 34]
( t- m, _' o7 U- c6 s
%菜单选项
5 [+ g8 ^$ x9 p! x2 M, q: m
MENUN=menu('选择方法','惯例Q值法','dhondt法','都用');
r- ]2 ^& [9 G
%惯例Q值法
8 P5 O$ X/ e& w/ Q4 l- ?7 r
if (MENUN==1||MENUN==3)==1
1 x; M9 t3 w. q+ s, ~$ _
n=size(P,2);
# b# e, {' c0 z# T
ps=sum(P);
# q* b7 |0 ?1 J! M8 `! t
RS=0;
1 R1 p$ I9 _$ t* k
for i=1:n
# N2 N1 O, u8 Z) u! [1 i5 l% ]1 h# C
R(i)=fix(P(i)/ps*wy);
2 { [! z A% E- G# n
RS=RS+R(i);
. Q* p$ S5 G; |5 q' @- E9 O F+ Q
end
- N' |5 j% v( G7 V" H9 p
wy1=wy-RS;
3 v( o2 j C! _6 d
for i=1:wy1
, T; J: g* H/ L9 j$ A
for j=1:n
# O4 E) y9 ?5 L
Q(j,i)=P(j).^2/(R(j)*(R(j)+1));
) _4 s ~% ^2 }$ R3 N* z3 r4 \5 k
end
8 b0 d: `; b- X2 g0 H- w3 R8 \# D$ s/ Y
t=Q(1,i);
- I! |9 A/ x! R2 V: Q
a=1;
" E' ~9 B. P* d$ M
for m=2:n
9 t$ F5 |! _. Z% ^2 N
if Q(m,i)>t
7 ?0 n6 {, K7 x# y% l
t=Q(m,i);
1 G6 T+ b$ A2 J3 r+ b }1 C
a=m;
3 A0 [. k- N a4 t: m' Q6 a
end
' d7 L6 `1 x. m; p# M
end
5 T3 e" I$ H& D [% T
R(a)=R(a)+1;
' Y# R! J, k( O( w% j
end
4 n k9 V/ Q7 Z8 x" u3 m
disp('惯例Q值法席位分配为：')
% b5 q* V0 `9 l* [& n
R
0 z% e2 g+ d' S& d7 `; H
end
& g G5 f' R1 k3 P
%d'hondt法………………………………………………………………………………$$
" s! V9 Q$ O6 H3 D/ h* T9 N
if (MENUN==2||MENUN==3)==1
! a4 V' j! g1 s# x
zs=P;
9 Z, _, }0 \3 G4 g
n=size(zs,2);
1 h* Q8 R- ~; t
for i=1:n
' |- c4 U- P _* |
for j=1:wy
1 y3 I2 |4 P1 F( c4 y
r(i,j)=zs(i)/j;
+ Y7 I: m. Q& z9 H% e
end
+ _0 ~, x3 a8 o5 P7 t
end
8 c) n( x0 s$ |3 r2 p9 @, {
t=r(
';
' m4 g. y) B# g% v7 z1 @
T=sort(t,'descend');
0 S& `& C9 `/ p5 C; l0 I
min=T(wy);
8 h# J- S/ f7 W. P, X% h
R1=zeros(1,n);
) ~/ {; K% ?2 ^$ h2 v
for i=1:n
4 J" f3 R& `8 F) X2 H0 R0 v
for j=1:wy
N) L# J1 T7 H0 r1 y9 d
if r(i,j)>min
8 U% @- J% R; c0 {* K* ?8 z3 J2 R
R1(i)=R1(i)+1;
. f' o9 R# s) M$ v3 R
end
6 o. P0 {0 e$ ^% i" ^+ R
end
, X6 S' [5 x* ~7 S! A0 S. R
end
& u3 [- Z( m# Z& W" R
rest=wy-sum(R1);%剩余席位
' R) q9 C3 L9 L; ~1 F
%对剩余席位的分配处理
1 b; |, U! [0 ~8 Y4 q6 L
k=0;
) }7 o# T, j) ~' w* e0 ]
for i=1:n
: z5 R3 q v, E, [. s* b% {
if k>=rest%席位分配完毕，结束循环
% o& p6 U0 L' r. i5 y9 N$ t0 D
break;
2 V6 r( Q+ E( z/ k* K
end
4 V! [; R9 b( j8 q! r" b
for j=1:wy
- z6 o6 L' n; c b, I2 X6 W
if r(i,j)==min
' o6 M. K: e' Y
R1(i)=R1(i)+1;
) { Y3 j$ W4 y% ]1 f" @
k=k+1;%对剩余席位分配计数
9 V" d, B7 o& @2 d4 Z+ @6 j
break;
/ m0 |( P. E+ M; Q$ E; P
end
3 H4 O& D/ b/ B
end
7 G7 P [0 P2 R( W$ P! u# l
end
' G! C! i9 A2 N
disp ('dhondt法席位分配：')
1 g0 _/ A6 i0 ?: S8 B0 t
R1
; @- L y* ~) G
end
" Q3 Q9 a; I7 j2 j1 M& C
% e$ _! n. e, V- E, t. X
0 ]+ y) i' Z% {6 V( g" P
6 _& C9 Z& A8 W* F4 Y$ [$ `- {& P
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
';