行列式模板

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

行列式~~~~

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);