perl区间交集

yinyuemi

回复 1# 清泉一边1

1. #/usr/bin/env perl
   
2. 
   
3. use strict;
   
4. use Data::Dumper;
   
5. 
   
6. my %h;
   
7. 
   
8. while(<DATA>){
   
9.     chomp;
   
10.     my @a = split;
    
11.     my $k = substr($a[0],0,1)."\t".$a[3];
    
12.     $h{$k}{min} = $a[4] if $h{$k}{min} > $a[4] or ! defined $h{$k}{min};
    
13.     $h{$k}{$a[4]} = $a[5];
    
14. }
    
15. my %max;
    
16. my %min;
    
17. my $start;
    
18. my $end;
    
19. foreach my $k(keys %h){
    
20.     foreach my $t (sort{$a<=>$b}grep{/\d/}keys %{$h{$k}}){
    
21.         if($t == $h{$k}{min} || $t > $max{$k}{$start}->[-1]){
    
22.             $max{$k}{$t} = [$t,$h{$k}{$t}];
    
23.             $min{$k}{$t} = [$h{$k}{$t}-$t,$t."-".$h{$k}{$t}];
    
24.             $start = $t;
    
25.         }elsif($t <= $max{$k}{$start}->[-1] ){  
    
26.             $min{$k}{$start} = [$end-$t,$t."-".$end] if $min{$k}{$start}->[0] > ($end-$t);
    
27.                         $max{$k}{$start} = [$start,$h{$k}{$t}] if $h{$k}{$t} > $max{$k}{$start}->[-1];
    
28.         }
    
29.                 $end = $h{$k}{$t};
    
30.     }
    
31. }
    
32. 
    
33. print "Max Gap\n";
    
34. foreach my $k (keys %max){
    
35.         print $k,"\t";
    
36.         print join "_",map{join "-",@{$max{$k}{$_}}}(sort{$a<=>$b}keys %{$max{$k}});
    
37.         print "\n";
    
38. }
    
39. print "Min Gap\n";
    
40. foreach my $k (keys %min){
    
41.         print $k,"\t";
    
42.         print join "_",map{@{$min{$k}{$_}}[1]}(sort{$a<=>$b}keys %{$min{$k}});
    
43.         print "\n";
    
44. }
    
45. 
    
46. 
    
47. __DATA__
    
48. a.1        1        3        1_b        2        8
    
49. a.2        2        5        1_b        5        9
    
50. a.3        3        6        1_b        4        7
    
51. a.4        2        4        1_b        12        15
    
52. c.1        3        8        1_d        1        7
    
53. c.2        4        8        1_d        6        8
    
54. 
    
55. -----output-----
    
56. Max Gap
    
57. a       1_b     2-9_12-15
    
58. c       1_d     1-8
    
59. Min Gap
    
60. a       1_b     5-7_12-15
    
61. c       1_d     6-7
    
62. 
    

复制代码