perl

MMMIX

@b114213903

MMMIX

为啥不自己先写一个呢? 新手就应该多练习嘛

MMMIX

回复 8# super_two


    不能运行也罢, 贴上来看看.

MMMIX

回复 11# super_two


    其实, 你只需要逐个读取序列, 对于每个 id, 只保留最长的那个序列即可.

MMMIX

回复 18# super_two


    真心建议你找些 Perl 编程的书好好看看, 不过估计你是听不进去的.

无论如何, 下述脚本演示了解决你的问题的一种方法. 另外,如果在你的输入文件中,id相同的序列顺序出现,那么也可以不用hash.

1. 
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. use v5.14;
   
6. 
   
7. my %record;
   
8. while (<DATA>) {
   
9.     my @f = split;
   
10.     $f[0] =~ s/(?:\.\d*)?$//;
    
11.     if (! $record{$f[0]} || $record{$f[0]}->[0] < length($f[1])) {
    
12.         $record{$f[0]} = [length $f[1], $f[1]];
    
13.     }
    
14. }
    
15. 
    
16. say "output version 1:";
    
17. while (my ($id, $seq) = each %record) {
    
18.     say "$id\t", $seq->[1];
    
19. }
    
20. 
    
21. say "output version 2:";
    
22. for my $id (sort keys %record) {
    
23.     printf "%-12s %s\n", $id, $record{$id}->[1];
    
24. }
    
25. 
    
26. __DATA__
    
27. qwer.1         AFRTYUIFGHJKLVBNM
    
28. qwer.2         BVXNVFGSFYEBCSHB
    
29. qwer.3         HDFKSHFGSERYFIEURHFSUFDSHVBSJEUABFUHABFCAHFBC
    
30. rtyuip00.1    AFHBVSFHUACFKUSHDBAKFHAKUFHSADKFUA
    
31. hhjkl.1         JDNVKHFBAKHFAKFAFJNSADFJAS
    
32. hhjk.2         HFSDHNFKANFAKFIJI
    
33. hhjk.3         JNFAJNFALSDFLMAD
    
34. hhjk.4         KJGSEGJOAKFFDSMFAPOKEF
    

复制代码

MMMIX

回复 20# 聆雨淋夜


    对于我不熟悉的语言来说, 我就没有发现可读性好的.

MMMIX

回复 24# super_two


    你是修改过人家的脚本吧?

MMMIX

回复 27# super_two


    你读懂那个脚本了么? 懂了的话就不会有疑问了.

MMMIX

回复 29# b114213903


    我觉得他是光看了脚本的输出, 但并没有检查脚本写入的文件.