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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
0 I; G- g1 e6 W4 z6 l3 Q1 u4 S                      4。看所有人离开房间需要多少时间。
6 O: c1 i+ b, Y  M1 w" t下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
               clc
1 @( {6 s" z# U; y( ]3 [7 h6 {* ]clear
nx=60;
ny=60;
/ s/ h4 F! M) X%定义固定障碍物的状态
: @7 }  o5 P  N, Wz=zeros(nx,ny);
peo=z;
%peoNew=z;
9 t# p2 m4 ?( Jwall=z;
wall(1,1:ny)=1;
" {% M2 e  d2 |: f# v6 \& Qwall(1:nx,1)=1;
+ w  `3 Y# P) A% _! cwall(1:nx,ny)=1;
wall(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;
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;
wall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
7 Q- X; A/ W2 x  e%定义人的状态
for x=1:1:12
    for y=1:1:12
5 c* L6 o/ O/ K        x1=round(rand(1)*59+1);
) S2 [- ]( S7 c9 z- K        x2=round(rand(1)*59+1);
        if wall(x1,x2)~=1
9 o) T: [& c) X3 z% ~- w           peo(x1,x2)=1;
0 B( ]+ A6 `9 G. g       end
; u/ u  O3 T  o2 o- j    end
end
# g* L& G! V' B( g9 f%定义安全矩阵，指引人群运动
safe=zeros(60,60);
- L* H) D# b9 ^+ msafe(nx,ny/2+1)=0;
9 X% k+ K% g  a) c9 \. Qsafe(nx,ny/2+2)=0;
safe(nx,ny/2)=0;
" y4 `" E1 a' U1 w6 g/ u. K) m- Jsafe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
for i=nx-1:-1:2
4 ^) s' t+ ~  U1 I3 K3 O) @    for j=ny/2+3:-1:2
        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;
) ?0 X# j! F, s0 C/ n3 h2 l" I% E7 y        safe(i,j)=safe(i,j+1)+1;
    end
+ T9 V+ o9 F, [# r$ send
* J( @  L" G2 J. t  e! O) J: ufor i=nx-1:-1:2
    for j=ny/2+3:1:ny-1
; h5 Z( S$ w9 K9 A0 Q1 Q7 f        safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
; ~  \9 Q: h9 e& b9 @2 o        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
( f# j7 O. @$ X, u        safe(i,j)=safe(i,j-1)+1;
    end
0 G+ d, `7 b6 X( A( ]end
* z- o/ W/ W5 v, W" Y        
imh=image(cat(3,z',peo',wall'));
set(imh,'erasemode', 'none')
axis equal
& }$ N0 l$ U7 A4 b+ g" Uaxis tight
%peoNew=peo;
5 n% q* ^+ R* I/ p- P% D3 X- H- U6 Yxind=2:59;
/ I6 D: c; }% L' Cyind=2:59;
for m=1:1000
9 ^5 f1 |+ ]3 Z3 y9 M/ K' @% ]    if peo(xind,yind)==1
2 [) q7 s7 b9 G( w" ]       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
1 u7 x9 d, o. B* l           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;
& x  f2 m8 j( `5 `& E5 C       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
           peo(xind-1,yind)=1;
; e; f- i) s0 L! U9 ]# h; @  l           peo(xind,yind)=0;
9 }' }1 S; @- m       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
           peo(xind+1,yind)=1;
           peo(xind,yind)=0;
( `2 v6 w) U. `2 D       end
         peo=peoNew;
1 W' U8 k4 S0 Z( U2 V         set(imh, 'cdata', cat(3,z',peo',wall') )
# |9 l8 \/ m) |5 f' g         drawnow
         end
end

etzhu

实际上不用那么麻烦，
  h7 B- U) o% j( G& v# F5 w2 k/ b首先这个问题必然是满足 贪心 前提的，所以可以首先计算出每个人移动到门口的时间，然后排序，从小到大出门，考虑相等时的拥塞，就可以了。
当然，也有一点不足，这里没有考虑一个人挡在另一个人的路上的情况。

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