{重新编辑，抱歉}求高手指点（比较有挑战的题目），使用perl提取两个文件中不同部分

kk861123

回复 10# yanglc2013

try:

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. use Data::Dumper;
   
5. use utf8;
   
6. 
   
7. my $fileA = "D:\\1\\图A.fasta";
   
8. my $fileB = "D:\\1\\图B";
   
9. open A, $fileA or die $!;
   
10. open B, $fileB or die $!;
    
11. 
    
12. my %hBMap;
    
13. while (<B>) {
    
14.     next if /^$/;
    
15.     chomp;
    
16.     $hBMap{$_}++;
    
17. }
    
18. #print Dumper \%hBMap;
    
19. 
    
20. {
    
21.     local $/ = '>';
    
22.     my $line = <A>;
    
23.     my ($title, $data);
    
24.     while ( defined($line = <A>) ) {
    
25.         chomp($line);
    
26.         $title = '>' . (split /\n/,$line,2)[0];
    
27.         print ">$line" unless exists $hBMap{$title};
    
28.     }
    
29.    
    
30. }
    
31.   
    
32. exit;

复制代码

yanglc2013

回复 11# kk861123

牛大哥！真谢谢您。您帮了我两次了。


   

kk861123

回复 12# yanglc2013


    不客气，暂时有时间，互相学习

yanglc2013

回复 13# kk861123
我已经学了7天的perl了,能加你QQ（515465501）吗，交个朋友。


   

笨鱼月色

完了，我看不到楼主说的神马意思，怎么破...

heritrix

恩恩，学习了

heritrix

帮顶

dqw233

1. use v5.14;
   
2. use warnings;
   
3. use autodie;
   
4. 
   
5. open(my $a_fd, '<', 'a.txt');
   
6. open(my $b_fd, '<', 'b.txt');
   
7. 
   
8. my %id_in_b;
   
9. 
   
10. while (<$b_fd>) {
    
11.     chomp;
    
12.     $id_in_b{$_} = 1;
    
13. }
    
14. 
    
15. my $need_print = 0;
    
16. while (<$a_fd>) {
    
17.     chomp;
    
18.     if (m/^>/xms) {
    
19.         $need_print = 0;
    
20.         if (not exists($id_in_b{$_})) {
    
21.             $need_print = 1;
    
22.         }
    
23.     }
    
24.     say if $need_print;
    
25. }

复制代码