篩選所有連續(xù)塊

sequencing

現(xiàn)在有如下格式的文件，第一列是線段名稱，第二、三列是線段的起點和終點，線段之間沒有重疊，文件中線段是無序的，位置在下圖中

1. b        9        12
   
2. h        76        83
   
3. a        5        7
   
4. d        25        35
   
5. e        39        43
   
6. f        60        63
   
7. j        97        99
   
8. i        84        87
   
9. c        15        19
   
10. k        102        107
    
11. l        109        120
    
12. g        67        72

復(fù)制代碼

想實現(xiàn)下面的目的：
如果連續(xù)的兩條線段之間距離小于6，就說他們相鄰，
如果相鄰線段數(shù)大于2就說這些線段形成連續(xù)塊，
希望找出所有連續(xù)的塊。圖中a-b-c、f-g-h-i和j-k-l是三個滿足條件的情況
請高手幫忙！

sequencing

1. a-b-c
   
2. f-g-h-i
   
3. j-k-l

復(fù)制代碼

這樣就行了

yinyuemi

回復(fù) 3# sequencing

1. echo '    b        9        12
   
2.     h        76        83
   
3.     a        5        7
   
4.     d        25        35
   
5.     e        39        43
   
6.     f        60        63
   
7.     j        97        99
   
8.     i        84        87
   
9.     c        15        19
   
10.     k        102        107
    
11.     l        109        120
    
12.     g        67        72' |awk '{b[$2+$3]=$1;s[$1]=$2;e[$1]=$3}END{for(i=1;i<=asort(b);i++){if(i<length(b)){if(s[b[i+1]]-e[b[i]]<6){t=(i==1?b[i]:t)"-"b[i+1]}else{printf gsub("-","-",t)>=2?t"\n":"";t=b[i+1]}}};printf gsub("-","-",t)>=2?t"\n":""}'
    
13. 
    
14. a-b-c
    
15. f-g-h-i
    
16. j-k-l
    

復(fù)制代碼

改正：

1. awk '{a[$1]=$2+$3;b[$2+$3]=$1;s[$1]=$2;e[$1]=$3}END{for(i=1;i<=asort(a);i++){print i,b[a[i]];if(i<length(b)){if(s[b[a[i+1]]]-e[b[a[i]]]<6){t=(i==1?b[a[i]]:t)"-"b[a[i+1]]}else{printf gsub("-","-",t)>=2?t"\n":"";t=b[a[i+1]]}}};printf gsub("-","-",t)>=2?t"\n":""}'
   

復(fù)制代碼

sequencing

回復(fù) 4# yinyuemi


    yinyuemi，你簡直就是神！

yinyuemi

回復(fù) 5# sequencing


    其實,按照你定的規(guī)則一步一步來就行

cjaizss

現(xiàn)在有如下格式的文件，第一列是線段名稱，第二、三列是線段的起點和終點，線段之間沒有重疊，文件中線段是 ...
sequencing 發(fā)表于 2011-11-22 14:51

    沒有重疊的線段對嗎?

1. 
   
2. awk 'NF==3{name[$2]=$1;left[$1]=strtonum($2);right[$1]=strtonum($3)}END{
   
3. x=asort(left);
   
4. asort(right);
   
5. y=1;
   
6. s=name[left[1]];
   
7. for(i=2;i<=x;i++) {
   
8.      if(left[i]-right[i-1]>=6) {
   
9.             if(i-1>y+1)
   
10.                    print s;
    
11.             y=i;
    
12.             s=name[left[i]];
    
13.      } else {
    
14.             s=s"-"name[left[i]];
    
15.      }
    
16. }
    
17. if(y<x-1)
    
18.      print s;
    
19. }'
    

復(fù)制代碼

cjaizss

回復(fù)  sequencing
yinyuemi 發(fā)表于 2011-11-22 15:32

我覺得還是按長度來排序吧,不按名字吧

sequencing

沒有重疊的線段對嗎?
cjaizss 發(fā)表于 2011-11-22 15:57

    對，沒有重疊的線段
版主給的方法也是對的，消化中。。。