列表计算问题

newfinder

各位大侠，小弟现在有一个对列表作计算的问题，不知道怎么写比较好

现在有这样一个列表data，如下：
O B M
1 - -
2 - -
3 1 -
4 1 2
5 3 4
6 1 4
7 5 6

想计算一下这些数字所指代的事物之间的内部关系，
具体的计算方法是：

对角线数字的计算方法：
1、当Bi或Mi不存在时，A(ii)=1，如A(11)=1, A(22)=1, A(33)=1
2、当Bi和Mi都存在时，A(ii)=1+0.5A(BiMi),如A(44)=1+0.5A(B4M4)=1+0.5A(12)=1+0=1,   A(55)=1+0.5A(B5M5)=1+0.5A(34)=1+0.5*0.25=1.125

非对角线数字的计算方法：
Aij=Aji
1、当Bj和Mj都不存在时，A(ij)=0,如A(12)=A(21)=0，此时的B2和M2都不存在
2、当Bj存在，Mj不存在，A(ij)=0.5A(iBj),如A(13)=A(31)=0.5A(1B3)=0.5A(11)=0.5*0.5=0.25,  A(23)=A(32)=0.5A(2B3)=0.5A(21)=0
3、当Bj不存在，Mj存在，A(ij)=0.5A(iMj)，计算方法同上类似
4、当Bj和Mj都存在时，A(ij)=0.5(A(iBj)+A(iMj)),如A(34)=0.5(A(3B4)+A(3M4))=0.5(A(31)+A(32))=0.5*0.5=0.25，此时B4和M4都存在

故这样计算的结果如下：
1 0 0.5 0.5  0.5   0.75   0.625
   1 0    0.5  0.25  0.25   0.25
      1   0.25 0.625 0.375  0.5
             1    0.625 0.75   0.6875
                   1.125 0.5625 0.84375
    对称                    1.25   0.90625
                                      1.28125

不知道这个程序应该怎么写比较好些呢？请大家支支招

huang6894

看晕了~坐等大神回复吧。。。。

我的理解是，建立一个二维数组，然后每读入一行存到hash（以位置为key），然后得到所有的组合后替换对应的数组元素~

pitonas

输出结果:

1. 1       0       0.5     0.5     0.5     0.75    0.625
   
2. 0       1       0       0.5     0.25    0.25    0.25
   
3. 0.5     0       1       0.25    0.625   0.375   0.5
   
4. 0.5     0.5     0.25    1       0.625   0.75    0.6875
   
5. 0.5     0.25    0.625   0.625   1.125   0.5625  0.84375
   
6. 0.75    0.25    0.375   0.75    0.5625  1.25    0.90625
   
7. 0.625   0.25    0.5     0.6875  0.84375 0.90625 1.28125
   

复制代码

如果所有的 B, M 比行号小 ~ {:2_172:}


回复 1# newfinder


   

pitonas

我的理解 ~ {:2_172:}

1. #!/usr/bin/perl
   
2. use 5.010;
   
3. my @A;
   
4. my @B = map [ @_ = split, ~~ grep /-/, @_[ 1, 2, 2 ] ], <DATA>;
   
5. 
   
6. for my $i ( 1 .. $#B ) {
   
7.     $A[$i][$i] = SA($i, $i);
   
8.     $A[$_][$i] = $A[$i][$_] = SA( $i, $_ ) for $i + 1 .. $#B;
   
9. }
   
10. 
    
11. say join "\t", splice @$_, 1 for splice @A, 1;
    
12. 
    
13. sub SA {
    
14.     my ( $i, $j, $k ) = ( @_, $B[ $_[-1] ][-1] );
    
15.     $i == $j
    
16.     ? ( $k ? 1 : 1 + $A[ $B[$i][1] ][ $B[$i][2] ] * 0.5 )
    
17.     : ( $k == 0 ? $A[$i][ $B[$j][1] ] + $A[$i][ $B[$j][2] ] :
    
18.         $k == 1 ? $A[$i][ $B[$j][2] ] :
    
19.         $k == 2 ? $A[$i][ $B[$j][1] ] : 0 ) * 0.5
    
20. }
    

复制代码

huang6894

回复 4# pitonas


    niubility~~~~~~

newfinder

回复 2# huang6894

……坐等
   

newfinder

回复 4# pitonas

好赞好赞！！
是这样子的，大侠威武啊！