/************************************************************************
功能: 用递归方法求行列式的值
作者: 黄志斌 于 广西河池
在自学《高等数学(二)》时，做行列式习题所得结果有些与教材所附
答案不一致，为了验证究竟答案的正确性，编写了这个程序。
其实用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;
  }
 }
}

发表于 @ 2006年07月10日 17:41:00|评论(loading...)

新一篇: [07-19] 解决灰鸽子新变种、Rootkit.Vanti.gen等及www.58111.com劫持(第5版) | 旧一篇: 一生悬命的程序员