练习题

yestreenstars

原帖：http://bbs.chinaunix.net/forum.php?mod=viewthread&tid=4115260&fromuid=26669282

处理前：

1.    name1 name2 name3 name4 name5
   
2. MG1 A G C D C A C G A D
   
3. MG2 A G D C A C G C D A
   
4. SH_1 G G C C C C 0 0 A A
   
5. SH_2 A A 0 0 A A C C D D
   
6. SH_3 G A D D 0 0 G G 0 0
   
7. SD1 A G C D C C C G A D
   
8. SD2 0 0 D C C A C C D A

复制代码

处理后：

1.         MG      SH      SD
   
2. name1   2,2     3,3     1,1
   
3. name2   2,2     2,2     2,2
   
4. name3   2,2     2,2     1,3
   
5. name4   2,2     2,2     3,1
   
6. name5   2,2     2,2     2,2

复制代码

处理规则：
每个"name"分别对应两列，现在想统计以MG、SD、SH开头的行分别含有字母的个数，字母只可能是A、C、D、G这四种，每两列中只含有其中的两种，但不一定是哪两种，请按格式处理文本（注意顺序）。

下面是我的方法：

游客，如果您要查看本帖隐藏内容请回复

li0924

1. awk '{s=substr($1,1,2);b[s];for(i=2;i<=NF;i++){num=int(i/2);a[s,num,$i]++}}END{split("ACDG",c,"")

复制代码

后续的排列真的很难！

aqbssh

非得回复？

huang6894

我要学习学习~

yestreenstars

回复 3# aqbssh
我发的帖都是这样的，我希望能答题者能独立思考后再回复，不喜勿回~

   

reyleon

顺序是不是应该是这样子的呐？

1. [root@centos6-1 ~]# awk -f f.awk file
   
2.         MG      SH      SD      
   
3. name1   2,2     3,3     1,1     
   
4. name2   2,2     2,2     2,2     
   
5. name3   2,2     2,2     1,3     
   
6. name4   2,2     2,2     3,1     
   
7. name5   2,2     2,2     2,2     
   
8. [root@centos6-1 ~]#

复制代码

Herowinter

回复 1# yestreenstars
折腾了一个中午，终于勉强写出来了，用了
三维数组，第一维1-5分别表示name1-5，
第二维 1-3 表示 MG SD SH
第三维表示 A C D G
代码如下，请大大指教。

1. awk '{if($1~/MG/)for(i=2;i<=NF;i++)a[int(i/2),1,$i]++;else if($1~/SD/)for(i=2;i<=NF;i++)a[int(i/2),2,$i]++;else if($1~/SH/)for(i=2;i<=NF;i++)a[int(i/2),3,$i]++}END{c[1]="A";c[2]="C";c[3]="D";c[4]="G";print "\tMG\tSD\tSH";for(i=1;i<=5;i++){printf "name"i;for(j=1;j<=3;j++){flag=0;for(k=1;k<=4;k++)if(a[i,j,c[k]]>0){if(!flag)printf "\t%d",a[i,j,c[k]];else printf ",%d",a[i,j,c[k]];flag=1}}print ""}}' i
   
2. 
   
3.                 MG         SD         SH
   
4. name1        2,2        1,1        3,3
   
5. name2        2,2        2,2        2,2
   
6. name3        2,2        1,3        2,2
   
7. name4        2,2        3,1        2,2
   
8. name5        2,2        2,2        2,2
   

复制代码

reyleon

1. #!/bin/awk
   
2. 
   
3. BEGIN{ n=split("ACDG",c,"");for(i=0;++i<=n;)e[c[i]]=1 }
   
4. NR==1{ l=split($0,a,FS);next }
   
5. 
   
6. {
   
7.         t=substr($1,1,2)
   
8.         if(!m[t]++)q[++seq]=t
   
9.         num=2
   
10.         for(x=1;x<=l;x++){
    
11.                 for(y=num;y<num+2;y++)if($y in e)f[a[x],t,$y]++
    
12.                 num+=2
    
13.         }
    
14. }
    
15. 
    
16. END {
    
17.         printf("%-8s","");for(i=1;i<=seq;i++)printf("%-8s",q[i]);printf "\n"
    
18.         for(i=1;i<=l;i++){
    
19.                 printf("%-8s",a[i])
    
20.                 for(j=1;j<=seq;j++){
    
21.                         for(k=1;k<=n;k++)if(f[a[i],q[j],c[k]])s=s?s","f[a[i],q[j],c[k]]:f[a[i],q[j],c[k]]
    
22.                         printf("%-8s",s)
    
23.                         s=""
    
24.                 }
    
25.                 printf "\n"
    
26.         }
    
27. }

复制代码

yestreenstars

回复 6# reyleon
我复制结果时复制错了~已改~

   

yinyuemi

回复 1# yestreenstars

1. awk 'NR==1{for(i=1;i<=NF;i++)a[i]=$i;next}{sub(/_?[0-9]+$/,"",$1);if(!c[$1]++){printf "\t" $1;d[++m]=$1};for(i=2;i<=NF;i++){if($i~/[A-Z]/)b[a[int(i/2)]][$1][$i]++}}END{print "";for(i in b){printf i;for(j=1;j<=m;j++){printf "\t";t=0;for(k in b[i][d[j]])printf ((!t++?"":",")b[i][d[j]][k])};print ""}}'

复制代码