perl区间交集问题（急！）

清泉一边1

第2列三列没用,那么用每一行的第5列和第6列与下面各列的进行比较，5列为起始坐标，6列为终止坐标。输出他们的最大交集（即如果坐标区间有重叠,则合并之并取其最小的起始坐标和最大的终止坐标）
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出现的次数
另外还需要输出另一个文件，即最小交集（如果坐标区间有重叠,则合并之并取其最大的起始坐标和最小的终止坐标）
理想结果是
a        b=1        5        7
a        b=2        12        15
c        d        6        7
求大神帮忙！！！

jjboy110

且不另你年少无知，当你得不到答案的时候多了，你才能知道其实提问是有提问的智慧的，问题最终都是自己解决，如果要别人帮你解决，请花钱。

清泉一边1

回复 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

复制代码

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

复制代码

结果和你的有些出入，不清楚你说的b出现的次数具体指什么，而且格式你需要自己调整下。

清泉一边1

太感谢了！！！回复 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

复制代码

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 ~]#

复制代码

清泉一边1

能改成这种形式吗？
本人菜鸟，看不太懂啊
理想输出
a        1_b        2-9_12-15
c        1_d        1        8
另外还需要输出另一个文件，即最小交集（如果坐标区间有重叠,则合并之并取其最大的起始坐标和最小的终止坐标）
理想结果是
a        1_b        5-7_12-15
c        1_d        6        7回复 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

这是我做的一个简单模型，实际数据比这个复杂些。回复 8# pitonas


   

pitonas

小伙伴们, 有没有? {:2_171:}
回复 9# 清泉一边1