如何用C语言求行列式子的值呢值呢?
肯定要用递归,肯定要用数组
废话少说,上代码
首先定义一个函数,求行列式去掉第一行某个元素的余子式
用一个二维数组来表示这个矩阵,但是参数传递的时候我们不能确定这个矩阵的行数列数,所以编译器会报错
所以定义一个一维数组来模拟这个行列式,而一维数组作为参数传递给函数时又不用指定大小,符合我们的需要。
同时我们定义一个结构体来表示这个余子式
typedef struct yuzi{
int *a;//表示指向余子式的数组
int n;//表示余子式的长宽
}yuzi;
yuzi *yu(int a[],int k,int j,int s){//a是要计算余子式的行列式,k是他的宽度,j是要求第一行哪个位置的余子式,s默认是零(表示求第一行的余子式)
int *b=new int[(k-1)*(k-1)];
int l=0;
for(int i=0;i<k;i++)
for(int m=0;m<k;m++){
if(m!=j&&i!=s)
b[l++]=a[k*i+m];//b[k*i+m]就相当于行列式的i行m列;
}
yuzi *c=new yuzi;
c->a=new int[(k-1)*(k-1)];
for(int i=0;i<k-1;i++)
for(int j=0;j<k-1;j++){
c->a[i*(k-1)+j]=b[i*(k-1)+j];
}
c->n=k-1;
return c;
}
然后定义一个递归函数,求一个行列式的值
思路:如果这个行列式 a[] 的宽度为二,直接返回a[0]*a[3]-a[1]*a[2];
否则 返回他的代数余子式的和
int calcu(int a[],int k)
{
if(k==2){
return(a[0]*a[3]-a[1]*a[2]);
}
else{
int sum=0;
for(int i=0;i<k;i++){
yuzi *y=yu(a,k,i,0);
int k=pow(-1,i)*calcu(y->a,y->n)*a[i];
sum+=k;
}
return sum;
}
}
然后是我们的main函数
int main(){
int n;
cin>>n;
int b;
int *a=new int[n*n];//动态申请,用多少,申请多少
for(int i=0;i<n;i++)
for(int j=0;j<n;j++){
cin>>a[n*i+j];
}
int c=calcu(a,n);
cout<<"最终答案"<<c<<endl;
return 0;
}
可以看出,这种方法非常愚蠢暴力,那么如何优化这个思路呢,我们以后在讲