求大神帮写一个去除交集的perl脚本！！

reyleon

回复 10# pitonas


    别开玩笑了,楼主不是把数据贴出来了嘛,你去整个代码吧,我是没看懂楼主的意思,SOGA,整不出来呀

pitonas

大牛你整个把图上数据扣出来代码嘛{:2_172:}

reyleon 发表于 2014-01-20 12:15
回复 10# pitonas
SOGA

jzp520520

回复 8# xiumu2280
大神！昨天那代码那么强，这还能难倒你？等你的精彩代码(一定要用贴图的方式哦)！

   

xiumu2280

1. open (IN,"C:\\Users\\Desktop\\Test\\2014_1_14\\Noname5.txt");
   
2. open (IN1,"C:\\Users\\Desktop\\Test\\2014_1_14\\Noname3.txt");
   
3. 
   
4. my %hash_goal;
   
5. map {chomp;push @{$hash_goal{(split)[0]}},((split)[1]..(split)[2])}<IN>;
   
6. ##################
   
7. map {
   
8.         chomp;
   
9.         my $line = $_;
   
10.         my ($zone,$start,$end) = split /\s+/,$line;
    
11.         my %hash_NB;
    
12.         @hash_NB{($start..$end)}=1;
    
13.         !(%hash_NB ~~ @{$hash_goal{$zone}}) && print "$line\n";
    
14. }<IN1>;

复制代码

超级耗内存版本

pitonas

看得我都激动了...好厉害的小伙伴们, 这个必须赞一个
{:2_172:}

回复 14# xiumu2280


   

jzp520520

回复 14# xiumu2280
大神，这一句好像有些问题：my ($zone,$start,$end) = split,$line;本人测试提示$line是处于空上下文的，没有作用。
若my ($zone,$start,$end) = （split,$line）;这样写，是没有语法错误，不过在本例中$line也是没有赋值给变量，也无任何作用。

   

xiumu2280

的确，那个split的写法是错误的···
感谢，大神~{:3_187:} 回复 16# jzp520520


   

felix0608

顶，求大神帮忙！！

pitonas

{:2_177:}头好晕 @@
{:2_172:}小伙伴们, 测试测试?
            # if [18656   18756] [18075   18656] is overlay
            # USE:
            # if ( $r2[$y][0] > $x->[1] ) {
            # elsif ( $r2[$y][1] < $x->[0] ) { next }
            
            # if [18656   18756] [18075   18656] is not overlay
            # USE:
            # if ( $r2[$y][0] >= $x->[1] ) {
            # elsif ( $r2[$y][1] <= $x->[0] ) { next }

1. #!/usr/bin/perl
   
2. open my $one, 'doc1';
   
3. open my $two, 'doc2';
   
4. my %doc1; map { @_ = split; push @{ $doc1{ shift @_ } }, [@_] } <$one>;
   
5. my %doc2; map { @_ = split; push @{ $doc2{ shift @_ } }, [@_] } <$two>;
   
6. 
   
7. while ( my ( $s, $r ) = each %doc1 ) {
   
8.     next unless $doc2{$s};
   
9.     my @r1 = @$r;
   
10.     my @r2 = @{ $doc2{$s} };
    
11.     my $i  = 0;
    
12.     for my $x (@r1) {
    
13.         my $print = 1;
    
14.         for my $y ( $i .. $#r2 ) {
    
15.             if ( $r2[$y][0] >= $x->[1] ) {
    
16.                 $i = $y;
    
17.                 last;
    
18.             }
    
19.             elsif ( $r2[$y][1] <= $x->[0] ) { next }
    
20.             else                           { $print = 0 }
    
21.         }
    
22.         print join( "\t", $s, @$x ), "\n" if $print;
    
23.     }
    
24. }
    
25. 
    
26. __DATA__
    
27. 
    

复制代码

pitonas

jason680
大神 ! 帮我看看代码， 它是否有问题 ?