QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3611|回复: 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 a# W. |/ Y& s: a
    我个人理解,g==1除非g是空的,否则g==1都会返回一个逻辑矩阵才对啊,不会变成空的。
    7 F' x7 @% P) X+ v7 mfunction [w] = bfsf1(g)
    ' M: h/ T$ Q' [8 \) r. l/ V7 ]: I%广度优先搜索
    6 w, i8 @9 r  T/ S' |0 V8 B%g:邻接矩阵。w顶点标号
    6 }  A0 j1 G. {% Z1 i7 qn=size(g,1);
    9 Z/ K4 O- K* k( F4 t  D. B3 W3 Dw=zeros(1,n);! v( Z, q; t6 w% R* X1 o0 S
    l=0;
    . Q9 \( g7 a$ n3 K2 l" N) pv=1;
    2 w$ D" W5 J& Y" P3 y8 e& q7 ma1=find(g(v,:)==1);8 x4 E0 m9 f$ Q  Y5 g4 ^2 \& ?
    g(v,a1)=2;
    1 H8 \  K% L/ l9 Kg(a1,v)=2;* ^$ m' |: w9 c. a
    w(a1)=l+1;+ V, O3 F$ s: Y& e; D
    s1=union(a1,v);8 R; C# t8 M. f
    l=l+1;
    1 H! d! B0 ~' D6 o' y8 f2 i) fwhile~isempty(g==1)' s6 a. a% P* H
       a1=find(g(s1,:)==1);* g$ L; U3 s" U, J0 O# p
        t=length(s1);
    # s( l0 E7 q4 }) ~4 T0 g8 k    d=[];
    1 p# _4 y% v, {) c! Y& x) e. Q    for i=1:length(a1)( M' _  o" j: A
            if a1(i)/t>floor(a1(i)/t)
    4 c" s- [3 \* Z; k. a. S8 o            t2=floor(a1(i)/t)+1;
    , s0 t6 f# L7 D5 B. K0 q        else
    : k. I! ~* U( a3 L( e3 p) M            t2=floor(a1(i)/t);0 k# A" I% Y# ]
            end8 X! o8 H0 O2 [8 ~" w: @
            if isempty(intersect(d,t2))
    2 r5 P2 o' u& ~, @! M2 n% t: m3 b            d=union(d,t2);
    ! A1 q( m% u, h: m+ ]) r: R  ^        end
    0 k8 J' U! H! G2 J. A3 ?9 C    end6 M* _/ Z# H3 o6 s2 W8 J, @# Y
        d1=setdiff(d,s1);
    . P- y- @* _1 v# i# R% S0 D: G( N2 {    if isempty(d1)3 |: s) X8 p  k3 j/ i6 z# C! O
            break;
      f4 A: ~2 y4 T& E/ W$ X    else
    7 j3 N1 L6 l4 R, j        w(d1)=l+1;. O2 r7 t" y4 c) j+ ^) ~
            g1=g(s1,:);) V9 J! z, b3 |) q8 m) {/ \
            g1(a1)=2;% j- Z3 q- b* l7 K1 {: X; U6 q
            g(s1,:)=g1;
    1 Y0 U7 a3 a; _0 I        g(:,s1)=g1';* U/ z% g$ q0 ]  R) e  W
            s1=union(s1,d1);
    7 W4 T4 h9 p! o1 i  j' W        l=l+1;
    ' g5 J( B5 D) ^4 x& q  x    end
    & z; A# e# F8 _  t: [end
    * Z8 l% M, k6 K( d9 @. v   
    3 H8 Y# l/ d% n" E4 q- m. p* Bend
    & X* W2 E" E. y/ M' b
    " Z' H8 q! I* o( B. n# {+ e
    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万

    主题

    1312

    听众

    5万

    积分

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

    [LV.Master]伴坛终老

    自我介绍
    数学中国站长

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

    群组数学建模培训课堂1

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

    群组Matlab讨论组

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

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

    ganquanlife 发表于 2013-3-2 17:04
    - R- z( D' M/ \+ V& p帖子里面的笑脸其实是: 和) 不知怎么搞的打出来了就变成了笑脸

    ( M# T1 c% J% a$ r; g8 }0 h这种情况在右侧点击“禁用表情”就没有了。
    数学建模社会化
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-6-6 15:15 , Processed in 0.415414 second(s), 68 queries .

    回顶部