C++计算行列式（函数）

C++计算行列式

思路：

1. 用二维数组存放行列式，然后将首地址传入函数，通过递归不断储存第一列的代数余子式，直到传入的是一个2*2的数组首地址
2. 目前阶段，我仅仅知道，在计算机中，二维数组名的本质类似一维常指针，a[i][j]等同于*(a+i*n+j)，n为a数组第二个方括号内的数，要想不限定传入函数的二维数组首地址指向的大小，可以将其强制转换成一维指针，然后通过指针操作

代码如下：

//det
double det(double *d,int n){
	double temp=0;
	if(2==n)
		temp= (*d) * (*(d+1*n+1)) -(*(d+1*n)) *(*(d+1));
	else if(1==n)temp=*d;
	else{
	//求出n个代数余子式
		double complement[n-1][n-1];
		for(int i=0;i<n;i++){//关于第一行的n列展开 
		
		//开辟了代数余子式空间
			for(int k=0;k<n-1;k++){
				for(int j=0;j<n-1;j++){
					complement[k][j]=(*(d+(k+1)*n+(j<i?j:j+1)));
				}
			}
			if(*(d+i)!=0)temp+=*(d+i)*det((double*)complement,n-1)*(((i+2)%2)==0?1:-1);
			//判断一下，再计算，提高效率 
		}
	}
	return temp;
}
int main(){
	double a[][3]={1,0,0,0,-2,0,0,0,3};
	cout<<det((double*)a,3);
	return 0;
} 

缺陷

1. 每增加一阶，由于递归的存在，时长成n倍增长，12阶相当于12个11阶行列式计算的时长，当阶数超过12，由于计算时间过长，没有必要用这个函数
2. 每次调用这个函数，都要先把二维数组名强制转化为一维double指针