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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
  K$ G/ [/ O/ A! K# v1 D: n                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
. a. H( |. s4 C/ m9 i- d* E                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
               clc
clear
nx=60;
. ~% h' O9 s2 u6 i4 o; Fny=60;
%定义固定障碍物的状态
! m0 s6 M4 L$ Mz=zeros(nx,ny);
peo=z;
%peoNew=z;
. u% L6 h  s+ o; ^% Y: i7 D9 u5 M; owall=z;
wall(1,1:ny)=1;
) K+ a9 [  Y' h' I+ _1 B6 Twall(1:nx,1)=1;
7 y3 ~6 U/ X4 i$ O2 S9 Ywall(1:nx,ny)=1;
+ V# c( e* C$ `0 jwall(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;
; V, I" h4 H* _- e% @wall(3*nx/4:3*nx/4+1,ny/4:ny/4+1)=1;
  e0 G- C: m( G# S' w; Z# mwall(nx/4:nx/4+1,3*ny/4:3*ny/4+1)=1;
% N8 N; @: w- G. d  hwall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
%定义人的状态
% K- ^! A  m: r" [, g3 qfor x=1:1:12
7 K( h* B$ S; \! V    for y=1:1:12
        x1=round(rand(1)*59+1);
: N3 c! \3 E7 k( H7 b, A5 ~# r        x2=round(rand(1)*59+1);
        if wall(x1,x2)~=1
           peo(x1,x2)=1;
% A  p- m3 M  |6 |4 Z- \, x5 \       end
    end
end
%定义安全矩阵，指引人群运动
safe=zeros(60,60);
& E, q9 z( p: Z1 Ssafe(nx,ny/2+1)=0;
safe(nx,ny/2+2)=0;
! I  O+ Z7 O* l4 n" X' ^! x8 U: wsafe(nx,ny/2)=0;
& @9 k9 ^1 w9 Z; s. u% P/ V& \safe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
1 v* T) C, A1 z* k2 A5 I% Y0 ~8 ifor i=nx-1:-1:2
    for j=ny/2+3:-1:2
* m% P7 K' v# v        safe(i,ny/2)=safe(i+1,ny/2)+1;
6 }5 A  Q5 N. I5 v7 S6 \4 `        safe(i,ny/2-1)=safe(i+1,ny/2-1)+1;
7 ]' r3 A, d$ Z0 j! i; K& w. d( D( }1 P5 y        safe(i,ny/2-2)=safe(i+1,ny/2-2)+1;
        safe(i,j)=safe(i,j+1)+1;
3 G9 C% J+ i8 n& X7 K    end
9 C" {1 K, i* W6 h$ ~6 M# u. x- Vend
for i=nx-1:-1:2
% ?" v2 j4 S# j6 e2 P    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;
        safe(i,j)=safe(i,j-1)+1;
    end
end
% i1 J3 M, a  G0 U9 D        
5 m& f' s, J& v2 Aimh=image(cat(3,z',peo',wall'));
set(imh,'erasemode', 'none')
. k  K4 _1 B7 J7 Oaxis equal
axis tight
%peoNew=peo;
2 n9 y  f3 V5 Rxind=2:59;
yind=2:59;
: W# I- ]: s9 d. F" m2 j1 I" w- j( _for m=1:1000
    if peo(xind,yind)==1
2 u) i% T% }4 n! Y4 D3 O5 ?$ Y       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
           peo(xind,yind-1)=1;
, _% o, D" e7 g' d+ S! x. z6 N           peo(xind,yind)=0;
( X3 b8 D) L7 u- t; t( r: h       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
$ d2 [% z& }* b# H0 {! [           peo(xind,yind+1)=1;
* Q7 c1 I7 e% S8 t7 U           peo(xind,yind)=0;
       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
           peo(xind-1,yind)=1;
- M6 f) b7 j+ i2 v           peo(xind,yind)=0;
, a$ ^" y# T0 z1 [8 y0 d% Z! ~       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
0 c7 ?/ p7 f1 d1 i% U' L" H. Y           peo(xind+1,yind)=1;
: `% ]* ]; \2 G$ ~4 l- `' ?3 w           peo(xind,yind)=0;
( c* {, x2 s" k       end
         peo=peoNew;
         set(imh, 'cdata', cat(3,z',peo',wall') )
- v. }9 ~! b% }0 K6 ^' D         drawnow
' b5 w$ K/ ~$ T: ]/ q0 i% \9 J# u         end
end

etzhu

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

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