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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
3 h* a8 A# ]# F! t我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
               clc
2 ]  n( y( N: \3 Uclear
" L' Y8 N% r: r' Snx=60;
/ y0 b1 A4 I, ^4 hny=60;
7 Z  U, P7 g0 j7 t  ~8 s% Y8 _$ c" l' t%定义固定障碍物的状态
" N1 X/ X* h% U' Oz=zeros(nx,ny);
7 `7 C. U! _: \, @+ @+ Ipeo=z;
%peoNew=z;
* a# R" I& v7 x) S  r0 }; vwall=z;
wall(1,1:ny)=1;
wall(1:nx,1)=1;
wall(1:nx,ny)=1;
; }- w  X" U1 k/ |  iwall(nx,1:ny/2-2)=1;
wall(nx,ny/2+2:ny)=1;
9 R" U4 G" ~! L( p7 x3 T4 wwall(nx/4:nx/4+1,ny/4:ny/4+1)=1;
8 S& c1 a- \9 M$ Bwall(3*nx/4:3*nx/4+1,ny/4:ny/4+1)=1;
& y+ M5 ~0 l9 fwall(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;
0 y- f" f% d5 @- z0 _5 _%定义人的状态
for x=1:1:12
# s9 a& C5 `1 O5 l    for y=1:1:12
        x1=round(rand(1)*59+1);
4 d) C% @; v! {2 q1 S4 j        x2=round(rand(1)*59+1);
: ?( b8 W  B/ g* C" A: f# d        if wall(x1,x2)~=1
7 n& P* T2 C+ {8 |0 P7 n           peo(x1,x2)=1;
& C- J4 o% K3 o: i       end
    end
/ Y' U3 E% m5 C  `7 z0 x6 Uend
  H, Z6 E& v+ v! T%定义安全矩阵，指引人群运动
4 c1 H5 [2 N9 e) ]$ x& a& }safe=zeros(60,60);
! s9 T. e! {2 l/ x! _% y$ l, Dsafe(nx,ny/2+1)=0;
8 H! i  G) m0 O9 K' H* osafe(nx,ny/2+2)=0;
safe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
. h3 O6 S7 F# k( isafe(nx,ny/2-2)=0;
0 L7 T# S8 ?/ gfor i=nx-1:-1:2
    for j=ny/2+3:-1:2
, q- @* G; w! ^0 d9 O        safe(i,ny/2)=safe(i+1,ny/2)+1;
. s4 u! X# v7 ]( k$ b/ f( x$ T        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
7 |9 V) O8 k; g+ Y7 rend
) N6 M# p# [: @7 a" zfor i=nx-1:-1:2
, t/ S% a+ t; J* Y  }/ T0 I    for j=ny/2+3:1:ny-1
; Q) w0 R& e) d" `& ~* g! e2 T        safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
% u' }1 ~. F: i7 t% {3 }/ y+ ^        safe(i,j)=safe(i,j-1)+1;
    end
end
        
imh=image(cat(3,z',peo',wall'));
2 M6 L% @& d* M) T# e& mset(imh,'erasemode', 'none')
/ ?+ u% E" U$ d( g0 faxis equal
- W" {* _( ?; R; k3 \axis tight
& [4 O8 G: G, `/ s1 ^8 h%peoNew=peo;
' Y# M, S0 R0 B9 N5 w( ?  gxind=2:59;
, [( P7 b$ _! Z  _- Jyind=2:59;
* Y. K  \- w$ E; l2 Jfor m=1:1000
    if peo(xind,yind)==1
       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
           peo(xind,yind-1)=1;
+ ~7 s! x8 Y$ c" Q           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
# k# Z; ?# q% r           peo(xind-1,yind)=1;
8 O; K5 z# ]6 K2 i$ |" K7 v/ V6 @           peo(xind,yind)=0;
       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
           peo(xind+1,yind)=1;
- `+ s1 z1 j1 f5 e/ J8 G           peo(xind,yind)=0;
       end
         peo=peoNew;
         set(imh, 'cdata', cat(3,z',peo',wall') )
         drawnow
         end
end

etzhu

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

suxpert

up，友情帮顶。

sphereballball

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

lxsoft110

2# etzhu
) \: @5 x2 M+ }/ u, U呵呵~兄弟。你还是没有说明我这个算法错误在哪里啊~

qsoft_cn

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

wxl1988

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

cldndx

bless明天遇到顺手的题目！

victoriajh

注意休息啊 楼主

terrance_ho

up，友情帮顶。