tr///

papagogogo

$dna1="2220122201022001110111122";


foreach (0..0){
  $D[$_]= substr($dna1,30*$_,30);
  $E[$_]=($D[$_]=~ tr/1/1/);
  $C[$_]=($D[$_]=~ tr/01/01/s);
  $B[$_]=($D[$_]=~ tr/1/1/);
  print $E[$_];
  print "\n";
  print $C[$_];
  print "\n";
  print $B[$_];
  
}
输出结果是
9
15
4


我希望计算01出现的个数，连续的01。而非0或1。可是如果我写成这样，最后print $B是我想要的结果，可是我明明写的是B是tr/1/1/啊。C我写成 tr/01/01/s反而得不到我要的结果反而输出的是0或者1的个数。。。希望大家给我解释一下我如何在这段字符串中匹配计算出01的个数也就是4

papagogogo

$dna1="100010000100000111111111111111111110111100";
$W=($dna1=~ tr/01/01/s);
print $dna1;
print "\n";


结果是1010101010

他是将原字符串100010000100000111111111111111111110111100中连续的0变成单个0，连续的1变成单个1，为何会有这种效果呢，如何解读，谢谢

ly5066113

回复 1# papagogogo


try:

1. perl -e '$_="2220122201022001110111122";print s/01/01/g;'

复制代码

chenhao392

回复 3# ly5066113


    似乎没有那么简单...
    Bioinfo 还有效率问题，tr/// 返回计数该是比s///快很多的...
    看楼主数据量有多大，能不能忍吧...

    我试想可以把多行merge到一行，然后 index 和reindex..
    不过没有benchmark，不知效率...

jason680

回复 1# papagogogo

$ perl -e '$_="2220122201022001110111122";$c++ while(m/01/g);print "$c\n"'
4

$ perl -e '$_="2220122201022001110111122";$c++ while(m/0/g);print "$c\n"'
6

$ perl -e '$_="2220122201022001110111122";$c++ while(m/1/g);print "$c\n"'
9

   

chenhao392

回复 5# jason680


    赞！学习了！

cat tr_test.pl

1. 
   
2. #!/usr/bin/perl
   
3. 
   
4. use strict;
   
5. use warnings;
   
6. 
   
7. my $data="01" x 10000000;
   
8. my $one=$data=~tr/01/01/;
   
9. print "$one\n";
   

复制代码

cat s_test.pl

1. 
   
2. #!/usr/bin/perl
   
3. 
   
4. use strict;
   
5. use warnings;
   
6. 
   
7. my $data="01" x 10000000;
   
8. my $one=$data=~s/(01)/(01)/g;
   
9. print "$one\n";
   

复制代码

cat m_test.pl

1. 
   
2. #!/usr/bin/perl
   
3. 
   
4. use strict;
   
5. use warnings;
   
6. 
   
7. my $data="01" x 10000000;
   
8. my $one=0;
   
9. while($data=~m/01/g){
   
10.         $one++;
    
11. }
    
12. 
    
13. print "$one\n";

复制代码

结果:

1. time perl tr_test.pl
   
2. 20000000
   
3. 
   
4. real        0m0.061s
   
5. user        0m0.032s
   
6. sys        0m0.026s
   
7. 
   
8. 
   
9. 
   
10. time perl s_test.pl
    
11. 10000000
    
12. 
    
13. real        0m3.279s
    
14. user        0m3.192s
    
15. sys        0m0.072s
    
16. 
    
17. 
    
18. 
    
19. time perl m_test.pl
    
20. 10000000
    
21. 
    
22. real        0m1.233s
    
23. user        0m1.189s
    
24. sys        0m0.027s
    

复制代码

substr函数

tr/012/012/
tr/0/0/ + tr/1/1/ + tr/2/2/