求逆矩阵（C）

Lispboy

1. #include<stdio.h>  
   
2. #define N 10  
   
3. int getA(int arcs[N][N],int n)//按第一行展开计算|A|  
   
4. {  
   
5.         if(n==1)  
   
6.         {  
   
7.                 return arcs[0][0];  
   
8.         }  
   
9.         int ans = 0;  
   
10.         int temp[N][N];  
    
11.         int i,j,k;  
    
12.         for(i=0;i<n;i++)  
    
13.         {  
    
14.                 for(j=0;j<n-1;j++)  
    
15.                 {  
    
16.                         for(k=0;k<n-1;k++)  
    
17.                         {  
    
18.                                 temp[j][k] = arcs[j+1][(k>=i)?k+1:k];  
    
19. 
    
20.                         }  
    
21.                 }  
    
22.                 int t = getA(temp,n-1);  
    
23.                 if(i%2==0)  
    
24.                 {  
    
25.                         ans += arcs[0]*t;  
    
26.                 }  
    
27.                 else  
    
28.                 {  
    
29.                         ans -=  arcs[0]*t;  
    
30.                 }  
    
31.         }  
    
32.         return ans;  
    
33. }  
    
34. void getAStart(int arcs[N][N],int n,int ans[N][N])//计算每一行每一列的每个元素所对应的余子式，组成A*  
    
35. {  
    
36.         if(n==1)  
    
37.         {  
    
38.                 ans[0][0] = 1;  
    
39.                 return;  
    
40.         }  
    
41.         int i,j,k,t;  
    
42.         int temp[N][N];  
    
43.         for(i=0;i<n;i++)  
    
44.         {  
    
45.                 for(j=0;j<n;j++)  
    
46.                 {  
    
47.                         for(k=0;k<n-1;k++)  
    
48.                         {  
    
49.                                 for(t=0;t<n-1;t++)  
    
50.                                 {  
    
51.                                         temp[k][t] = arcs[k>=i?k+1:k][t>=j?t+1:t];  
    
52.                                 }  
    
53.                         }  
    
54. 
    
55. 
    
56.                         ans[j]  =  getA(temp,n-1);  
    
57.                         if((i+j)%2 == 1)  
    
58.                         {  
    
59.                                 ans[j] = - ans[j];  
    
60.                         }  
    
61.                 }  
    
62.         }  
    
63. }  
    
64. 
    
65. int main()  
    
66. {  
    
67.         int arcs[N][N];  
    
68.         int astar[N][N];  
    
69.         int i,j;  
    
70.         int n;  
    
71.         while(scanf("%d",&n)!=EOF && n)  
    
72.         {  
    
73.                 for(i=0;i<n;i++)  
    
74.                 {  
    
75.                         for(j=0;j<n;j++)  
    
76.                         {  
    
77.                                 scanf("%d",&arcs[j]);  
    
78.                         }  
    
79.                 }  
    
80. 
    
81.                 int a = getA(arcs,n);  
    
82.                 if(a==0)  
    
83.                 {  
    
84.                         printf("can not transform!\n");  
    
85.                 }  
    
86.                 else  
    
87.                 {  
    
88.                         getAStart(arcs,n,astar);  
    
89.                         for(i=0;i<n;i++)  
    
90.                         {  
    
91.                                 for(j=0;j<n;j++)  
    
92.                                 {  
    
93.                                         printf("%.3lf ",(double)astar[j]/a);  
    
94.                                 }  
    
95.                                 printf("\n");  
    
96.                         }  
    
97.                 }  
    
98.                 printf("\n");  
    
99. 
    
100.         }  
     
101.         return 0;  
     
102. }
     

#include<stdio.h>  <br /> #define N 10  <br /> int getA(int arcs[N][N],int n)//按第一行展开计算|A|  <br /> {  <br />         if(n==1)  <br />         {  <br />                 return arcs[0][0];  <br />         }  <br />         int ans = 0;  <br />         int temp[N][N];  <br />         int i,j,k;  <br />         for(i=0;i<n;i++)  <br />         {  <br />                 for(j=0;j<n-1;j++)  <br />                 {  <br />                         for(k=0;k<n-1;k++)  <br />                         {  <br />                                 temp[j][k] = arcs[j+1][(k>=i)?k+1:k];  <br /> <br />                         }  <br />                 }  <br />                 int t = getA(temp,n-1);  <br />                 if(i%2==0)  <br />                 {  <br />                         ans += arcs[0]<i>*t;  <br />                 }  <br />                 else  <br />                 {  <br />                         ans -=  arcs[0]<i>*t;  <br />                 }  <br />         }  <br />         return ans;  <br /> }  <br /> void getAStart(int arcs[N][N],int n,int ans[N][N])//计算每一行每一列的每个元素所对应的余子式，组成A*  <br /> {  <br />         if(n==1)  <br />         {  <br />                 ans[0][0] = 1;  <br />                 return;  <br />         }  <br />         int i,j,k,t;  <br />         int temp[N][N];  <br />         for(i=0;i<n;i++)  <br />         {  <br />                 for(j=0;j<n;j++)  <br />                 {  <br />                         for(k=0;k<n-1;k++)  <br />                         {  <br />                                 for(t=0;t<n-1;t++)  <br />                                 {  <br />                                         temp[k][t] = arcs[k>=i?k+1:k][t>=j?t+1:t];  <br />                                 }  <br />                         }  <br /> <br /> <br />                         ans[j]<i>  =  getA(temp,n-1);  <br />                         if((i+j)%2 == 1)  <br />                         {  <br />                                 ans[j]<i> = - ans[j]<i>;  <br />                         }  <br />                 }  <br />         }  <br /> }  <br /> <br /> int main()  <br /> {  <br />         int arcs[N][N];  <br />         int astar[N][N];  <br />         int i,j;  <br />         int n;  <br />         while(scanf("%d",&n)!=EOF && n)  <br />         {  <br />                 for(i=0;i<n;i++)  <br />                 {  <br />                         for(j=0;j<n;j++)  <br />                         {  <br />                                 scanf("%d",&arcs<i>[j]);  <br />                         }  <br />                 }  <br /> <br />                 int a = getA(arcs,n);  <br />                 if(a==0)  <br />                 {  <br />                         printf("can not transform!\n");  <br />                 }  <br />                 else  <br />                 {  <br />                         getAStart(arcs,n,astar);  <br />                         for(i=0;i<n;i++)  <br />                         {  <br />                                 for(j=0;j<n;j++)  <br />                                 {  <br />                                         printf("%.3lf ",(double)astar<i>[j]/a);  <br />                                 }  <br />                                 printf("\n");  <br />                         }  <br />                 }  <br />                 printf("\n");  <br /> <br />         }  <br />         return 0;  <br /> } <br />