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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
# x. R2 @# R% j1 A                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
8 c8 D( `7 s) C* G                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
' p6 `7 E7 K) A下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
+ Z) c; ]1 E/ ]4 G+ \               clc
clear
nx=60;
ny=60;
' f5 t1 c7 [8 Z%定义固定障碍物的状态
z=zeros(nx,ny);
, u% D2 N2 {' V# r) a3 a  opeo=z;
%peoNew=z;
wall=z;
5 j6 [4 v1 c- ]  O/ J! I4 C* ]wall(1,1:ny)=1;
8 W$ k: N- B1 }9 @wall(1:nx,1)=1;
wall(1:nx,ny)=1;
wall(nx,1:ny/2-2)=1;
wall(nx,ny/2+2:ny)=1;
wall(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;
wall(nx/4:nx/4+1,3*ny/4:3*ny/4+1)=1;
$ \6 j2 g- K1 N2 y8 L# K/ Ywall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
0 q& ]1 X' n9 K$ }( a2 z%定义人的状态
for x=1:1:12
3 G/ j6 }% X+ q* s3 P/ g    for y=1:1:12
; J0 y2 i' C( }% N- S4 o        x1=round(rand(1)*59+1);
        x2=round(rand(1)*59+1);
        if wall(x1,x2)~=1
# ?& Y1 U6 ^5 _, X( e: L& L' u) H           peo(x1,x2)=1;
       end
    end
+ d& c1 B5 S  t) X  d4 }( ^+ L! aend
%定义安全矩阵，指引人群运动
safe=zeros(60,60);
safe(nx,ny/2+1)=0;
safe(nx,ny/2+2)=0;
safe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
4 u6 c* K* H& }& v; M2 c! pfor i=nx-1:-1:2
    for j=ny/2+3:-1:2
2 ^- M5 E9 R3 Y) P& ~        safe(i,ny/2)=safe(i+1,ny/2)+1;
        safe(i,ny/2-1)=safe(i+1,ny/2-1)+1;
" \$ E3 C+ B. n  C: n        safe(i,ny/2-2)=safe(i+1,ny/2-2)+1;
        safe(i,j)=safe(i,j+1)+1;
0 v: b# z3 ^; J; z    end
: ?+ z, c, R0 E, ]0 cend
for i=nx-1:-1:2
    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;
# s& T: S4 A6 V( R1 ^" E9 _        safe(i,j)=safe(i,j-1)+1;
# X% @2 @* _2 P* E1 c. P8 R    end
8 \8 q2 z* S9 cend
        
imh=image(cat(3,z',peo',wall'));
  ?8 M( v4 A$ w- x' O$ X+ F7 lset(imh,'erasemode', 'none')
' l; d8 G; s5 |. P/ E: zaxis equal
axis tight
4 O  E& F" D0 O%peoNew=peo;
& p4 I; t# z# k0 `; l: z5 i1 m& \5 Bxind=2:59;
yind=2:59;
' {2 n* _: K( a4 y. a) afor m=1:1000
8 L* l: ~9 `( ]4 D    if peo(xind,yind)==1
# c/ z+ Y- w  G+ O, w5 r. E       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
9 H- [4 ?3 l% x2 K  k           peo(xind,yind-1)=1;
           peo(xind,yind)=0;
% }, c- k( P+ m3 d( j       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
' W  @" m7 c0 V           peo(xind,yind+1)=1;
           peo(xind,yind)=0;
       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
& j) _: [# w% ^1 L- X           peo(xind-1,yind)=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;
1 _( |* E2 ^' h% e       end
; _. Q% V, v4 C1 F         peo=peoNew;
$ H5 ~) e- l& A, h4 D& A! |         set(imh, 'cdata', cat(3,z',peo',wall') )
         drawnow
7 T3 X- ?; d3 K  ~         end
end

etzhu

实际上不用那么麻烦，
首先这个问题必然是满足 贪心 前提的，所以可以首先计算出每个人移动到门口的时间，然后排序，从小到大出门，考虑相等时的拥塞，就可以了。
当然，也有一点不足，这里没有考虑一个人挡在另一个人的路上的情况。

suxpert

up，友情帮顶。

sphereballball

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

lxsoft110

2# etzhu
  U- }$ n- L/ v5 G9 B3 T8 i6 S呵呵~兄弟。你还是没有说明我这个算法错误在哪里啊~

qsoft_cn

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

wxl1988

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

cldndx

bless明天遇到顺手的题目！

victoriajh

注意休息啊 楼主

terrance_ho

up，友情帮顶。