【luogu3389】【模板】高斯消元法 [高斯消元]

P3389 【模板】高斯消元法

 

#include<bits/stdc++.h>
using namespace std;
#define rg register
const int N=100+5;
const double eps=1e-20;
int n;
double a[N][N],ans[N];

template <class t>void rd(t &x)
{
    x=0;int w=0;char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    x=w?-x:x;
}
//fabs求浮点数x的绝对值
int main()
{
    rd(n);
    for(rg int i=1;i<=n;++i)
    for(rg int j=1;j<=n+1;++j) scanf("%lf",&a[i][j]);
    for(rg int i=1;i<=n;++i)//枚举定量 
    {
        int now=i;
        for(rg int j=i+1;j<=n;++j)
        if(fabs(a[now][i])<fabs(a[j][i])) now=j;//找最大的第i主元
        if(fabs(a[now][i])<eps) {printf("No Solution");exit(0);}
        if(now!=i) swap(a[now],a[i]);
    //    for(int j=1;j<=n+1;++j) swap(a[now][j],a[i][j]);//交换 
        double div=a[i][i];
        for(rg int j=i;j<=n+1;++j) a[i][j]/=div;//将后面的序数都除div 
        for(rg int j=i+1;j<=n;++j)
        {
            div=a[j][i];
            for(rg int k=i;k<=n+1;++k) a[j][k]-=a[i][k]*div;
        }
    }
    ans[n]=a[n][n+1];
    for(rg int i=n-1;i>=1;--i)
    {
        ans[i]=a[i][n+1];
        for(rg int j=i+1;j<=n;++j) ans[i]-=a[i][j]*ans[j];
    }
    for(rg int i=1;i<=n;++i) printf("%.2f\n",ans[i]);
    return 0;
}