【洛谷 P3389】高斯消元法（高斯消元模板）

最新推荐文章于 2021-07-16 11:11:04 发布

原创最新推荐文章于 2021-07-16 11:11:04 发布 · 411 阅读

0 ·

CC 4.0 BY-SA版权

模板专栏收录该内容

7 篇文章

订阅专栏

本文提供了一个高斯消元法的实现模板，通过该模板可以解决线性方程组问题。代码中详细展示了如何使用高斯消元法进行求解，包括矩阵的初始化、主元选择、行列式的标准化等步骤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门~

解题思路

留着自己看的高斯消元模板。
代码：

#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;
}