QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3842|回复: 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))是同一个意思?
    4 E/ F, i. t; S' U; k+ x, j7 ~' x我个人理解,g==1除非g是空的,否则g==1都会返回一个逻辑矩阵才对啊,不会变成空的。
    ) b: a$ J( F2 F8 Ifunction [w] = bfsf1(g)& G8 j9 n: ~" M* k
    %广度优先搜索4 g8 H" Z# q( p
    %g:邻接矩阵。w顶点标号
    * L( e8 ]; a7 X8 C! Cn=size(g,1);/ `/ A* {9 g4 z) a( u# U6 C; t
    w=zeros(1,n);
    ( W1 I# C& o& F% u( z* y9 {l=0;8 n; K/ S* a3 y" ~& D* f
    v=1;2 [$ s. {& t$ y" V0 F; H
    a1=find(g(v,:)==1);5 N4 P2 i- [1 q" G0 D; }2 f+ y
    g(v,a1)=2;
    ; D# D$ x! m7 O% Eg(a1,v)=2;
    " {  a! V, W" o% V* xw(a1)=l+1;, `! r% T3 s5 }
    s1=union(a1,v);
    6 _1 j5 u& o, _+ H1 El=l+1;( t9 K+ s( Q, f& M6 z( p
    while~isempty(g==1)+ g' s  L( ~+ g0 c- E$ O
       a1=find(g(s1,:)==1);
    ! U$ A. [9 T4 Q    t=length(s1);
    / X6 v4 t; A, B    d=[];6 F+ U) X3 L! ^( U5 o2 i7 U
        for i=1:length(a1)
    , M+ W, }: m7 r7 j, u6 H% N        if a1(i)/t>floor(a1(i)/t)
    * g. `+ ~& {9 a9 _: L            t2=floor(a1(i)/t)+1;
    . z& t6 \0 n% ]2 N3 \0 W' `        else6 Y7 j9 o! x9 P0 O2 z/ f8 f
                t2=floor(a1(i)/t);
    % s8 D0 e& L1 K7 Z        end
    : l- Z, U& Y: d! }1 J4 u& ^( G        if isempty(intersect(d,t2))& a; m3 y; x# n
                d=union(d,t2);
    0 P0 y! t% f8 Y, A/ v        end  D( o' q( M! n8 s, U/ {1 l% V
        end
    9 H6 o' I* o1 G' y6 x    d1=setdiff(d,s1);
    9 T6 M, k8 H0 Y- Q    if isempty(d1)
    # u1 r6 y5 o" V# g# ]9 _        break;
      }& [9 K/ u; u0 k6 X( m; {    else- u7 u7 f8 P8 M. m
            w(d1)=l+1;
    2 w4 F5 n6 a+ |5 `1 g! m4 y        g1=g(s1,:);7 B* W/ C# j  W  ~0 A! T* K
            g1(a1)=2;
    9 o" R$ J4 ]  q. g* ~1 }        g(s1,:)=g1;
    ' R* L/ R$ B+ t3 ^$ m        g(:,s1)=g1';
    9 i( J3 k( Y* d4 A- S        s1=union(s1,d1);
    ; b. ?$ C* ?% `! j! e- g3 m        l=l+1;
    ! b. }- W4 C+ o" v/ F* @7 q    end
    # P3 S- ]1 T, C/ R) u! s# Aend
    " h( P" x( V- S( i9 Q7 h( X   
    $ {3 d6 [1 h% p+ z) U: Gend
    % @1 o2 W4 t( M9 z$ r% A
    # Y7 F+ L1 N8 |$ H5 A% r2 A' |
    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 & _, ?2 C) |+ J) z6 m- b6 E5 j$ X
    帖子里面的笑脸其实是: 和) 不知怎么搞的打出来了就变成了笑脸

    * j8 P* F9 j: b/ I2 g) D) O这种情况在右侧点击“禁用表情”就没有了。
    数学建模社会化
    回复

    使用道具 举报

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

    fastpost qq
    收缩
    • 电话咨询

    • 04714969085

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-10-31 19:09 , Processed in 0.284524 second(s), 65 queries .

    回顶部