用C语言编程递归求行列式的值

原创 2007年09月17日 10:43:00
/************************************************************************
功能: 用递归方法求行列式的值
作者: 黄志斌 于 广西河池
在自学《高等数学(二)》时,做行列式习题所得结果有些与教材所附
答案不一致,为了验证究竟答案的正确性,编写了这个程序。
其实用MS Office中Excel也可以验证,不过自己 DIY  感觉更舒服。
    希望这个程序能对您有所帮助。
日     期     概  要
-----------------------------------------------------
2001.3.23     完成运算功能
2002.3.17     按照匈牙利式命名规则和
              林锐博士的《高质量C/C++编程指南》
              增强程序可读性
2002.3.19     停止使用全局变量iRowOrCol,
              将原先的调试部分改为预处理
*********************************************************************/
#include "stdio.h"
#include "mem.h"
#define IsDebug 0

/*int iRowOrCol;*/
void ReadMatrix(int *ipMatrixRead, int iRowOrCol4Read)
{
 int iRow4Read, iCol4Read;
 for (iRow4Read = 0; iRow4Read < iRowOrCol4Read; iRow4Read ++)
 {
  for (iCol4Read = 0; iCol4Read < iRowOrCol4Read; iCol4Read ++)
  {
   printf("matrix[%d, %d]=", iRow4Read, iCol4Read);
   scanf("%d", (ipMatrixRead + iRow4Read * iRowOrCol4Read + iCol4Read));
  }
 }
}


void WriteMatrix(int *ipMatrixWrite, int p_iRowOrCol)
{
 int iRow4Write, iCol4Write;
 putch('┌');

 for (iRow4Write = 0; iRow4Write < (p_iRowOrCol + 1) * 3 - 1; iRow4Write ++)
 {
  putch(' ');
 }
 printf("%c/n", '┐');

 for (iRow4Write = 0; iRow4Write < p_iRowOrCol; iRow4Write ++)
 {
  putch('│');
  for (iCol4Write = 0; iCol4Write < p_iRowOrCol; iCol4Write ++)
  {
   printf("%3d", *(ipMatrixWrite + iRow4Write * p_iRowOrCol + iCol4Write));
  }
  printf("%c/n", '│');
 }

 putch('└');
 for (iRow4Write = 0; iRow4Write < (p_iRowOrCol + 1) * 3 -1; iRow4Write ++)
 {
  putch(' ');
 }
 printf("%c/n", '┘');

 #if IsDebug == 1
  getch();
 #endif
}


int ComputMatrix(int *ipMatrixComput, int iColComput)
{
 int iRow4Comput, iCol4Comput, iRowTmp, *ipNewMatrix, CurrentRow, iResult;

 if (1 == iColComput)
 {
  iResult = *ipMatrixComput;
 }
 else
 {
  iResult = 0;
  for (CurrentRow = 0; CurrentRow < iColComput; CurrentRow ++)
  {
   if (*(ipMatrixComput + CurrentRow * iColCom! put))
   {
    if (ipNewMatrix = (int *)malloc((iColComput - 1) * (iColComput - 1) * sizeof(int)))
    {
     #if IsDebug == 1
      printf("/n%d", *(ipMatrixComput + CurrentRow * iColComput));
     #endif

     for (iRowTmp = 0, iRow4Comput = 0; iRow4Comput < iColComput; iRow4Comput ++)
     {
      if (iRow4Comput != CurrentRow)
      {
       for (iCol4Comput = 0; iCol4Comput < iColComput -1; iCol4Comput ++)
       {
        *(ipNewMatrix + iRowTmp*(iColComput - 1) + iCol4Comput) = *(ipMatrixComput + iRow4Comput * iColComput + iCol4Comput + 1);
       }
       iRowTmp ++;
      }
     }

     #if IsDebug == 1
      for (iRow4Comput = 0; iRow4Comput < iColComput - 1; iRow4Comput ++)
      {
       for (iCol4Comput = 0; iCol4Comput < iColComput - 1; iCol4Comput ++)
       {
        printf(" %d ", *(ipNewMatrix + iRow4Comput * (iColComput - 1) + iCol4Comput));
       }
      }
      printf("]/n************/n");
      WriteMatrix(ipNewMatrix, iColComput - 1);
     #endif

     iResult += *(ipMatrixComput + CurrentRow * iColComput) * ComputMatrix(ipNewMatrix, iColComput - 1) * ((CurrentRow + 1) % 2 ? 1 : -1);

     #if IsDebug == 1
      printf("result=%d/n", iResult);
      getch();
     #endif

     free(ipNewMatrix);
    }
   }
  }
 }
 return (iResult);
}


void main(void)
{
 int *ipMatrix, iRowOrCol;
 char chSure;
 while (1)
 {
  printf("/n/nEnter the row or col of matrix:");
  scanf("%d", &iRowOrCol);
  if (iRowOrCol > 0)
  {
   if (ipMatrix = (int *)malloc(iRowOrCol * iRowOrCol * sizeof(int)))
   {
    while (1)
    {
     ReadMatrix(ipMatrix, iRowOrCol);
     WriteMatrix(ipMatrix, iRowOrCol);
     printf("Are you sure(Y/N)?");
     chSure = getche();

     if (toupper(chSure) == 'Y')
     {
      break;
     }
     else
     {
      printf("/n/n");
     }
    }
    printf("/n/n= %d/n", ComputMatrix(ipMatrix, iRowOrCol));
    free(ipMatrix);
   }
   else
   {
    printf("Error:No enough memory!/n");
    break;
   }
  }
  else
  {
   printf("Invalid col value!/n/n");
  }
  printf("/nCompute another Matrix(Y/N)?");
  chSure = getche();
  if (toupper(chSure) != 'Y')
  {
   break;
  }
 }
}
 

n阶行列式计算----c语言实现(完结)

花了半天时间,写了这个n阶行列式计算的程序,应该算是比较优美吧,有很多地方多次做了优化,程序占用内存不是很大,要是说小吧,也不合适,因为里边有一个递归,而且递归的深度还比较深。时间复杂度具体没有细看,...

c语言实现求行列式值

double hanglieshi(double *h,int i,int j) { int k,l; if(i!=j) { printf("Error"); return 0; } ...

递归求解N阶行列式

#include "stdafx.h"#include #include #include #include using namespace std; int sum; int sum1,su...

递归法求解矩阵行列式

方阵行列式的求解是线性代数学习的重要内容,行列式能显示矩阵的某些重要特性,如:行列式的值不为0,表示方阵可逆。此外,行列式在几何上表示空间几何体的体积,这里就不详细展开论述,具体可以参见相关线性代数教...

n阶行列式计算Python和C语言实现

行列式在数学中,是一个函数,其定义域为det的矩阵A,取值为一个标量,写作det(A)或 | A | 。无论是在线性代数、多项式理论,还是在微积分学中(比如说换元积分法中),行列式作为基本的数学工具,...
  • xyisv
  • xyisv
  • 2017年09月11日 23:43
  • 760

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

C++递归方法计算n阶行列式

#include #include using namespace std; double Fun(int n,double *ab) { double *aa = new double[n*...

用递归方法实现行列式求值

// Array.cpp : 定义控制台应用程序的入口点。 // #include "stdafx.h" #include #include using namespace std; int** pA...
  • the_k1
  • the_k1
  • 2016年11月11日 22:25
  • 381

Hanoi塔问题 栈与递归C语言编程实现

Hanoi塔 栈与递归C编程实现 参考书 严蔚敏  数据结构 #include #include #include typedef int ElemType; typedef int ...

行列式的计算-c编写-算法递归

  • 2011年01月05日 12:25
  • 1KB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:用C语言编程递归求行列式的值
举报原因:
原因补充:

(最多只允许输入30个字)