谁能真正解决这个问题，送20个矩阵币

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
" A& `* K6 N2 @9 J                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
, J8 D$ u# k7 B                      4。看所有人离开房间需要多少时间。
& i3 ~' n# Z  w# z下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
9 I6 D) y/ i  c; w$ }               clc
5 k* c4 M. a1 z4 K( _- dclear
" e/ ]. x" H4 b! A* P6 Hnx=60;
ny=60;
! w( {3 p( x- K8 v  z1 z8 I%定义固定障碍物的状态
z=zeros(nx,ny);
; `/ q6 {, W1 F% zpeo=z;
%peoNew=z;
wall=z;
wall(1,1:ny)=1;
wall(1:nx,1)=1;
6 e$ Z1 m" r& V! r) Mwall(1:nx,ny)=1;
2 I8 l: w7 p( a, [# H* w8 wwall(nx,1:ny/2-2)=1;
9 ^' ~/ r- v) F$ Mwall(nx,ny/2+2:ny)=1;
: [2 {# U% G- Dwall(nx/4:nx/4+1,ny/4:ny/4+1)=1;
wall(3*nx/4:3*nx/4+1,ny/4:ny/4+1)=1;
7 c5 ?) P4 X3 d; p, uwall(nx/4:nx/4+1,3*ny/4:3*ny/4+1)=1;
wall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
' w: l# o. G" s' C%定义人的状态
7 @$ I/ s' D6 b! Y/ O0 ]' Xfor x=1:1:12
    for y=1:1:12
        x1=round(rand(1)*59+1);
& o# f8 K- h( \) F6 T        x2=round(rand(1)*59+1);
        if wall(x1,x2)~=1
! F* E/ P) `8 G0 {           peo(x1,x2)=1;
       end
/ k4 y% m/ a# ^    end
end
! y8 o2 i! y  M%定义安全矩阵，指引人群运动
6 z0 j0 m3 N/ c# L" N0 asafe=zeros(60,60);
safe(nx,ny/2+1)=0;
safe(nx,ny/2+2)=0;
1 E9 S' D9 g( d, d% W1 ~- ?safe(nx,ny/2)=0;
  @2 L4 {' G! M( `safe(nx,ny/2-1)=0;
" B: h0 [9 C. ^' s! q% J* O3 Dsafe(nx,ny/2-2)=0;
. E6 \3 m* X5 N: R$ r- T1 {1 Gfor i=nx-1:-1:2
' K3 `8 h& p; O/ U/ n    for j=ny/2+3:-1:2
+ F" {9 n" U( M+ @5 `  S        safe(i,ny/2)=safe(i+1,ny/2)+1;
        safe(i,ny/2-1)=safe(i+1,ny/2-1)+1;
( @1 y% Q# Y# `4 M        safe(i,ny/2-2)=safe(i+1,ny/2-2)+1;
  m7 @  e- p* X/ c3 n        safe(i,j)=safe(i,j+1)+1;
% `# V- }/ U! H0 B2 o3 K1 k0 I0 B    end
3 Z3 d% i7 ~; K+ @5 ?  @end
for i=nx-1:-1:2
, u8 a# j1 m" J% m) ]) P' b    for j=ny/2+3:1:ny-1
        safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
% z7 l' O% T6 H7 b/ @. u1 \% _7 Q0 n        safe(i,j)=safe(i,j-1)+1;
    end
6 C" u8 l  e, w; Jend
6 U0 I) n. g( D" L- T/ [0 w        
* k! T9 ^1 H5 t6 I/ Fimh=image(cat(3,z',peo',wall'));
set(imh,'erasemode', 'none')
0 P9 R! @/ q9 |4 s) iaxis equal
axis tight
6 k4 l# \3 l' H9 L0 k%peoNew=peo;
! S9 `* N5 N! ^; C. cxind=2:59;
yind=2:59;
1 R. g3 ?* |# W; k' O* V. ]for m=1:1000
1 F& e+ X5 Q0 p$ [1 |8 e% J    if peo(xind,yind)==1
       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
* c: v4 {, h0 M$ x( W           peo(xind,yind-1)=1;
           peo(xind,yind)=0;
       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
           peo(xind,yind+1)=1;
           peo(xind,yind)=0;
       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
           peo(xind-1,yind)=1;
           peo(xind,yind)=0;
; `' t7 r$ `2 q; E  R       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
$ P7 k3 i( c! V$ i) h- L           peo(xind+1,yind)=1;
3 f" {5 M8 a" v           peo(xind,yind)=0;
# s: K4 x' a  w4 P" d       end
' e1 l+ B9 H# U9 E+ H         peo=peoNew;
         set(imh, 'cdata', cat(3,z',peo',wall') )
9 H9 x, D" f& {+ T         drawnow
8 I* R1 M6 |1 y8 J% A         end
end

etzhu

实际上不用那么麻烦，
8 o! l. i* Y  ~- d8 Y9 D5 u首先这个问题必然是满足 贪心 前提的，所以可以首先计算出每个人移动到门口的时间，然后排序，从小到大出门，考虑相等时的拥塞，就可以了。
2 E7 ?+ J: C- y* ]/ @" `; h/ L: Z当然，也有一点不足，这里没有考虑一个人挡在另一个人的路上的情况。

suxpert

up，友情帮顶。

sphereballball

现在是休息重要，别再动脑了

lxsoft110

2# etzhu
呵呵~兄弟。你还是没有说明我这个算法错误在哪里啊~

qsoft_cn

你调试一下可以发现，，if peo(xind,yind)==1 这个有问题，for循环里面的一次都没有执行。。。我正在想，应该有if find(peo(xind,yind)==1),，，厄，现在我该休息了。没时间再想了。元胞刚刚学习。。。。

wxl1988

我也要休息了~不看了~~好成绩·

cldndx

bless明天遇到顺手的题目！

victoriajh

注意休息啊 楼主

terrance_ho

up，友情帮顶。