计算动态矩阵行列式的(C++)函数

#include <iostream>
#include <vector>

// 计算矩阵行列式的函数
double calculateDeterminant(std::vector<std::vector<double>>& matrix) {
    int size = matrix.size();
    double determinant = 0.0;

    // 递归基本情况：当矩阵大小为1时，直接返回矩阵元素值
    if (size == 1) {
        return matrix[0][0];
    }

    // 递归情况：通过展开第一行计算行列式
    for (int i = 0; i < size; i++) {
        // 创建子矩阵
        std::vector<std::vector<double>> subMatrix(size - 1, std::vector<double>(size - 1));

        // 填充子矩阵
        for (int j = 1; j < size; j++) {
            for (int k = 0; k < size; k++) {
                if (k < i) {
                    subMatrix[j - 1][k] = matrix[j][k];
                } else if (k > i) {
                    subMatrix[j - 1][k - 1] = matrix[j][k];
                }
            }
        }

        // 递归计算子矩阵的行列式
        double subDeterminant = calculateDeterminant(subMatrix);

        // 根据符号和子矩阵行列式值更新总行列式值
        determinant += (i % 2 == 0 ? 1 : -1) * matrix[0][i] * subDeterminant;
    }

    return determinant;
}

int main() {
    // 示例用法
    std::vector<std::vector<double>> matrix = {{1, 2, 3}, {4, 5, 6}, {7, 8, 9}};
    double determinant = calculateDeterminant(matrix);
    std::cout << "矩阵行列式的值为: " << determinant << std::endl;

    return 0;
}