【线性代数 | C++】行列式按行（列）展开

一、定义（余子式、代数余子式）

• 在$n$阶行列式中，把$(i,j)$元$a_{ij}$所在的第$i$行和第$j$列划去后，留下来的$n-1$阶行列式叫做$(i,j)$元$a_{ij}$的余子式，记作$M_{ij}$.
• 记$A_{ij}=(-1{)}^{i+j}{M}_{ij}，$$A_{ij}$叫做$(i,j)$元的代数余子式.

• 例如四阶行列式$$D=\left|\begin{array}{cccc}{a}_{11}& a_{12}& a_{13}& a_{14}\\ a_{21}& a_{22}& a_{23}& a_{24}\\ a_{31}& a_{32}& a_{33}& a_{34}\\ a_{41}& a_{42}& a_{43}& a_{44}\end{array}\right|$$中$(3,2)$元的$a_{32}$的余子式和代数余子式分别为$$M_{32}=\left|\begin{array}{ccc}{a}_{11}& a_{13}& a_{14}\\ a_{21}& a_{23}& a_{24}\\ a_{41}& a_{43}& a_{44}\end{array}\right|，$$ $$A_{32}=(-1{)}^{3+2}{M}_{32}=-M_{32}.$$

二、引理

• 一个$n$阶行列式，如果其中第$i$行所有元素除$(i,j)$元$a_{ij}$外都为零，那么这行列式等于$a_{ij}$与它的代数余子式的乘积，即$D=a_{ij}{A}_{ij}$.

三、定理3（行列式按行（列）展开法则）

• 行列式等于它的任一行（列）的各元素与其对应的代数余子式乘积之和，即$$D=a_{i1}{A}_{i1}+a_{i2}{A}_{i2}+\cdots +a_{in}{A}_{in}(i=1,2,\cdots ,n)，$$或$$D=a_{1j}{A}_{1j}+a_{2j}{A}_{2j}+\cdots +a_{nj}{A}_{nj}(j=1,2,\cdots ,n).$$

四、推论

• 行列式某一行（列）的元素与另一行（列）的对应元素的代数余子式乘积之和等于零.即$$a_{i1}{A}_{j1}+a_{i2}{A}_{j2}+\cdots +a_{in}{A}_{jn}=0,i\ne j,$$或$$a_{1i}{A}_{1j}+a_{2i}{A}_{2j}+\cdots +a_{ni}{A}_{nj}=0,i\ne j$$.

五、定理3的C++代码

#include <iostream>
#include <vector>
#include <cmath>

using namespace std;

//根据定理3求行列式值的函数
int GetDetValue(const vector<vector<int>> &vvDet)
{
    int iDetVal = 0;
    
    //当行列式内没有元素时，返回0
    if(0 == vvDet.size())
    {
        return 0;
    }
    //行列式仅有1个元素时，返回该元素的值
    if(1 == vvDet.size() && 1 == vvDet[0].size())
    {
        return vvDet[0][0];
    }
    else//行列式元素多于1个时，将行列式按第1行展开
    {
        for(unsigned int i = 0; i < vvDet[0].size(); i++)
        {
            //求余子式
            vector<vector<int>> vvCofactor(vvDet);
            vvCofactor.erase(vvCofactor.cbegin());//删除(0,i)元素所在行，即第1行
            for(unsigned int j = 0; j < vvCofactor.size(); j++) //删除（0，i）元素所在列
            {
                vvCofactor[j].erase(vvCofactor[j].cbegin() + i);
            }
            
            //将第1行元素与其代数余子式的乘积进行累加，即为行列式的值
            //行列式的阶数大于1时，进行递归
            iDetVal += (vvDet[0][i] * pow(-1, i) * GetDetValue(vvCofactor));
        }
    }
    return iDetVal;
}

//主函数
int main() 
{
    vector<vector<int>> vvDet{
                               {3,1,-1,2},
                               {-5,1,3,-4},
                               {2,0,1,-1},
                               {1,-5,3,-3}
                              };
    cout << "行列式的值 = " << GetDetValue(vvDet) << endl;
    return 0;
}

—-
引用文献：《工程数学 线性代数（第五版）》同济大学数学系编，高等教育出版社