行列式求值

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

靖-Drei

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

阅读量876

点赞数

分类专栏：南邮acm

本文链接：https://blog.csdn.net/sundujing/article/details/46939519

版权

南邮acm 专栏收录该内容

109 篇文章 8 订阅

订阅专栏

Description

求 N 阶方阵对应的行列式值。

Input

输入数据包含多组测试案例。每组测试案例由N（0

#include<stdio.h>
#define MAXN 11

#define fabs(x) ((x)>0?(x):-(x))
#define zero(x) (fabs(x)<1e-10)

struct mat{
//    int n,m;
    int n;
    double data[MAXN][MAXN];
};
double det(const mat& a){
    int i,j,k,sign=0;
    double b[MAXN][MAXN],ret=1,t;
//    if (a.n!=a.m)
//        return 0;
    for (i=0;i<a.n;i++)
        for (j=0;j<a.n;j++)
            b[i][j]=a.data[i][j];
    for (i=0;i<a.n;i++){
        if (zero(b[i][i])){
            for (j=i+1;j<a.n;j++)
                if (!zero(b[j][i]))
                    break;
            if (j==a.n)
                return 0;
            for (k=i;k<a.n;k++)
                t=b[i][k],b[i][k]=b[j][k],b[j][k]=t;
            sign++;
        }
        ret*=b[i][i];
        for (k=i+1;k<a.n;k++)
            b[i][k]/=b[i][i];
        for (j=i+1;j<a.n;j++)
            for (k=i+1;k<a.n;k++)
                b[j][k]-=b[j][i]*b[i][k];
    }
    if (sign&1)
        ret=-ret;
    return ret;
}


int main()
{
    mat a;
    while(scanf("%d",&a.n)!=EOF)
    {
        for(int i=0;i<a.n;i++)
            for(int j=0;j<a.n;j++)
                scanf("%lf",&a.data[i][j]);
        double det1;
        det1=det(a);
        printf("%.2f\n",det1);
    }
    return 0;
}