ええと、これは私に本当の頭痛を与えています。NxN 行列式を計算するための行列式関数を作成しており、再帰を使用しています。ロジックは正しく機能していますが、最終的な値を正しく計算できません。
マトリックス行列式のコードは次のとおりです。
public static double determinant(double[,]array){
double det=0;
double total = 0;
double[,] tempArr = new double[array.GetLength(0) - 1, array.GetLength(1) - 1];
if(array.GetLength(0)==2)
{
det = array[0, 0] * array[1, 1] - array[0, 1] * array[1, 0];
}
else {
for (int i = 0; i <1; i++)
{
for (int j = 0; j < array.GetLength(1); j++)
{
if (j % 2 != 0) array[i, j] = array[i, j] * -1;
tempArr= fillNewArr(array, i, j);
det+=determinant(tempArr);
total =total + (det * array[i, j]);
}
}
}
return det;
}
また、fillNewArr メソッドについては、配列をトリムするための単なるメソッドであり、メソッドは次のとおりです。
ublic static double[,] fillNewArr(double[,] originalArr, int row, int col)
{
double[,] tempArray = new double[originalArr.GetLength(0) - 1, originalArr.GetLength(1) - 1];
for (int i = 0, newRow = 0; i < originalArr.GetLength(0); i++)
{
if (i == row)
continue;
for (int j = 0, newCol=0; j < originalArr.GetLength(1); j++)
{
if ( j == col) continue;
tempArray[newRow, newCol] = originalArr[i, j];
newCol++;
}
newRow++;
}
return tempArray;
}
メソッドは「想定どおり」に機能していますが、最終結果が正しい方法で計算されていません。なぜでしょうか?!
4x4 配列の例:
{2 6 6 2}
{2 7 3 6}
{1 5 0 1}
{3 7 0 7}
最終結果は -168 になるはずですが、私の場合は 104 です!