行列式求值

最新推荐文章于 2022-09-17 18:32:59 发布

MrTinTin

最新推荐文章于 2022-09-17 18:32:59 发布

阅读量393

点赞数

分类专栏：数学

本文链接：https://blog.csdn.net/MrTinTin/article/details/77652084

版权

数学专栏收录该内容

28 篇文章 0 订阅

订阅专栏

O(n^3)的算法，利用了行列式的性质。

#include<cstdio>
#include<iostream>
#define N 19
using namespace std;
double a[N][N];
int n;
double calcudet()
{
	int out;
	int i,j;
	double k,result=1;
    for(out=0;out<n;out++){                                   //外层阶数循环 
    	for(i=out;i<n;i++){                                   //寻找第一项不为零的行
    		if(a[i][out]==0) continue;
    	    for(j=out;j<n;j++){                            //第一项不为零的行与顶行交换 
    	    	k=a[i][j];a[i][j]=a[out][j];a[out][j]=k;
    	    }
    	    if(i==out){
    	    	result*=a[out][out]; 
    	    }
    	    else{
				result*=-a[out][out]; 
			}                   //结果变号 
    	    break;
		}
    	if(i==n) return 0;
    	
    	if(a[out][out]!=1)
    		for(j=out+1;j<n;j++)                                 //将第一项变为1 
    			a[out][j]/=a[out][out];
        a[out][out]=1;
    	
    	for(i=out+1;i<n;i++)                                     //将该列其余项消为0 
    		for(j=out+1;j<n;j++)
    			a[i][j]-=a[i][out]*a[out][j];
    }
    return result;
}
int main()
{
	int i,j;
	scanf("%d",&n);
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			scanf("%lf",&a[i][j]);
		}
	}
	cout<<calcudet();
	return 0;
}