QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3618|回复: 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))是同一个意思?( L! K: S  ^' f1 ]1 b6 o: ^; R( n
    我个人理解,g==1除非g是空的,否则g==1都会返回一个逻辑矩阵才对啊,不会变成空的。& D# d% X  _) A
    function [w] = bfsf1(g)
    % i1 n9 h9 h4 v) ?5 o' J%广度优先搜索9 r9 Z* x' s9 \4 K% f: A! j* g* `% B
    %g:邻接矩阵。w顶点标号
    * ]( p6 V8 U' @& r) Vn=size(g,1);- Z, S5 b+ k. ~- ]* K  b
    w=zeros(1,n);/ c: D4 v' `  n- j* @6 W
    l=0;
    3 p( n3 l- x$ X6 e5 W; O7 f/ l$ y" wv=1;
    $ v7 g5 o1 h% R5 Z: H$ x& Ka1=find(g(v,:)==1);
    ; l1 D# B+ l5 Gg(v,a1)=2;
    ) d0 }1 b7 _* Z5 F7 X+ `g(a1,v)=2;4 {4 y% k* X  ?3 A2 P9 H3 e7 J; p
    w(a1)=l+1;
    3 d/ I! A4 |/ G1 ~s1=union(a1,v);
    # ~! W5 D. A$ R# ?' U5 W% bl=l+1;8 f  u* E. s% p9 Z- Q5 I
    while~isempty(g==1)2 q: ?. Y) |: h! p7 W
       a1=find(g(s1,:)==1);) t% \. x# t4 `$ Y6 {" m' M( e
        t=length(s1);
    ! l5 v- w- x' A- b: R" h8 d5 ~/ J    d=[];
    " u5 k; b& Z& g% l+ q, a# l2 _5 F7 h    for i=1:length(a1)
    ; v& d% v/ z, {1 X7 K        if a1(i)/t>floor(a1(i)/t)
    8 t1 X+ K/ a3 t, [            t2=floor(a1(i)/t)+1;
    , b+ w2 F) l! f3 Z6 [5 c1 C        else6 S( a: D. h, X7 u+ G- [5 b( z
                t2=floor(a1(i)/t);
    / }# S2 Y, T: S* X  ]$ H1 ?        end4 }. a  B& J: i: f- M$ R' p7 @& s
            if isempty(intersect(d,t2)): h$ W: ^6 c% q$ I7 I. O$ Q
                d=union(d,t2);% S) n  Z" {1 R' r7 H1 |7 p+ g
            end& X/ t% P9 S- z7 x& |7 ^
        end
    . l* o: ^. P) s" O* Z, _    d1=setdiff(d,s1);( e( D5 A$ n( m
        if isempty(d1)3 l+ {- f5 Z; v: X
            break;- a" V! o0 o3 e0 x6 U) u# w/ t; ?
        else9 b* Z& S2 [. {6 e$ Q4 j
            w(d1)=l+1;& u3 ^6 b2 k, Y9 Q6 b8 O
            g1=g(s1,:);5 L/ P; u* p( a/ _+ U; N: d; r- w
            g1(a1)=2;$ q/ P3 j/ W! `% r! g
            g(s1,:)=g1;+ H7 j: p( P  @) G
            g(:,s1)=g1';
    - p6 l+ E$ a7 W% d        s1=union(s1,d1);
    1 ]5 N  K) O* `7 J. s        l=l+1;/ _0 ?+ p9 `& G! f/ s' Z* |9 E0 T
        end
    7 P1 g# |7 }& \+ W/ ~* A* v- rend
    $ Q0 p! S9 D4 b9 ~, I  j   
    ! E* e" |- x' K  V' V4 Lend) G8 q6 B) {8 I: M
      G9 i' O1 O3 n7 o! r9 _
    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 ; t% a! N/ F& t# B, L/ T7 I
    帖子里面的笑脸其实是: 和) 不知怎么搞的打出来了就变成了笑脸
    " H/ Q1 S; K* ^+ Z- S
    这种情况在右侧点击“禁用表情”就没有了。
    数学建模社会化
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-6-7 23:00 , Processed in 0.536292 second(s), 68 queries .

    回顶部