用perl除去有交集的行

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

复制代码