行列式求值

来源：互联网发布：linux重新分区编辑：程序博客网时间：2024/04/30 15:48

简单的行列式求值的模版

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

下面是代码：

#include <stdio.h>#define maxn 11#define e 1e-6double 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;}