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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
$ e9 `$ _# t/ O0 G  ^我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
: n4 x% V  i! W& X8 c% y( T                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
  N& Z: q' H+ B8 S6 s6 r* [                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
8 e. ^6 u& I1 b, `                      4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
$ T! L! [+ m8 Z8 V- k) ^               clc
clear
7 D+ s! d  c# N5 Lnx=60;
! w" i- X; p! p) B7 Nny=60;
%定义固定障碍物的状态
+ i0 o# b; U5 D4 ]! Kz=zeros(nx,ny);
peo=z;
%peoNew=z;
wall=z;
wall(1,1:ny)=1;
" c8 F" [* M3 d8 k1 s7 c$ Fwall(1:nx,1)=1;
8 x: l. p& |7 s$ _9 C- swall(1:nx,ny)=1;
wall(nx,1:ny/2-2)=1;
. Z/ U4 P2 L: Owall(nx,ny/2+2:ny)=1;
wall(nx/4:nx/4+1,ny/4:ny/4+1)=1;
5 J" t( @1 v' c( V- Q$ uwall(3*nx/4:3*nx/4+1,ny/4:ny/4+1)=1;
1 d* i! t4 m8 O$ e; Vwall(nx/4:nx/4+1,3*ny/4:3*ny/4+1)=1;
) F: y1 {: T1 I: J9 kwall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
. m6 @$ ?- q& M4 N%定义人的状态
for x=1:1:12
9 I' g+ @: B% O0 x! X' k' c    for y=1:1:12
        x1=round(rand(1)*59+1);
- |1 Z* O/ l1 W! q' V$ l, g        x2=round(rand(1)*59+1);
$ L3 u8 _2 }3 h. s4 L0 Y        if wall(x1,x2)~=1
+ P( P, z" k1 G1 a           peo(x1,x2)=1;
       end
6 Y% R4 R3 b: H    end
end
%定义安全矩阵，指引人群运动
/ E5 i" I: b8 ~; Ssafe=zeros(60,60);
safe(nx,ny/2+1)=0;
safe(nx,ny/2+2)=0;
9 {# o4 w6 @+ D0 E; C: n7 Wsafe(nx,ny/2)=0;
( u2 O/ z& g$ _safe(nx,ny/2-1)=0;
; O2 |5 ]# n8 @- X: L  ssafe(nx,ny/2-2)=0;
for i=nx-1:-1:2
    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;
        safe(i,j)=safe(i,j+1)+1;
, n3 P. p- L" }( j    end
% V, n, [( v# Bend
for i=nx-1:-1:2
& ]$ J/ i2 e9 T' {, c    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;
& S$ \. o7 P$ q! o) e    end
- W& P4 [5 |9 }7 ]+ x9 Uend
        
imh=image(cat(3,z',peo',wall'));
% g1 f5 r/ i5 R9 y, \4 `9 ?7 l# {) bset(imh,'erasemode', 'none')
! Z3 u2 f* }8 D; {; [axis equal
axis tight
%peoNew=peo;
xind=2:59;
( m) ?# E% O$ m  y# `yind=2:59;
for m=1:1000
    if peo(xind,yind)==1
3 A. s/ U7 ^, N' j) w; M- H4 |       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
           peo(xind,yind-1)=1;
' `$ Z0 ]7 M$ `; I5 L: b6 s: m           peo(xind,yind)=0;
       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
, \* |- B$ P$ ~7 T7 U           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 n+ p8 J( U6 Q           peo(xind,yind)=0;
' [2 S4 A, z) X2 l* q       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
           peo(xind+1,yind)=1;
           peo(xind,yind)=0;
       end
+ ]1 s  m/ y, u8 ]6 n" m  E         peo=peoNew;
7 Q. M: p( I9 ~8 m, ~         set(imh, 'cdata', cat(3,z',peo',wall') )
- s" P9 V: ^! ~6 F& `& J- A         drawnow
* v; t; m" h2 C9 ]. J8 u7 R         end
- w& l5 L) Z, P* Dend

etzhu

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

suxpert

up，友情帮顶。

sphereballball

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

lxsoft110

2# etzhu
6 F: |- t2 A4 p呵呵~兄弟。你还是没有说明我这个算法错误在哪里啊~

qsoft_cn

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

wxl1988

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

cldndx

bless明天遇到顺手的题目！

victoriajh

注意休息啊 楼主

terrance_ho

up，友情帮顶。