朴素的矩阵行列式计算之C++实现

来源:互联网 发布:联通网络客服人工电话 编辑:程序博客网 时间:2024/04/30 09:05

在求n行n列矩阵的行列式时,可以根据选择比较多0的行或者比较多0的列来划掉进行下一步(依稀还记得大一上线性代数的时候的一些知识点),

我这里是简单的暴力版本,也没有优化过的,每次只划掉当前矩阵的头行来进行递归计算:

int mat[100][100] = {{1,-1,3,2},{-2,-3,3,2},{6,-1,-1,2},{-4,-1,3,-1}};int det(vector<vector<int> > vec, int n){if(n == 1) return vec[0][0];//递归终止条件else{int sum = 0 ;for(int j = 0 ; j < n; j++){//得到从当前矩阵中划去第0行和第j列的所有元素后得到的矩阵vector<vector<int> > s;for(int r = 1 ; r < n ; r++){vector<int> temp;for(int c = 0 ; c < n ; c++){if(c != j) temp.push_back(vec[r][c]);}s.push_back(temp);} sum += pow(-1,0-j) * vec[0][j] * det(s,n-1);}return sum;} }int main(){int n; //方阵的行数,列数cin>>n;vector<vector<int> > vec; //放入vector中,方便处理 for(int i = 0 ; i < n ; i++){vector<int> temp;for(int j = 0 ; j < n; j++)temp.push_back(mat[i][j]);vec.push_back(temp);}printf("矩阵的行列式是%d\n",det(vec,n));}


原创粉丝点击