求助：有关递归调用的问题！急用~~

arev326

数据.zip (236.03 KB, 下载次数: 33)

2011-09-21 20:19 上传

点击文件名下载附件

本人想做的是这样一件事：
文件seeds_list.txt为初始节点，文件Adjacency_List.txt为一个以邻接表形式存储的图。根据seeds_list.txt中的初始节点，利用邻接矩阵实现通路查找，先找到与初始节点相同的点，再从此节点开始随机选取与其相连的下一个节点，再将其视为初始节点重复前一步骤，直到无法再找到下一个相连的节点，将此过程中进过的顶点记录下来。最后得到的是N条（N的大小与seeds_list.txt文件中顶点个数相同）通路。

1. use warnings;
   
2. use strict;
   
3. 
   
4. my ($seeds,$j,$k);
   
5. my (@path,@protein);
   
6. 
   
7. open(my $path_list,"Adjacency_List.txt")or die $!;
   
8. while(my $path = <$path_list>){
   
9.         chomp($path);
   
10.         my @temp = split/\t/,$path;
    
11.         my @add = [$temp[0],$temp[1]];
    
12.         push @protein,@add;
    
13. }
    
14. close($path_list);
    
15. 
    
16. open(my $seeds_list,"seeds_list.txt")or die $!;
    
17. open(OUTPUT,">>path_list.txt")or die $!;
    
18. 
    
19. while($seeds = <$seeds_list>){
    
20.         chomp($seeds);
    
21.         $j = 0;
    
22.         &find_next_protein($seeds);
    
23.         # for my $i(0 .. $#path){  #该for循环用于输出数组path
    
24.                 # print OUTPUT "$path[$i]\t";
    
25.         # }
    
26.         print OUTPUT "\n";
    
27.         # @path = ();
    
28. }
    
29. close();
    
30. 
    
31. sub find_next_protein{
    
32.         if($j <= $#protein){
    
33.                 if($protein[$j][0] eq $seeds){
    
34.                         #push @path,$seeds;
    
35.                         print OUTPUT "$seeds\t";
    
36.                         my @next_temp = split/ /,$protein[$j][1];
    
37.                         my $choic_number = int(rand($#next_temp));
    
38.                         $seeds = $next_temp[$choic_number];
    
39.                         if($protein[$j][0] ne $seeds){
    
40.                                 &find_next_protein($seeds);
    
41.                         }
    
42.                 }else{
    
43.                         # if($j <= $#protein){
    
44.                                 $j++;
    
45.                                 &find_next_protein($seeds);
    
46.                         # }
    
47.                 }
    
48.         }
    
49. }

复制代码

捣腾了很久，还是没搞定，急死了~~~
麻烦各位大侠看看，不吝赐教！

Perlvim

一定要用递归调用才行吗？
用其他办法不可以吗？
你说的是不是一个树状结构的遍历循环的问题？
实际应用中节点，通路什么的代表什么数学模型呢？

arev326

回复 2# Perlvim


     用递归应该要容易实现一些。
   我是想在一个无权无向图中找一些通路，起始点就是已经给出的种子节点，然后根据种子节点来随机查找一条通路。
   种子节点在seeds_list.txt中给出，图由邻接表数据文件Adjacency_List.txt给出。数据中的每个节点代表一个蛋白质，边代表蛋白质间的相互作用，而形成的网络就是蛋白质相互作用网。
  各位大仙，帮小弟看看，到底是什么地方错了，或是有什么其他方法实现该操作~~~

zhlong8

这种应用还是换FP来写吧，Perl 格式化数据还可以

arev326

回复 4# zhlong8


    哎~~本人菜鸟一个，只会一点基本的perl~~