BA.rar_BA模型_C++中BA模型_ba_复杂网络_复杂网络 程序

#include <iostream> using namespace std; void main() { int m=5; /*初始态为全耦合网络，m=5,初始时有m个点，每次连接m条边*/ int N=20; /* N——网络的大小*/ int Nodes[20][20]={0}; /*Nodes——N*N矩阵，网络中各点的连接情况*/ int i ,j,v,w ; int t ; int n,p,q,b; int sum=20; int list[300]={1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5}; /*list——另外生成的一个向量， 表示了各点连接的情况，随机从该向量中取点，相当于按照度的概率取点*/ for(i=0;i<5;i++) /*生成的初始网络*/ { Nodes[i][i]=1-m; for(j=0;j<m;j++) { if(j!=i) {Nodes[i][j]=1; Nodes[j][i]=1; } } } cout<<"初始矩阵："<<endl; for(q=0;q<5;q++) for(int b=0;b<5;b++) { cout<<Nodes[q][b]<<" "; if(b==4) cout<<endl; } for(n=5;n<N;n++) { for(i=0;i<5;i++) { list[sum]=n+1; sum++; } for(i=0;i<5;i++) { w=1; v=sum-5; while(w<2) { p=rand(); if (p>=0&&p<v) { if(p==0) { list[sum++]=1; if(Nodes[n][0]==0) { Nodes[0][n]=1; Nodes[n][0]=1; Nodes[n][n]=Nodes[n][n]-1; Nodes[0][0]=Nodes[0][0]-1; } } else { if(Nodes[n][list[p]-1]==0) { if(list[p]-1!=n) { list[sum++]=list[p]; Nodes[n][list[p]-1]=1; Nodes[list[p]-1][n]=1; Nodes[list[p]-1][list[p]-1]=Nodes[list[p]-1][list[p]-1]-1; Nodes[n][n]=Nodes[n][n]-1; w=3; } } else w=1; } } else w=1; } } } cout<<"最后矩阵:"<<endl; for(q=0;q<N;q++) { for(int b=0;b<N;b++) { cout<<Nodes[q][b]<<" " ; if(b==(N-1)) cout<<endl; } } cout<<endl<<"LIST:"<<endl; for(i=0;i<200;i++) cout<<list[i]<<" "; cout<<"结束"; }