解题思路
留着自己看的高斯消元模板。
代码:
#include<algorithm>
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<string>
#include<cstdlib>
using namespace std;
const double eps=1e-8;
double f[110][110];
int n;
void gauss(){
for(int i=1;i<=n;i++){
int pre=i;
for(int j=i+1;j<=n;j++)
if(fabs(f[j][i]-f[pre][i])<=eps)
pre=j;
for(int j=1;j<=n+1;j++) swap(f[i][j],f[pre][j]);
if(fabs(f[i][i])<=eps){printf("No Solution");exit(0);}
for(int j=n+1;j>=i;j--) f[i][j]/=f[i][i];
for(int j=1;j<=n;j++)
if(i!=j){
for(int k=i+1;k<=n+1;k++) f[j][k]-=f[j][i]*f[i][k];
f[j][i]=0;
}
}
}
int main(){
scanf("%d",&n);
for(int i=1;i<=n;i++)
for(int j=1;j<=n+1;j++)
scanf("%lf",&f[i][j]);
gauss();
for(int i=1;i<=n;i++) printf("%.2lf\n",f[i][n+1]);
return 0;
}