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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
  S" y3 b7 u6 U8 B; {                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
# f8 h  y% Z6 [$ _# e下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
               clc
clear
nx=60;
ny=60;
%定义固定障碍物的状态
+ q" \4 T* X9 Nz=zeros(nx,ny);
peo=z;
( L- Q" S" I- R) }. e%peoNew=z;
wall=z;
wall(1,1:ny)=1;
wall(1:nx,1)=1;
5 x  L/ I/ A' \) x" @- k' a+ Twall(1:nx,ny)=1;
4 ^7 G" w, a6 twall(nx,1:ny/2-2)=1;
wall(nx,ny/2+2:ny)=1;
# T' W4 N! t; G- L% m) uwall(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;
2 c- S5 d0 @! Q, twall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
- W' ]7 F  o' L; E! S" A3 W%定义人的状态
for x=1:1:12
    for y=1:1:12
        x1=round(rand(1)*59+1);
        x2=round(rand(1)*59+1);
4 M* K* |" q2 Y9 D( U        if wall(x1,x2)~=1
           peo(x1,x2)=1;
: j- K' F% V% H/ L       end
2 ^, L$ f1 v: ^: b- A6 c    end
end
2 h! K! I7 P( j) @% d+ v  H: P%定义安全矩阵，指引人群运动
$ t* p) p2 T+ l7 k4 F4 T; Rsafe=zeros(60,60);
1 t8 k- a1 }0 G; Y3 G# Hsafe(nx,ny/2+1)=0;
safe(nx,ny/2+2)=0;
* g2 V8 q8 F8 L! ], nsafe(nx,ny/2)=0;
7 y; j2 r9 W( L% bsafe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
# k+ H" \( w) f/ sfor i=nx-1:-1:2
& |1 s/ ~3 f9 }& x) f4 u2 @' \    for j=ny/2+3:-1:2
9 T, a: g' S1 [7 ~+ A& ?3 m' ]; G        safe(i,ny/2)=safe(i+1,ny/2)+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;
* k3 \* S! {: O- `' I        safe(i,j)=safe(i,j+1)+1;
    end
end
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;
8 c: Q3 Y# l6 I4 B7 ~3 U+ z0 u4 i        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
- _9 O' T; x3 }. K  L% _        safe(i,j)=safe(i,j-1)+1;
5 V* L. w; J  t' K! }6 q2 N/ [    end
end
4 c" a- z" v- b/ M* D8 E2 a        
; ?1 l* V8 `& ~2 }% pimh=image(cat(3,z',peo',wall'));
! `# M& W3 y  R1 F- Q% L. W/ X/ m! V9 zset(imh,'erasemode', 'none')
axis equal
axis tight
* X8 \& E, X5 U8 S" u. m%peoNew=peo;
! t9 {# q0 @4 ^* f2 Y6 Axind=2:59;
yind=2:59;
0 ^2 U' p, r0 Pfor m=1:1000
' G( l4 u7 I  P0 G    if peo(xind,yind)==1
       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
' ?- g! F# q1 ]6 O9 R" Z3 @           peo(xind,yind-1)=1;
7 N. g& X1 F: F' |           peo(xind,yind)=0;
       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
" `! q: j6 T, y4 v$ t$ V           peo(xind,yind+1)=1;
, H$ ^: H% i) w0 h& b  |6 M! U           peo(xind,yind)=0;
/ U+ |( y4 O9 \" h# a- }9 k       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
3 x& ?% ~+ x  _9 r! [           peo(xind-1,yind)=1;
           peo(xind,yind)=0;
       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
- A" S0 t! U- I0 B! B           peo(xind+1,yind)=1;
- B1 \4 Z" i# @2 Z4 n7 v           peo(xind,yind)=0;
       end
* ?) s0 q8 l! p) T2 ?         peo=peoNew;
         set(imh, 'cdata', cat(3,z',peo',wall') )
         drawnow
         end
end

etzhu

实际上不用那么麻烦，
, P. J" s( A8 q; {% T# c% 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，友情帮顶。