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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
; e" \& J& d- n7 ?" \- i                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
$ H- j7 D2 C9 I" d下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
               clc
+ `" l) n! N( s/ w% v! X+ Eclear
nx=60;
ny=60;
# F. z2 F9 u- F%定义固定障碍物的状态
6 I! c1 y/ r' G- ?z=zeros(nx,ny);
: Z9 A& j1 y6 ^  dpeo=z;
) c. {9 }: I8 I8 _%peoNew=z;
( G: }* }% f& W  ], x. |  Cwall=z;
wall(1,1:ny)=1;
wall(1:nx,1)=1;
wall(1:nx,ny)=1;
( z9 Y* s* T. I' D7 r! I! Owall(nx,1:ny/2-2)=1;
9 K* S7 X( n  k: e4 U. Dwall(nx,ny/2+2:ny)=1;
' c) d& p' U' x( c, m/ Zwall(nx/4:nx/4+1,ny/4:ny/4+1)=1;
3 k! y. F+ h3 p! v- v% E# S# ]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;
: ?4 m( h& K& Z6 k: m5 fwall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
1 f+ \5 v- f2 Y5 j%定义人的状态
9 S0 X: F9 l( [8 u4 y( q+ i9 bfor x=1:1:12
" I) W7 t2 A8 ~    for y=1:1:12
        x1=round(rand(1)*59+1);
        x2=round(rand(1)*59+1);
        if wall(x1,x2)~=1
           peo(x1,x2)=1;
       end
    end
1 m* Q, u: z" D: Hend
" R3 N* L( m2 ]' z$ t- d' P%定义安全矩阵，指引人群运动
( z. m+ p. M4 e$ O7 b: L( ~safe=zeros(60,60);
% Y& W1 D2 G+ ~! J  Osafe(nx,ny/2+1)=0;
( p' _: {8 f. isafe(nx,ny/2+2)=0;
: P, M, n  C0 I1 [# k: ]9 b$ z% @4 Jsafe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
for i=nx-1:-1:2
( B1 r( @; N( c: z! T7 M    for j=ny/2+3:-1:2
7 r4 a+ p2 q. ]: ]4 u        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;
        safe(i,j)=safe(i,j+1)+1;
* C; @6 k* r9 Y, G0 [2 W2 L2 l, K    end
- t. B. l8 R: {3 P! Hend
% V5 Q, V! x8 B9 Kfor i=nx-1:-1:2
8 @' `; I1 c; D    for j=ny/2+3:1:ny-1
* d0 [8 T; v/ N& H0 J! \( B$ [        safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
0 C, ~( Y; {- w' P& F        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
        safe(i,j)=safe(i,j-1)+1;
% h  r* I9 W1 }3 U$ x: W    end
( A& k7 f) X8 y) I4 Nend
. S  r+ j' I% u+ V' J' i        
imh=image(cat(3,z',peo',wall'));
3 U" j) w: s; Y* g7 Kset(imh,'erasemode', 'none')
" y, E# y! w5 y9 `) j3 v: w1 Z( jaxis equal
axis tight
- Z+ V- c1 ~4 d, N6 g%peoNew=peo;
. {/ P( e6 C% U2 G1 D( ]# Oxind=2:59;
yind=2:59;
. u3 U1 G2 v" I3 a9 C0 r$ O; Cfor m=1:1000
    if peo(xind,yind)==1
- L3 E  o7 ]; w8 [       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
- f' \* F: F/ Z; p, t: _           peo(xind,yind-1)=1;
$ n7 u' o. f' F           peo(xind,yind)=0;
       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
/ _' o" \* f8 B* }           peo(xind,yind+1)=1;
           peo(xind,yind)=0;
       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
           peo(xind-1,yind)=1;
5 z) s0 C; `7 Y. K6 M           peo(xind,yind)=0;
       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
           peo(xind+1,yind)=1;
# _5 r, |* k2 q0 G$ Z           peo(xind,yind)=0;
( k7 z' F8 N8 S8 G4 [5 o' [       end
! R! G9 V. ^6 s7 N  h4 ]  m& L8 d         peo=peoNew;
! _: p6 e% p. t% Y/ \* M3 l; c         set(imh, 'cdata', cat(3,z',peo',wall') )
, S* U6 t4 h  ~' [5 D  t         drawnow
         end
end

etzhu

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

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