行列式求值

最新推荐文章于 2025-01-17 10:22:56 发布

mr_zys

最新推荐文章于 2025-01-17 10:22:56 发布

阅读量3.3k

点赞数

文章标签：线性代数行列式

本文链接：https://blog.youkuaiyun.com/mr_zys/article/details/16952749

版权

博客内容介绍了如何求解线性代数中的行列式，提供了具体的题目链接和求解代码，帮助理解行列式的计算方法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

简单的行列式求值的模版

题目：http://acm.njupt.edu.cn/acmhome/problemdetail.do?&method=showdetail&id=1035

下面是代码：

#include <stdio.h>
#define maxn 11
#define e 1e-6
double c[maxn][maxn];

bool zero(double a)
{
  return a > -e && a < e;
}

double det(int n)
{
  double q,ans = 1;
  int i,j,k;
  for(i = 0; i < n; i++)
  {
    if(zero(c[i][i]))
    {
      for(j = i + 1; j < n; j++)
      {
	if(!zero(c[j][i])) break;
      }
      if(j == k) return 0;
      ans *= -1;
      for(k = i; k < n; k++)
      {
	double temp = c[j][k];
	c[j][k] = c[i][k];
	c[i][k] = temp;
      }
    }
    ans *= c[i][i];
    for(j = i + 1; j < n; j++)
    {
      if(!zero(c[j][i]))
      {
	q = c[j][i] / c[i][i];
	for(k = i; k < n; k++) c[j][k] -= c[i][k] * q;
      }
    }
  }
    return ans;
}
int main()
{
  int n;
  while(~scanf("%d",&n))
  {
    for(int i = 0; i < n; i++)
    {
      for(int j = 0; j < n; j++) scanf("%lf",&c[i][j]);
    }
    printf("%0.2lf\n",det(n));
  }
  return 0;
}