chemperl无向图中找简单闭合回路perl的实现(DOC).docVIP

chemperl无向图中找简单闭合回路perl的实现(DOC)

这里我们使用的是广度优先算法，所以我们得到的是图（b）.广度优先生成树BFS

通过BFS算法，把边的权重改为-1；

------------------------------------------------------------------------------------------------------------------------------------------------------------------------

结果当然要放到正文里了。

------------------------------------------------

下面是连通信息

12

15

18

23

213

36

314

46

47

415

57

511

79

816

910

918

1012

1017

1112

1121

1222

1316

1323

1324

1625

1626

1720

1727

1819

1828

1920

1929

2030

------------------------------------------------------

这里是连通信息表达的分子式

-----------------------------------------------------------------------------

这里是我用perl写的脚本寻找环的编号，

#!/usr/bin/perl-w

usestrict;

#N=30E=3330个点，33条边

#构建邻接矩阵

my$x;

my$y;

my@matrix;

for($x=1;$x=30;$x++)

{

for($y=1;$y=30;$y++)

{

$matrix[$x][$y]=0;

}

}

openFH,F:\\C\\bond;

while(FH)

{

/(\d+)\s+(\d+)/;

$matrix[$1][$2]=1;

$matrix[$2][$1]=1;

}

closeFH;

#for($x=1;$x=30;$x++)

#{

#for($y=1;$y=30;$y++)

#{

#print$matrix[$x][$y];

#}

#print\n;

#}

#打印邻接矩阵

my@color;#对每个顶点u∈V，其色彩存储于变量color中.

my@parent;#结点u的父母存于变量π[u]=parent中

#如果u没有父母(例如u=s或u还没有被检索到)，则π[u]=NIL

##初始化

#把所有的点初始化为白色

for(1..30)

{

$color[$_]=w;

}

my@box;#用来存放灰色的顶点

$parent[1]=0;#表示1号原子没有父母，也就是说任命1号原子为root，祖宗

$color[1]=g;#把1号原子的颜色变为灰色ggrey，

#黑色bBlack;灰色greyg

push@box,1;#把灰色的1号原子放到盒子中

#写一个子例程，用来返回一个原子的相连原子

#printjoin\n,adj(2);

subadj

{

my@adjatoms=();

my$adj_atom=$_[0];

openFH,F:\\C\\bond;

while(FH)

{

if($_=~/\b$adj_atom\b\s+(\d+)/||$_=~/(\d+)\s+\b$adj_atom\b/)

{

push@adjatoms,$1;

}

}

return@adjatoms;

}

my@adj;#Adj[u]表示图中和u相邻的节点集合

while(@box!=0)#当盒子中没有灰色的原子时，停止运行

{

my$greyatom=$box[0];

my@atoms=adj($greyatom);

#print\@a