标签：洛谷 int 高斯消 include 模板 P3389 110

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

洛谷传送门

题目背景

Gauss消元

题目描述

给定一个线性方程组，对其求解

输入格式

第一行，一个正整数 nn

第二至 n+1n+1行，每行 n+1n+1 个整数，为a_1, a_2 \cdots a_na1,a2⋯a**n 和 bb，代表一组方程。

输出格式

共n行，每行一个数，第 ii行为 x_ix**i （保留2位小数）

如果不存在唯一解，在第一行输出"No Solution".

题解：

高斯消元模板题。

关于高斯消元，请走：

浅谈高斯消元

代码：

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const double eps=1e-6;
int n;
double a[110][110],b[110];
bool flag;
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
            scanf("%lf",&a[i][j]);
        scanf("%lf",&b[i]);
    }
    for(int i=1;i<=n;i++)
    {
        flag=0;
        for(int j=i;j<=n;j++)
        {
            if(flag)
                break;
            if(fabs(a[j][i])>eps)
            {
                for(int k=1;k<=n;k++)
                    swap(a[i][k],a[j][k]);
                swap(b[i],b[j]);
                flag=1;
            }
        }
        if(!flag)
        {
            puts("No Solution");
            return 0;
        }
        for(int j=1;j<=n;j++)
        {
            if(i==j)
                continue;
            double rate=a[j][i]/a[i][i];
            for(int k=i;k<=n;k++)
                a[j][k]-=a[i][k]*rate;
            b[j]-=b[i]*rate;
        }
    }
    for(int i=1;i<=n;i++)
        printf("%.2lf\n",1.0*b[i]/a[i][i]);
    return 0;
}

标签：洛谷,int,高斯消,include,模板,P3389,110
来源： https://www.cnblogs.com/fusiwei/p/13848151.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。