[问题求助] 求助关于求图割点的代码，哪里错了？怎么改正啊？

[复制链接]

字体大小: 正常放大

ganquanlife

9 主题	7 听众	188 积分

升级 44%

TA的每日心情

	慵懒 2013-8-23 15:15

签到天数: 84 天

[LV.6]常住居民II

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

群组: 西安交大数学建模

群组: 数学建摸协会

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

电梯直达

发表于 2013-3-3 14:54 |只看该作者 |正序浏览

|招呼Ta 关注Ta

2体力

function [nc] = ncutf(g)
%求割点的算法 g为邻接矩阵 nc为割点的集合
n=size(g,1);
if n>=3
a=sum(g);
b=sum(a==2);
if b==n
      fprintf('本图为圈，无割点。n')
      nc=0;
end
else
[w,k]=dfs3(g);
%nc=[];
nc=isncf(w,k);
n=size(g,1);
for i=1:n
      for j=1:n
         if w(i,j)>1
            if k(i)>k(j)
                  g(i,j)=2;
            else
                  g(i,j)=3;
            end
         end
      end
end

for i=1:n
      f1=find(g(i,

==2);
f2=find(g(i,

==3);
      f=union(f1,f2);
      l(i)=min([k(f) k(i)]);
end

for i=1:n
      for j=1:n
         if g(i,j)==3 & k(i)>1&l(j)>=k(i)
            nc=union(i,nc);
         end
      end
end
end
end

function nc=isncf(w,k)
nc=[];
t=zeros(size(w));
n=size(w,1);
a=find(w~=0);
for i=1:length(a)
      d(i)=w(a(i));
      if a(i)/n>floor(a(i)/n)
      t(i)=floor(a(i)/n)+1;
      else
         t(i)=floor(a(i)/n);
      end
      t1(i)=mod(a(i),n);
      if t1(i)==0
         t1(i)=n;
      end
end
[b,c]=sort(d);
p=[1];pc=0;
for i=1:length(a)
      if k(t1(c(i)))<k(t(c(i)))
         p=union(p,t(c(i)));
         t(t1(c(i)),t(c(i)))=3;
      end
      if pc==0
         tc=isempty(setdiff([1:n],p));
         if tc
            t0=sum(t(1,