perl區(qū)間交集問題（急�。�

清泉一邊1

第2列三列沒用,那么用每一行的第5列和第6列與下面各列的進(jìn)行比較，5列為起始坐標(biāo)，6列為終止坐標(biāo)。輸出他們的最大交集（即如果坐標(biāo)區(qū)間有重疊,則合并之并取其最小的起始坐標(biāo)和最大的終止坐標(biāo)）
a.1        1        3        b        2        8
a.2        2        5        b        5        9
a.3        3        6        b        4        7
a.4        2        4        b        12        15
c.1        3        8        d        1        7
c.2        4        8        d        6        8
理想輸出
a        b=1        2        9
a        b=2        12        15
c        d        1        8
"="代表b出現(xiàn)的次數(shù)
另外還需要輸出另一個文件，即最小交集（如果坐標(biāo)區(qū)間有重疊,則合并之并取其最大的起始坐標(biāo)和最小的終止坐標(biāo)）
理想結(jié)果是
a        b=1        5        7
a        b=2        12        15
c        d        6        7
求大神幫忙�。�！

jjboy110

且不另你年少無知，當(dāng)你得不到答案的時候多了，你才能知道其實提問是有提問的智慧的，問題最終都是自己解決，如果要別人幫你解決，請花錢。

清泉一邊1

回復(fù) 2# jjboy110
一個腳本哦。。。。

   

xiumu2280

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my %hash_data;
   
6. map {chomp;$_=~/([a-z]+)./;push @{$hash_data{$1}},[split,1,0,0]} <DATA>;
   
7. 
   
8. foreach my $k (sort {$a cmp $b}keys %hash_data) {       
   
9.         my @ori = shift@{$hash_data{$k}};
   
10.         $ori[0][7]=$ori[0][4];
    
11.         $ori[0][8]=$ori[0][5];
    
12.         foreach my $line (@{$hash_data{$k}}) {
    
13.                 $line->[7]=$line->[4];
    
14.                 $line->[8]=$line->[5];
    
15.                 if ($line->[4] <= $ori[-1][5]) {
    
16.                          if ($ori[-1][5] <= $line->[5]) {
    
17.                                 $ori[-1][5] = $line->[5];
    
18.                                 $ori[-1][6]++;
    
19.                                 if ($ori[-1][7] <= $line->[7]) {
    
20.                                         $ori[-1][7] = $line->[7];
    
21.                                 }
    
22.                                 if ($ori[-1][8] >= $line->[8]) {
    
23.                                         $ori[-1][8] = $line->[8];
    
24.                                 }
    
25.                         }else{
    
26.                                 if ($ori[-1][7] <= $line->[7]) {
    
27.                                         $ori[-1][7] = $line->[7];
    
28.                                 }
    
29.                                 if ($ori[-1][8] >= $line->[8]) {
    
30.                                         $ori[-1][8] = $line->[8];
    
31.                                 }
    
32.                                 $ori[-1][6]++;
    
33.                         }
    
34.                 }else{
    
35.                         push @ori,$line;
    
36.                 }
    
37.         }
    
38.         foreach my $li (@ori) {
    
39.                 print "$k\t$li->[3]=$li->[6]\t$li->[4]\t$li->[5]\n";
    
40.         }
    
41.         foreach my $li (@ori) {
    
42.                 print "$k\t$li->[3]=$li->[6]\t$li->[7]\t$li->[8]\n";
    
43.         }
    
44. }
    
45. 
    
46. __DATA__
    
47. a.1        1        3        b        2        8
    
48. a.2        2        5        b        5        9
    
49. a.3        3        6        b        4        7
    
50. a.4        2        4        b        12        15
    
51. c.1        3        8        d        1        7
    
52. c.2        4        8        d        6        8

復(fù)制代碼

1. a       b=3     2       9
   
2. a       b=1     12      15
   
3. a       b=3     5       7
   
4. a       b=1     12      15
   
5. c       d=2     1       8
   
6. c       d=2     6       7

復(fù)制代碼

結(jié)果和你的有些出入，不清楚你說的b出現(xiàn)的次數(shù)具體指什么，而且格式你需要自己調(diào)整下。

清泉一邊1

太感謝了�。�！回復(fù) 4# xiumu2280


   

yestreenstars

稍微長了點~

1. #!/usr/bin/perl
   
2. use 5.010;
   
3. my ( %hash, %max, %min );
   
4. while (<DATA>) {
   
5.         @_ = split;
   
6.         my ($head) = /([^.]+)/;
   
7.         if ( ! $hash{$head} ) {
   
8.                 $hash{$head} = $_[-3];
   
9.                 $n=1;
   
10.                 $max{$head}{$n} = $_[-1];
    
11.                 $min{$head}{$n} = $_[-2];
    
12.         } else {
    
13.                 if ( $_[-2] >= $min{$head}{$n} && $_[-2] <= $max{$head}{$n} && $_[-1] > $max{$head}{$n} ) {
    
14.                         $max{$head}{$n} = $_[-1];
    
15.                 } elsif ( $_[-2] < $min{$head}{$n} && $_[-1] >= $min{$head}{$n} && $_[-1] <= $max{$head}{$n} ) {
    
16.                         $min{$head}{$n} = $_[-2];
    
17.                 } elsif ( $_[-2] >= $min{$head}{$n} && $_[-1] <= $max{$head}{$n} ) {
    
18.                 } else {
    
19.                         ++$n;
    
20.                         $max{$head}{$n} = $_[-1];
    
21.                         $min{$head}{$n} = $_[-2];
    
22.                 }
    
23. 
    
24.         }
    
25. }
    
26. for my $key ( sort keys %hash ) {
    
27.         for ( keys %{$max{$key}} ) {
    
28.                 say join "\t", $key, "$hash{$key}=$_", $min{$key}{$_}, $max{$key}{$_};
    
29.         }
    
30. }
    
31. __DATA__
    
32. a.1        1        3        b        2        8
    
33. a.2        2        5        b        5        9
    
34. a.3        3        6        b        4        7
    
35. a.4        2        4        b        12        15
    
36. c.1        3        8        d        1        7
    
37. c.2        4        8        d        6        8

復(fù)制代碼

1. [root@localhost ~]# ./p
   
2. a       b=1     2       9
   
3. a       b=2     12      15
   
4. c       d=1     1       8
   
5. [root@localhost ~]#

復(fù)制代碼

清泉一邊1

能改成這種形式嗎？
本人菜鳥，看不太懂啊
理想輸出
a        1_b        2-9_12-15
c        1_d        1        8
另外還需要輸出另一個文件，即最小交集（如果坐標(biāo)區(qū)間有重疊,則合并之并取其最大的起始坐標(biāo)和最小的終止坐標(biāo)）
理想結(jié)果是
a        1_b        5-7_12-15
c        1_d        6        7回復(fù) 6# yestreenstars


   

pitonas

{:2_172:}小伙伴們, 這個文件有沒有?

a.1        1        3        b        2        8
a.2        2        5        b        5        9
a.3        3        6        b        4        7
a.4        2        4        b        12        15
a.5        3        6        c        4        7
a.6        2        4        c        10        13
c.1        3        8        d        1        7
c.2        4        8        d        6        8
c.3        4        8        e        8        18
c.4        4        8        f        9        14

清泉一邊1

這是我做的一個簡單模型，實際數(shù)據(jù)比這個復(fù)雜些。回復(fù) 8# pitonas


   

pitonas

小伙伴們, 有沒有? {:2_171:}
回復(fù) 9# 清泉一邊1