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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
- z, ^& c6 e2 a2 W; Y1 ]5 ^; c               clc
2 Y5 I9 z/ k: u- r# n6 e% Uclear
, l6 ]  K" U. i4 S+ _! T0 Gnx=60;
ny=60;
$ u2 a4 u7 l0 b0 X& f; f: o%定义固定障碍物的状态
9 O& V7 I' |. k- bz=zeros(nx,ny);
peo=z;
% n5 _  z! g3 ~5 I/ Y9 ~' S# R%peoNew=z;
wall=z;
$ P/ X$ t+ O' v+ K0 p! R9 Ywall(1,1:ny)=1;
7 \2 \* E, N, |( l: c3 A' Twall(1:nx,1)=1;
( B" o4 Q  y4 l7 t$ {' Dwall(1:nx,ny)=1;
/ R$ t, `5 q' G6 E1 o3 vwall(nx,1:ny/2-2)=1;
wall(nx,ny/2+2:ny)=1;
0 E! _/ O8 l7 M! b' k* |0 pwall(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;
$ y5 k+ Q' W, X) f! }wall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
6 g/ y' D* \2 J% f* e: A%定义人的状态
/ F0 H' g( r8 i# X$ ]% D6 N6 }for x=1:1:12
) S# M' V7 J, C2 R1 B    for y=1:1:12
        x1=round(rand(1)*59+1);
- h) E! h/ [& X: n! v        x2=round(rand(1)*59+1);
% i1 w4 g! `6 T- F, S5 l  o        if wall(x1,x2)~=1
           peo(x1,x2)=1;
       end
: G& f* B6 |7 B1 D, ^! d    end
7 W) Z* C& U$ }. U/ g  T; Uend
0 |; }8 G  C5 z; h. R4 p: Q% W0 u%定义安全矩阵，指引人群运动
9 _0 U( E: G* x* w' ]% s# |safe=zeros(60,60);
8 F# D! d  u0 I3 l8 k3 rsafe(nx,ny/2+1)=0;
% R- ^/ Q/ @- M/ w+ g+ g& b* G3 hsafe(nx,ny/2+2)=0;
safe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
8 i' b7 d: E0 Z: ]+ J% Msafe(nx,ny/2-2)=0;
$ f7 K$ K6 [/ h  ?0 P1 u; Z) C$ W) `for i=nx-1:-1:2
7 E5 y( [4 F  I( ^( Z3 D, E    for j=ny/2+3:-1:2
1 `3 m5 ?+ t2 m; i. }8 @        safe(i,ny/2)=safe(i+1,ny/2)+1;
+ e. \) a. ?' p5 L        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;
    end
end
  x9 N7 u1 y* t" K. O% Z  ], ~+ L' _for i=nx-1:-1:2
- Z, k& I, a! T+ J8 v) z: F" C: J    for j=ny/2+3:1:ny-1
8 z9 H& e: A: ]) v* ~        safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
  y2 ?) J: `* F3 S1 Y9 K" F$ o        safe(i,j)=safe(i,j-1)+1;
, A/ O, k8 z% j- m) d9 v    end
end
4 W0 @/ q9 G" Z3 V* J        
imh=image(cat(3,z',peo',wall'));
set(imh,'erasemode', 'none')
axis equal
axis tight
%peoNew=peo;
xind=2:59;
yind=2:59;
for m=1:1000
) ]) P9 q: x8 K. H# W0 C    if peo(xind,yind)==1
       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
) g1 N% \9 m& A. |' p           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
; O# {* i4 f7 B5 q; P           peo(xind-1,yind)=1;
           peo(xind,yind)=0;
       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
( x0 q4 i3 k( y: q. u           peo(xind+1,yind)=1;
" i" M. U1 i2 ]. v5 r           peo(xind,yind)=0;
/ \2 v% m1 ~8 p5 O: S; r; R       end
         peo=peoNew;
: Y9 I1 P9 R/ v+ g7 W         set(imh, 'cdata', cat(3,z',peo',wall') )
( g4 U7 ?1 L+ _1 {# |         drawnow
% m' n  c+ f  a1 }$ D1 F         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，友情帮顶。