Matrix（1133）

今天看了算法如下：1.求一个二维数组i和j互质的情况，就是求i和j的公约数只能为1.采用求最大公约数的算法实现。

2.求向量的相乘，具体实现算法：需要判断两个数组长度是否一致。

3.求两个二维数组的乘积，主要是注意行等于列的情况，中间值k一定为第一个数组的列或者是第二个数组的行。

1.向量乘积：

result += a[i]*b[i];

具体代码：

public static double dot( double[] x , double[] y){ //向量点乘

double result = 0 ;

if(x. length == y. length){

for ( int i = 0 ; i < x. length ; i ++)

result += x[i] * y[i] ; //不用加括号，优先级是不一样的

} else {

StdOut. println ( "Error!" ) ;

return result ;

}

return result ;

}

2.向量相乘：MxN * NxP， M与P不一定相等，因为超出的部分为零。

其中设置共同的部分N为k值。

代码如下：

public static double[][] mult( double[][] a , double[][] b){ //矩阵与矩阵之积 [M][N] * [N][P] = [M][P]

if (a[ 0]. length == b. length){

double [][] result = new double [a. length ][b[ 0 ]. length ] ; //第一个是行，第二个是列

for ( int i = 0 ; i < result. length ; i ++)

for ( int j = 0 ; j < result[i]. length ; j ++)

for ( int k = 0 ; k < b. length ; k ++) //这里b.length 可以换为 a[0].length

result[i][j] += a[i][k] * b[k][j] ;

return result ;

} else {

StdOut. println ( "Error!!!" ) ;

return null;

}

}

3.转置矩阵：

就是从新设置一个数组，让他的大小满足转置后的条件，比如[i][j] 变为[j][i] .将原来数组的内容复制过去就可以了。

代码如下：

public static double[][] transpose( double[][] a){ //转置矩阵

double[][] result = new double[a[ 0]. length][a. length] ;

for( int i = 0 ; i < a. length ; i ++){

for ( int j = 0 ; j < a[ 0 ]. length ; j ++){

result[i][j] = a[j][i] ;

}

}

return result ;

}