行列式模板

最新推荐文章于 2022-09-04 09:00:20 发布

natsuyu

最新推荐文章于 2022-09-04 09:00:20 发布

阅读量1.2k

点赞数

分类专栏：数论。

本文链接：https://blog.csdn.net/natsuyu/article/details/39523565

版权

数论。专栏收录该内容

19 篇文章 0 订阅

订阅专栏

那啥，计算行列式的。不过对于余子式就要稍微处理一下。

行列式~~~~

for (i=0,j=0;i<N&&j<N;i++,j++)
 {
  if (a[i][j]==0)
  {
   for (m=i;a[m][j]==0;m++);
  /* if (m==N)
   {
    sn=0;
    printf("detA=%lf\n",sn);
    exit(0);
   }
   else*/
    for (n=j;n<N;n++)
    {
     c=a[i][n];
     a[i][n]=a[m][n];
     a[m][n]=c;
    }
    k*=(-1);
  }
  for (s=N-1;s>i;s--)
  {
   x=a[s][j];
   for (t=j;t<N;t++)
    a[s][t]-=a[i][t]*(x/a[i][j]);
  }
 }
 for (i=0;i<N;i++)
  f*=a[i][i];
 sn=k*f;
 printf ("%.lf\n",sn);