行列式求解【线性代数】

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cmath>
#include<map>
#include<string>
#include<functional>
using namespace std;
const int maxn = 12;
struct Dmt {
    int x;
    int y;
    int a[maxn][maxn];
};
int determinant_1(Dmt tmp)
{
    int a_1 = tmp.a[1][1];
    return a_1;
}
int determinant_2(Dmt tmp)
{
    int a_1 = tmp.a[1][1] * tmp.a[2][2];
    int a_2 = tmp.a[1][2] * tmp.a[2][1];
    return a_1 - a_2;
}
int determinant_3(Dmt tmp)
{
    int a_1 = tmp.a[1][1] * tmp.a[2][2] * tmp.a[3][3];
    int a_2 = tmp.a[1][2] * tmp.a[2][3] * tmp.a[3][1];
    int a_3 = tmp.a[1][3] * tmp.a[2][1] * tmp.a[3][2];
    int a_4 = tmp.a[1][3] * tmp.a[2][2] * tmp.a[3][1];
    int a_5 = tmp.a[1][1] * tmp.a[2][3] * tmp.a[3][2];
    int a_6 = tmp.a[1][2] * tmp.a[2][1] * tmp.a[3][3];
    return a_1 + a_2 + a_3 - a_4 - a_5 - a_6;
}
int get_determinant(Dmt tmp,int ans)//ans——当前行列式的值
{
    Dmt newtmp;
    if (tmp.x == 1 && tmp.y == 1)
        return determinant_1(tmp);
    if (tmp.x == 2 && tmp.y == 2)
        return determinant_2(tmp);
    if (tmp.x == 3 && tmp.y == 3)
        return determinant_3(tmp);
    //按照最后一行的代数余子式展开~
    newtmp.x = tmp.x - 1;
    //cout << newtmp.x << endl;
    newtmp.y = tmp.y - 1;
    //cout << newtmp.y << endl;
    int cnt = ans;
    for (int i = 1; i <= tmp.x; i++)//按照最后一行展开成代数余子式!
    {
        for (int j = 1; j <= tmp.x-1; j++)
        {
            int cs = 1;
            for (int k = 1; k <= tmp.y; k++)
            {
                if (k == i)
                    continue;
                newtmp.a[j][cs] = tmp.a[j][k]; 
                cs++;
            }

        }
        cnt+=pow(-1,tmp.y+i)*tmp.a[tmp.y][i]*get_determinant(newtmp,0);
    }
    return cnt;
}
int main()
{
    Dmt test;
    int x; int y;
    while (cin >> test.x >> test.y)
    {
        for (int i = 1; i <= test.x; i++)
        {
            for (int j = 1; j <= test.y; j++)
            {
                cin >> test.a[i][j];
            }
        }
        cout<<get_determinant(test, 0)<<endl;
    }
    system("pause");
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值