QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3440|回复: 2
打印 上一主题 下一主题

[问题求助] 关于图论的广度优化搜索,语句while~isempty(g==1)正确吗?

[复制链接]
字体大小: 正常 放大

9

主题

7

听众

188

积分

升级  44%

  • TA的每日心情
    慵懒
    2013-8-23 15:15
  • 签到天数: 84 天

    [LV.6]常住居民II

    自我介绍
    我想参加明年的数模竞赛

    群组西安交大数学建模

    群组数学建摸协会

    群组英语科技论文写作实训

    跳转到指定楼层
    1#
    发表于 2013-3-2 17:01 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    此程序的14行,while~isempty(g==1),为什么行呢,什么意思?难道这个和while~isempty(find(g==1))是同一个意思?
    8 h. U: \. W) `  Y) T我个人理解,g==1除非g是空的,否则g==1都会返回一个逻辑矩阵才对啊,不会变成空的。
    ) S: z/ o7 w6 s2 Gfunction [w] = bfsf1(g)3 a& c3 g( G  d' l" |/ I
    %广度优先搜索
    5 K9 @1 u: l' L7 t  v  `2 u%g:邻接矩阵。w顶点标号% F6 |5 F0 O9 ?* h5 V+ j5 A1 F9 W4 o
    n=size(g,1);% `! P# y. o; c2 ?% G; `
    w=zeros(1,n);1 J: m; u2 c: S* ], W% \3 E
    l=0;' q/ h# P: h8 X8 R# b
    v=1;# ^# N* k  p0 s1 T9 c* I
    a1=find(g(v,:)==1);) T- F; @+ Q3 W$ n! V/ I2 D. W  T
    g(v,a1)=2;0 P1 E& g  t4 h$ Y  D  k" t3 f2 P
    g(a1,v)=2;
    " @/ N# _1 `8 j' N, Rw(a1)=l+1;/ l5 N. w3 S: a. W: E
    s1=union(a1,v);" V  t+ J/ ]* V% C& R
    l=l+1;
    + A/ Q) \4 n6 o5 w( Xwhile~isempty(g==1). R' c8 F  E, S  g0 c
       a1=find(g(s1,:)==1);
      N1 e2 L/ K7 }0 [, K    t=length(s1);
    ! k$ [! p# m; |5 c8 G    d=[];  W' X9 H# B8 ?( q8 T
        for i=1:length(a1)% G' r& B' m0 a4 q# a
            if a1(i)/t>floor(a1(i)/t)" Q1 @% Z& w' g) Y7 c1 q9 y4 E
                t2=floor(a1(i)/t)+1;9 `) l6 e3 {6 j1 S& x
            else5 h9 B1 }" x* ^. F, G6 }) t
                t2=floor(a1(i)/t);
    3 I$ U5 i8 ]$ G" t' ~        end' }& j4 ?5 z, D
            if isempty(intersect(d,t2)), J0 z# |' M# e! o5 P9 C9 T, `9 ^
                d=union(d,t2);; [0 W; x8 A- |
            end8 i2 W) C2 Q4 u9 v; K& S% L
        end
    4 c" I8 L% y% F% \  p- j0 m    d1=setdiff(d,s1);9 y/ u' r5 Q, Z' Z
        if isempty(d1)
    2 D; \+ ~9 n, L; {. V+ k9 V( F$ B        break;9 \& S% d  i2 ~, o$ J1 U0 y
        else
    7 O  O" f  X% F# x# k        w(d1)=l+1;
    4 }- a6 ^5 d! c% H2 c; W        g1=g(s1,:);
    7 g0 i( p" H' p& @8 j8 a: `# w        g1(a1)=2;7 }1 f) s. T% k: i/ q
            g(s1,:)=g1;
    3 w/ ]% a/ A; i: O* |  v- c        g(:,s1)=g1';
    4 o" I2 [, {+ [* j8 D, i  `        s1=union(s1,d1);
    2 a$ n3 a0 _+ N# @        l=l+1;
    7 f; r; M% _* n. {    end' G( R3 \1 |; t& d. O# e8 m
    end
    9 g! J; B' `* G" |$ j    1 t+ O/ U; k( \' \6 v5 g0 J
    end
    * H4 d7 @" H% s9 \# X5 {( c8 n! a
    ' t1 g4 I7 E  j) A7 w
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

    9

    主题

    7

    听众

    188

    积分

    升级  44%

  • TA的每日心情
    慵懒
    2013-8-23 15:15
  • 签到天数: 84 天

    [LV.6]常住居民II

    自我介绍
    我想参加明年的数模竞赛

    群组西安交大数学建模

    群组数学建摸协会

    群组英语科技论文写作实训

    帖子里面的笑脸其实是: 和) 不知怎么搞的打出来了就变成了笑脸

    点评

    madio  这种情况在右侧点击“禁用表情”就没有了。  详情 回复 发表于 2013-3-2 21:34
    回复

    使用道具 举报

    madio        

    3万

    主题

    1311

    听众

    5万

    积分

  • TA的每日心情
    奋斗
    2024-7-1 22:21
  • 签到天数: 2014 天

    [LV.Master]伴坛终老

    自我介绍
    数学中国站长

    社区QQ达人 邮箱绑定达人 优秀斑竹奖 发帖功臣 风雨历程奖 新人进步奖 最具活力勋章

    群组数学建模培训课堂1

    群组数学中国美赛辅助报名

    群组Matlab讨论组

    群组2013认证赛A题讨论群组

    群组2013认证赛C题讨论群组

    ganquanlife 发表于 2013-3-2 17:04 4 _3 g3 _/ P9 [$ D3 V( s" U
    帖子里面的笑脸其实是: 和) 不知怎么搞的打出来了就变成了笑脸
    . Z- W' H! E) E) K* `& |
    这种情况在右侧点击“禁用表情”就没有了。
    数学建模社会化
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2025-10-3 14:41 , Processed in 0.555123 second(s), 67 queries .

    回顶部