改进下程序

zhagnqiang829

我用的是5.14
弱弱问下改成5.14，是不把say改成print就可以。
运行有结果，就输出格式不对。
#!/usr/bin/perl

#use 5.14;
use threads;

open my $d, '1.txt';
open OUT, '>3.txt';
my $s = <$d>;
chomp $s;
my ( $p, $ta, $x, $in, @th ) = ( 20, '', 0 );
my $i = int( length($s) / $p ) + 1;

for ( 1 .. $p ) {
    my $sp = substr( $s, $x, $i );
    $x += $i;
    $in = $ta . $sp;
    push @th, threads->create( \&that );
    $ta = substr( $sp, -7, 7 );
}

my ( @r, %k );
push @r, $_->join for @th;

for (@r) {
    my %h = %$_;
    $k{$_} += $h{$_} for keys %h;
}

#say "$_\t$k{$_}" for sort keys %k;
print OUT "$_\t$k{$_}" for sort keys %k
sub that {
    my %h;
    for ( 0 .. length($in) - 8 ) {
        $h{ substr( $in, $_, 8 ) }++;
    }
    return \%h;
}
回复 10# rubyish


   

rubyish

1. print OUT "$_\t$k{$_}\n" for sort keys %k

复制代码

zhagnqiang829

优化过的代码，运行起来还是
Out of Memory”（内存不足）错误信息。
文件大小1.5G左右，C++能跑起来，就是太慢了。
回复 12# rubyish


   

rubyish

回复 13# zhagnqiang829

1. my ( $p, $ta, $x, $in, @th ) = ( 10, '', 0 );

复制代码

zhagnqiang829

谢谢你，
还是不行。
可能是机子配置问题吧，我自己想办法吧。回复 14# rubyish


   

咏咏672418539

你好，我在算一批序列的六联体个数和频率，但是不太会编程，可以把你的8联体的给我吗？先谢谢啦。在WINDOWS下运行吗？

substr函数

比如下面序列：文件data

1. gtatagaaatcatcgatttttcaggaaaagcaagaaagatcggtttaaaaaattatgaaattttattttaaaaatcgaaaaatcaacaaaaaaaatcgaaattttagaatttttctaaaatcggtagttttcgaaataaaataaaaatgttattttaaaaattattttaaaaaatttcgagaaattcaaaaatcgaacagaaatgttcgattttttgctacagtgcattttttttccacttctacgactttaaaggggggcgcatttatgcggaagggtcttgccgtgcacatagtcatcatttttaacaaagtctgctcaaaattcgcggagatcacgcgtagatcacaattaatttgtgaaaaagcaacagaaatacgaaatatgggataaatattttgtgatcttcatgtgatctacgcgaattttacacagaaactgctaaaaatgatgactaaatgcgcggcaagacccttccgcataaatgctcccccctttaaagtcgtagaagtaggaaaaaatgcactgtatggtgaatttgaaaattttctaaaaatccgaataatttagtatattggggtttcgacgaaaaatgccaatttttcagtaatttctggtcttttctaaaaaaatctacatatttcagatctcttcctccttttccgagttttttcgttattttttctaaactcacagtatattagcgaattttcaacaaattcgacaaaaaattctgactgttttctagaaaaaaaaaaaggcggcatgcgcaagtttcaaaatttgccgcacacctctgggcattaccttgaaagtatcaatggaaaagtgtcgcaagatttaattttctgatattccccccccccctcccgttttacttgtttctaggtttgagtttacttaaattaataacaaaactgggagatatatgctcgaaaattagacagaatctctgcctggacgatttctaggctacgccgtcattcccatagttttcataatctttcggctcagtgatcgaa
   

复制代码

可以用上pypy加速

1. #!/usr/bin/python2
   
2. F = 'data'
   
3. S = open(F).read().rstrip()
   
4. D = {}
   
5. 
   
6. for i in xrange(0, len(S) - 8 + 1):
   
7.     K = S[i:i + 8]
   
8.     D.setdefault(K, 0)
   
9.     D[K] += 1
   
10. 
    
11. for k in sorted(D):
    
12.     print k, D[k]
    

复制代码

最终这个
结果的形式

1. aaaaaaaa 5
   
2. aaaaaaag 1
   
3. aaaaaaat 2
   
4. aaaaaagg 1
   
5. aaaaaatc 2
   
6. aaaaaatg 1
   
7. aaaaaatt 3
   
8. aaaaagca 1
   
9. aaaaaggc 1
   
10. aaaaatca 1
    
11. aaaaatcc 1
    
12. aaaaatcg 3
    
13. ............

复制代码