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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
/ E& r6 r& A+ D! P( ^我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
& r* O/ p. {& Y" u( K# M! F                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
               clc
clear
nx=60;
ny=60;
1 g. n) j* W, p%定义固定障碍物的状态
z=zeros(nx,ny);
peo=z;
  Z, t  _2 K. x1 k%peoNew=z;
1 g! a2 a% ^& y2 l, Ewall=z;
wall(1,1:ny)=1;
/ k$ t: b2 o+ r' |& r% Vwall(1:nx,1)=1;
wall(1:nx,ny)=1;
" i( N4 q' l8 d( @6 J6 C, Z  M$ Pwall(nx,1:ny/2-2)=1;
wall(nx,ny/2+2:ny)=1;
* x$ s. ]" ~( q( Hwall(nx/4:nx/4+1,ny/4:ny/4+1)=1;
( E! R' e$ @. y/ Nwall(3*nx/4:3*nx/4+1,ny/4:ny/4+1)=1;
" B$ M& ^4 E4 L3 D4 j% Y% j" T) 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;
%定义人的状态
for x=1:1:12
. C- D8 u: w: x5 a* h    for y=1:1:12
        x1=round(rand(1)*59+1);
        x2=round(rand(1)*59+1);
        if wall(x1,x2)~=1
6 M$ H, C2 D2 W9 w6 }$ h: K           peo(x1,x2)=1;
       end
    end
end
%定义安全矩阵，指引人群运动
* ~  [8 s  i* b$ v9 dsafe=zeros(60,60);
" i6 [" a8 t; vsafe(nx,ny/2+1)=0;
safe(nx,ny/2+2)=0;
3 E8 E% s4 D7 d2 ~& \# ~2 \safe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
# f/ x: _% L4 l( nsafe(nx,ny/2-2)=0;
for i=nx-1:-1:2
! X( g6 l; Z- {! ~& f    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;
- z, p4 U, V' d, D; h        safe(i,ny/2-2)=safe(i+1,ny/2-2)+1;
        safe(i,j)=safe(i,j+1)+1;
    end
end
; X8 I! l1 z. [! Z) V$ h2 `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;
- Y$ n1 C! F. G) `        safe(i,j)=safe(i,j-1)+1;
+ K( n) a! B) {4 T    end
end
' O8 R' k. j  K# `, V3 {        
imh=image(cat(3,z',peo',wall'));
set(imh,'erasemode', 'none')
' @" X) ~0 w0 X; K" Q: F* a) Z& haxis equal
axis tight
%peoNew=peo;
. [7 h: F& z3 a2 H3 R# @! h! xxind=2:59;
7 T/ E# u5 _: xyind=2:59;
for m=1:1000
# u6 @3 v) V) |7 v) {    if peo(xind,yind)==1
5 H% I8 k* g, L       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
5 y, |' Z" B# v$ U* o/ R7 {           peo(xind,yind-1)=1;
, k4 s% e# d' K: C( z! c; V- t           peo(xind,yind)=0;
9 p# L  K6 ]) \' Z1 ~: h/ D       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
           peo(xind,yind+1)=1;
           peo(xind,yind)=0;
& W' `) ^4 u; A/ h* W7 A9 A7 ~, S# p       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
, W9 m- G. ~4 u' u           peo(xind-1,yind)=1;
7 c0 E8 x& M  i3 K+ u           peo(xind,yind)=0;
       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
+ }% |: W! u" j# v* W/ v! |           peo(xind+1,yind)=1;
, c0 Z0 H& @& p( z5 A% g8 ?0 i           peo(xind,yind)=0;
       end
         peo=peoNew;
* m2 e  K# H( o! n         set(imh, 'cdata', cat(3,z',peo',wall') )
         drawnow
         end
end

etzhu

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

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，友情帮顶。