矩阵逆的计算

第 2328 题  ：  矩阵逆的计算
(时间限制为：1000毫秒)

 

    请你设计一个计算器，计算出给你的一个矩阵的逆矩阵。

输入：

输入的第一行表示需要你计算的矩阵个数N(N<50)，每个矩阵的第一行为该矩阵的阶数n，接下来的n行，每行为由空格分隔的n个实数。

输出：

如果给定的矩阵有逆，就使用一行输出“Matrix 1 ：”，接下来使用n行输出计算的结果，每行输出n个数，每个数之间使用一个空格隔开,输出结果数据保留6位小数；如果没有逆，直接输出“Matrix 1 : NO.”；

 

样例输入：

2

3

1 2 3

2 2 1

3 4 3

3

2 3 -1

-1 3 5

1 6 4

 

样例输出：

Matrix 1 :
1.000000 3.000000 -2.000000
-1.500000 -3.000000 2.500000
1.000000 1.000000 -1.000000
Matrix 2 : NO.

#include<iostream>#include<cmath>#include<iomanip>using namespace std;int n,N,count1;double mat[51][102];//高斯消元法解方程bool jiefangchen(int n,int x_num){//n表示方程数目，x_num表示未知数变量数目，mat的第x_num+1列存储常数bi值。int i,j,m;for(i=0;i<x_num/2;i++){//i列for(j=i;j<n;j++)if(fabs(mat[j][i])>0.000001)break;//找到第i行第i列为非零if(j==n)return false;if(j!=i){for(int t=0;t<x_num;t++)//交换两行swap(mat[i][t],mat[j][t]); }double a=mat[i][i];//对角线系数化为1for(m=i;m<x_num;m++)mat[i][m]/=a;//消元法for(m=0;m<2*n;m++){if(m==i)continue;double b;b=mat[m][i];for(int j=i;j<x_num;j++)mat[m][j]-=b*mat[i][j];}}return true;}int main(){int i,j;count1=0;cin>>N;while(N--){count1++;cin>>n;for(i=0;i<n;i++){for(j=0;j<n;j++)cin>>mat[i][j];}for(i=0;i<n;i++){for(j=n;j<2*n;j++){if(j-n==i){mat[i][j]=1;}else{mat[i][j]=0;}}}if(jiefangchen(n,2*n)){cout<<"Matrix "<<count1<<" :"<<endl;for(j=0;j<n;j++){for(i=n;i<2*n;i++)printf("%f ",mat[j][i]);//cout<<setiosflags(ios::fixed)<<setprecision(6)<<mat[j][i]<<" ";cout<<endl;}}else{cout<<"Matrix "<<count1<<" : NO."<<endl;}}return 0;}