查看: 5964|回复: 3

[代码资源] k-shell代码

字体大小: 正常放大

2983 主题	142 听众	9762 积分

升级 95.24%

TA的每日心情

	开心 2017-1-9 14:34

签到天数: 272 天

[LV.8]以坛为家I

自我介绍: 吃吃吃

群组: 乐考无忧

群组: 2014国赛优秀论文解析

群组: 2016美赛冲刺培训

群组: 2016国赛优秀论文解析

群组: 2016国赛备战群组

电梯直达

1^#

发表于 2015-12-26 10:32 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

本帖最后由吃苹果的梨于 2015-12-26 10:37 编辑

%function [ret_m]=k_shell(mixedsig)
%邻接矩阵
%  A=load('cdbBA_4000_5_0_.txt'); %里边的数据结构好像是每2个一组，例如 3 4 表示一组，表示点3和节点4有连接。
%
%
%  TT=A(:, 1:2)+1;
%  maxtt=max(max(TT));
%  mixedsig=zeros(max(max(TT)));
%  len=length(TT);
%  for i=1:len
%    mixedsig(TT(i,1),TT(i,2))=1;
%    mixedsig(TT(i,2),TT(i,1))=1;
%  end
%  disp('已转化为邻接矩阵！')

[XY]=distri(mixedsig);
mixedsig=[0,1,1,1;1,0,0,1;1,0,0,1;1,1,1,0]
function [ret_m]=kshell(mixedsig)
mixedsig=[0,1,0,1;1,0,0,1;0,0,0,1;1,1,1,0];
mixedsig=[0,1,0,1,0;1,0,0,1,0;0,0,0,1,1;1,1,1,0,0;0,0,1,0,0];

zzz=0;
j=1; %第j层
t=1;  %
tad=mixedsig;  %邻接矩阵
k=0;
ii=1;
ret_m=[];  %返回值；  i，j  ：第i层：所有节点
len=length(mixedsig);  %矩阵长度
tt=1;
n=1;
nn=1;
m_t=[];
m_tt=[];
m_ttt=[];
p=1;
while (tt==1) %控制最高层结束
sss=sum(sum(tad))/(maxtt*maxtt)
if sum(sum(tad))==0  %所有元素为0，则退出
      break;
end
t=1;  %控制第j层  计算
while(t==1)  %  每循环一次，tad改变，去掉度数小于j的节点；直到没有度数小于j的节点
      t=0;
      ii=1; %第j层第ii个节点
      for i=1:len  %从矩阵1至len行,  去掉小于j的节点
         lll=i/len
         k=sum(tad(i, : )); %计算i行度数
         if k==0  %度数为0，下一个i值
            % t=1 ;  %

            continue;
         elseif k<=j  %度数小于j层
            t=1; %控制下一次还要循环
            tad(i, : )=0;  %i节点加到j层，将度数至为0，所有i行值至为0
            tad(:,i)=0; %相应i列至为0

            ret_m(j,ii)=i; %将i节点加到j层

            m_ttt=union(ret_m(j,ii),m_ttt);
            for n=1:len %判断其他行有无因为将tad(:,i)元素设为0 而所有行变0
                  if sum(tad(n, : ))==0

                     ii=ii+1;
                     ret_m(j,ii)=n;
                     m_t=n;
                     m_tt = intersect(m_ttt,m_t);
                     if  length(m_tt)==0  %length(m_t)~=0 &&

                        m_ttt=union(ret_m(j,ii),m_ttt);
                        m_t=[];
                     end
                  end

            end %for n=1:len
         end % if k==0
         ii=ii+1;
      end %end of  i=1:len

end  % end of while(t==1)
j=j+1;
ii=1;
end % end of while(tt==1)
disp(ret_m);