匹配所有字母前的数字运算--大工程~

huang6894

回复 38# damcool


    牛人！

huang6894

回复 39# yestreenstars


    谢谢您！

damcool

huang6894 发表于 2013-11-15 18:42
回复 38# damcool

希望不是为了研究转基因东东来害我们。。。。

yestreenstars

用perl琢磨了一下，我对perl不是很熟，只能写成这样了，第一题：

1. #!/usr/bin/perl
   
2. use 5.010;
   
3. open(FILE,"use.txt") or die;
   
4. $len = 10;
   
5. while (<FILE>) {
   
6.         split;
   
7.         $sum = 0;
   
8.         $N = $_[8] > 0 ? 1 : 2;
   
9.         while ($_[5] =~ /(\d+)S/g) {
   
10.                 push(@numbers,$1);
    
11.         }
    
12.         $length = @numbers;
    
13.         if ($length > 1) {
    
14.                 if ($numbers[0] >= $len) {
    
15.                         printf "%s_%s\n%s\n+\n%s\n\n",$_[0],$N,substr($_[9],0,$numbers[0]),substr($_[10],0,$numbers[0]);
    
16.                 }
    
17.                 if ($numbers[1] >= $len) {
    
18.                         printf "%s_%s\n%s\n+\n%s\n\n",$_[0],$N,substr($_[9],length($_[9])-$numbers[1]),substr($_[10],length($_[10])-$numbers[1]);
    
19.                 }
    
20.         } else {
    
21.                 if ($numbers[0] >= $len) {
    
22.                         if ($_[5] =~ /S$/) {
    
23.                                 while ($_[5] =~ /(\d+)[MI]/g) {
    
24.                                         $sum += $1;
    
25.                                 }
    
26.                         }
    
27.                         printf "%s_%s\n%s\n+\n%s\n\n",$_[0],$N,substr($_[9],$sum,$numbers[0]),substr($_[10],$sum,$numbers[0]);
    
28.                 }
    
29.         }
    
30.         undef @numbers;
    
31. }

复制代码

yestreenstars

优化了一下第一题的awk：

1. #!/usr/bin/awk
   
2. {
   
3.         s=0;t=$6;N=$9>0?1:2;
   
4.         gsub(/[0-9]*[DMI]/,"",t);
   
5.         l=split(t,a,"S");
   
6.         if(l>2){
   
7.                 if(a[1]>=10)printf "%s_%s\n%s\n+\n%s\n\n",$1,N,substr($10,1,a[1]),substr($11,1,a[1]);
   
8.                 if(a[2]>=10)printf "%s_%s\n%s\n+\n%s\n\n",$1,N,substr($10,length($10)-a[2]+1),substr($11,length($11)-a[2]+1)
   
9.         } else {
   
10.                 if(a[1]>=10){
    
11.                         gsub(/[0-9]*(D|S.*)/,"",$6);
    
12.                         split($6,b,"[MI]");
    
13.                         for(i in b)s+=b[i];
    
14.                         printf "%s_%s\n%s\n+\n%s\n\n",$1,N,substr($10,s+1,a[1]),substr($11,s+1,a[1])
    
15.                 }
    
16.         }
    
17. }

复制代码

huang6894

@damcool谢谢您~@yestreenstars由衷感谢大大的帮助~

pitonas

damcool 发表于 2013-11-15 14:42
希望不是为了研究转基因东东来害我们。。。。

希望不是{:2_179:}