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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
! Q, r4 u4 r5 e) u: {/ X                      4。看所有人离开房间需要多少时间。
1 j5 p2 D" ?; ?. j* d下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
: d1 V6 H9 V* B& C1 V/ a               clc
) O+ c) L0 [  E1 p  \clear
7 A' g1 U4 g) k4 _nx=60;
2 k" }* Q7 J6 E, Bny=60;
%定义固定障碍物的状态
& y: r, W6 G- ~2 k4 h  h8 Tz=zeros(nx,ny);
& [, N+ O9 N7 q6 o7 ]5 ^6 _peo=z;
+ Q; q- H. s$ a/ ~4 N4 z: p6 j8 ]%peoNew=z;
& I( V1 C, m6 i  T8 awall=z;
wall(1,1:ny)=1;
wall(1:nx,1)=1;
! w0 U  M% e( i* d1 Ywall(1:nx,ny)=1;
wall(nx,1:ny/2-2)=1;
wall(nx,ny/2+2:ny)=1;
+ D$ M/ B, j; ~1 s. lwall(nx/4:nx/4+1,ny/4:ny/4+1)=1;
# Q! \0 O& k! ^' Lwall(3*nx/4:3*nx/4+1,ny/4:ny/4+1)=1;
- `1 t3 R: `) }. Awall(nx/4:nx/4+1,3*ny/4:3*ny/4+1)=1;
/ Y  N- Y7 f' U. F( Qwall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
%定义人的状态
+ N5 }1 O# f7 u+ f% Nfor x=1:1:12
1 V4 U6 O# N# V    for y=1:1:12
+ e/ e8 K2 p$ d6 e' i% L        x1=round(rand(1)*59+1);
- @$ n1 _7 y) f! F        x2=round(rand(1)*59+1);
) s* j' ~' H- j2 L0 W( h6 O        if wall(x1,x2)~=1
% t  [) S; O0 i9 D  l  K7 c( N           peo(x1,x2)=1;
       end
* T' V" g5 S% {" x3 d/ F7 N' E    end
& ~* p7 W, p' I5 Kend
%定义安全矩阵，指引人群运动
0 ^7 \& V4 j- T& ysafe=zeros(60,60);
safe(nx,ny/2+1)=0;
: S% Y$ D0 O4 e1 `safe(nx,ny/2+2)=0;
safe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
for i=nx-1:-1:2
    for j=ny/2+3:-1:2
        safe(i,ny/2)=safe(i+1,ny/2)+1;
        safe(i,ny/2-1)=safe(i+1,ny/2-1)+1;
- \/ A9 o' E! s0 ~        safe(i,ny/2-2)=safe(i+1,ny/2-2)+1;
6 v/ _$ w" H" l1 z        safe(i,j)=safe(i,j+1)+1;
    end
end
8 b7 @3 x* ]/ [3 N) h1 n" Kfor i=nx-1:-1:2
    for j=ny/2+3:1:ny-1
1 u! W, F8 w0 S) t5 y8 _        safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
        safe(i,j)=safe(i,j-1)+1;
3 q+ x! B4 r) P    end
end
        
imh=image(cat(3,z',peo',wall'));
set(imh,'erasemode', 'none')
axis equal
$ s' x# w  l5 P5 ~5 C! ]axis tight
%peoNew=peo;
xind=2:59;
yind=2:59;
for m=1:1000
1 m% c4 [8 z3 L8 ]! P$ V' l    if peo(xind,yind)==1
$ h, A$ B2 h8 R; m4 u. q5 k) ^       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
           peo(xind,yind-1)=1;
           peo(xind,yind)=0;
       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
. s! w1 r8 ~/ p& }/ ^           peo(xind,yind+1)=1;
           peo(xind,yind)=0;
3 z( w7 B; ?6 C7 L5 d( f       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
           peo(xind-1,yind)=1;
           peo(xind,yind)=0;
+ z3 o3 w8 G. I' }$ B& l  J8 M% U       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
, U/ Q8 U5 c" {           peo(xind+1,yind)=1;
           peo(xind,yind)=0;
7 @! Z# F$ J- Q' j0 t       end
2 o% L, ~5 B6 [0 Q, _3 n9 x         peo=peoNew;
         set(imh, 'cdata', cat(3,z',peo',wall') )
, D5 b. \9 c* Q! S! D1 b  C7 o- a         drawnow
( K# [% \1 {! F+ w         end
" y7 i8 u8 F1 c9 `" g9 uend

etzhu

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

suxpert

up，友情帮顶。

sphereballball

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

lxsoft110

2# etzhu
/ O( t) |4 Z- W2 E- f; i呵呵~兄弟。你还是没有说明我这个算法错误在哪里啊~

qsoft_cn

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

wxl1988

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

cldndx

bless明天遇到顺手的题目！

victoriajh

注意休息啊 楼主

terrance_ho

up，友情帮顶。