【洛谷 P3389】高斯消元法(高斯消元模板)

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

摘要生成于 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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值