用perl除去有交集的行

清泉一边

           start        end
a        1        6
a        2        7
a        15        19
a        3        9
a        2        13
a        21        25
b        2        8
b        32        7
b        5        11
b        13        17
b        19        22
..
..
我想输出不含交集的行，即只要数据中行与行的start——end之间有交集，就舍去那些行。想得到如下结果：
a        15        19
a        21        25
b        13        17
b        18        22
..
..
等大神的代码！看能不能尽量省内存~
万谢！！！

清泉一边

为什么没有给解决的大神啊！！！？ 回复 1# 清泉一边


   

Monox

要省内存的话代码可以参考http://bbs.chinaunix.net/thread-4119878-1-2.html这个帖子里我的回复（55楼：http://bbs.chinaunix.net/forum.p ... &fromuid=663671），当然他的需求跟你的不完全一样，你根据我的代码稍稍修改一下就可以了。你要是说对Perl不熟啊什么的（这个论坛很多人这样说），就等其他好心人帮你修改我的代码或者提供新的代码。

ba_du_co

有行 start > end
b        32        7

start < end 的话可以参考

1. my ( $key, @data ) = '#!/usr/bin/perl';
   
2. while (<DATA>) {
   
3.     my ( $head, @tail ) = split;
   
4.     &Print if @data and ( $head ne $key );
   
5.     ( $key, $data[@data] ) = ( $head, [@tail] );
   
6. }
   
7. 
   
8. &Print;
   
9. 
   
10. sub Print {
    
11.     @data = sort { $a->[0] <=> $b->[0] } @data;
    
12.     my @print = shift @data;
    
13.     for my $r ( splice @data ) {
    
14.         $r->[0] < $print[-1][-1]
    
15.           ? push @{ $print[-1] },
    
16.               ( $r->[1], $print[-1][-1] )[ $r->[1] < $print[-1][-1] ]
    
17.           : @{ $print[-1] } > 2 ? ( $print[-1] = $r ) : push @print, $r;
    
18.     }
    
19.     print join( "\t", $key, @$_ ), "\n" for @print;
    
20. }
    
21. 
    
22. __DATA__
    
23. a        1        6
    
24. a        2        7
    
25. a        15        19
    
26. a        3        9
    
27. a        2        13
    
28. a        21        25
    
29. b        2        8
    
30. b        32        37
    
31. b        5        11
    
32. b        13        17
    
33. b        19        22

复制代码

清泉一边

谢谢！！！回复 4# ba_du_co


   

expert1

这应该是perl版的月经问题了，最少我看过N次

expert1

回复 1# 清泉一边

1. [root@localhost tmp]# cat a
   
2. a        1        6
   
3. a        2        7
   
4. a        15        19
   
5. b        7        32
   
6. a        3        9
   
7. a        21        25
   
8. b        2        8
   
9. b        5        11
   
10. a        2        13
    
11. b        13        17
    
12. b        19        22
    
13. c        3         5
    
14. c        5         7
    
15. c        200       400
    
16. a        100       120
    
17. d        5         9
    
18. d        3         5
    
19. m        2         5
    
20. 
    
21. awk 'NR == FNR{ a[$2" "$3]=$1 }NR>FNR{for(i in a){ flag=0
    
22. >
    
23. >                if(a[i]==$1){  split(i , m ," ")
    
24. >
    
25. >                               if ( ($2 >  m[1] && $2 <  m[2]) || ($3 > m[1] && $3 < m[2]) || ($2 <  m[1] && $3 >  m[2]) )
    
26. >
    
27. >                                     flag = 1     }
    
28. >
    
29. >                    if (flag) delete a[i]
    
30. >                  }}END { for(i in a)print a[i]" "i}' a a
    
31. a 100 120
    
32. c 200 400
    
33. a 15 19
    
34. m 2 5
    
35. d 3 5
    
36. a 21 25
    
37. c 5 7
    
38. d 5 9
    

复制代码

如果取决区间无重合，比如1 3和3 5/ 6 8 和8 10 /  7 9和 7 9 ，这个脚本可以解决你的问题。
当然前提是$2<$3,如果不是，也不是什么难题，加个条件就能过滤了。

expert1

awk毫无压力，写个区间相交的，比如1 3 和 3 5以及 8 9 和9 10这种。

该脚本前提是第二列<=第三列，若最有这种数据，请自行处理，是舍弃还是交换，比如32 7是舍弃了还是转换为7 32。

1. [root@221 ~]# cat f
   
2. a        1        6
   
3. a        2        7
   
4. a        15        19
   
5. b        27        32
   
6. a        3        9
   
7. a        21        25
   
8. b        2        8
   
9. b        5        11
   
10. a        2        13
    
11. b        13        17
    
12. b        19        22
    
13. c        3         5
    
14. c        5         7
    
15. c        200       400
    
16. a        100       120
    
17. d        5         9
    
18. d        3         5
    
19. b        15        20
    
20. d        500     600
    
21. m        2         5
    
22. [root@221 ~]# awk 'NR==FNR{a[$1" "$2" "$3]=$1;b[$1]++}NR>FNR{
    
23. >
    
24. >         counter=0
    
25. >      
    
26. >         for(i in a)
    
27. >
    
28. >         {
    
29. >
    
30. >                 if($1==a[i])
    
31. >
    
32. >                    {
    
33. >
    
34. >                      split(i,m," ")
    
35. >
    
36. >
    
37. >                      if($3<m[2] || $2>m[3] || (m[2]==$2 && m[3]==$3))
    
38. >
    
39. >
    
40. >                           counter++
    
41. >                                   
    
42. >                      else next
    
43. >                        
    
44. >                           
    
45. >                      if(counter == b[$1])
    
46. >
    
47. >                         print $0
    
48. >
    
49. >                      }
    
50. >
    
51. >         }
    
52. >
    
53. > }' f f
    
54. a        15        19
    
55. b        27        32
    
56. a        21        25
    
57. c        200       400
    
58. a        100       120
    
59. d        500     600
    
60. m        2         5
    
61. [root@221 ~]#

复制代码

思嘉不丽

猜猜我是谁

清泉一边

I don't know~回复 9# 思嘉不丽