#include<stdio.h>
#include<math.h>
double Z[4][4]; /*定义一个全局二维数组用来存放N-1阶余子式,因为A的伴随矩阵除以|A|时会产生小数，因此定义成float而非int*/
double R[4][4];
double daishu(double *p,int n);
void bansui(double Z[4][4],double R[4][4]);
/////////////////////////////////////////////
void main()
{
double w[4][4];
int m=4,i,j;
double r;
memset(w,0,sizeof(w));
memset(R,0,sizeof(R));
Z[0][0]=1;Z[0][1]=1;Z[0][2]=1;Z[0][3]=1;
Z[1][0]=0;Z[1][1]=1;Z[1][2]=1;Z[1][3]=1;
Z[2][0]=0;Z[2][1]=0;Z[2][2]=1;Z[2][3]=1;
Z[3][0]=0;Z[3][1]=0;Z[3][2]=0;Z[3][3]=1;
r=daishu(&Z[0][0],m);
// printf("%f\n",r);
bansui(Z,R); /*调用n_1()函数*/
for(i=0;i<m;i++) /*打印A的逆矩阵*/
{
for(j=0;j<m;j++)
{
w[i][j]=R[i][j]/r;
// printf("%f\t",w[i][j]);
}
// printf("\n");
}
}
double daishu(double *p,int n) /*计算行列式的函数*/
{
int k=0,i,s2=0,s1=0,j,s,t;
for(j=0;j<n;j++)
{
k=j;t=1;
for(i=1;i<=n;i++) /*对角线元素乘积*/
{
t=t*p[k];
if ((k+1)%n==0) k=k+1;
else k=k+n+1;
}
s1=s1+t; /*对角线乘积求和*/
}
for(j=0;j<n;j++)
{
k=j;t=1;
for(i=1;i<=n;i++)
{
t=t*p[k];
if(k%n==0) k=k+(2*n-1);
else k=k+(n-1);
}
s2=s2+t;
}
s=s1-s2;
return s;
}
/////////////////////////////////////////////////////
void bansui(double Z[4][4],double R[4][4]) /*把除第i行j列后的N-1阶矩阵的每个元素赋给一维数组d[]的函数*/
{
int i,j,e,f,l,m,k,t;
double q[3][3];
double y[4][4];
double tt;
//double daishu(int *p,int n);
memset(q,0,sizeof(q));
memset(y,0,sizeof(y));
for(i=0;i<4;i++)
{
l=i;
for(j=0;j<4;j++)
{
m=j;
k = 0;
for(e=0;e<4;e++)
{
t = 0;
if(e==l)
continue;
for(f=0;f<4;f++)
{
if(f==m)
continue;
q[k][t] = Z[e][f];
t++;
}
k++;
}
/*每得到一个一维数组d[]的值，便调用JS()函数,得到除第i行j列后的N-1阶矩阵的行列式的值,也即是余子式M[j]*/
tt=daishu(&q[0][0],3);
R[m][l]=tt*pow((-1),(l+m)); /*并把余子式的值存放在全局二维数组z[]中*/
}
}
}
