求行列式的值

如何用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;
}

可以看出,这种方法非常愚蠢暴力,那么如何优化这个思路呢,我们以后在讲

测试

测试样例

  • 1
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值