#include<bits/stdc++.h>
using namespace std;
const int N=110;
double a[N][N];
int n;
void gauss()
{
for(int i=1;i<=n;i++){
int id=i;
for(int j=i+1;j<=n;j++)
if(fabs(a[j][i])>fabs(a[id][i]))
id=j;
if(fabs(a[id][i])<1e-6){
puts("No Solution");
return ;
}
for(int j=1;j<=n+1;j++)
swap(a[i][j],a[id][j]);
for(int j=1;j<=n;j++){
if(j!=i)
{
double temp=a[j][i]/a[i][i];
for(int k=i+1;k<=n+1;k++)
a[j][k]-=temp*a[i][k];
}
}
}
for(int i=1;i<=n;i++)
{
printf("%.2lf\n",a[i][n+1]/a[i][i]);
}
}
int main(){
cin>>n;
for(int i=1;i<=n;i++)
for(int j=1;j<=n+1;j++)
cin>>a[i][j];
gauss();
}
高斯消元模板
最新推荐文章于 2022-04-30 17:17:42 发布