P3389 【模板】高斯消元法

发布时间：2020-12-14 05:15:37 所属栏目：大数据来源：网络整理

导读：P3389?【模板】高斯消元法注意每次选择当前列绝对值最大元素去消除其他的行是误差尽可能的小，本题特有的性质是，bi全都大于0，所以如果出现对角线上为0的情况是不合理的? 0x = bi #includebits/stdc++.husing namespace std;#define eps 1e-8#define maxn

P3389?【模板】高斯消元法

注意每次选择当前列绝对值最大元素去消除其他的行是误差尽可能的小，

本题特有的性质是，bi全都大于0，所以如果出现对角线上为0的情况是不合理的? 0x = bi

#include<bits/stdc++.h>
using namespace std;
#define eps 1e-8
#define maxn 123
double A[maxn][maxn];
int n,tp;
int main()
{
    scanf("%d",&n);
    for(int i=0; i<n; i++)
        for(int j=0; j<=n; j++)
            scanf("%lf",&A[i][j]);
    for(int i=0; i<n; i++)
    {
        tp=i;
        for(int j=i; j<n; j++)
            if(fabs(A[j][i])-fabs(A[tp][i])>eps)tp=j;
        if(tp!=i)
            for(int j=0; j<=n; j++)
                swap(A[i][j],A[tp][j]);
        if(fabs(A[i][i])<=eps)
        {
            printf("No Solutionn");
            return 0;
        }
        for(int j=i+1; j<=n; j++)
            A[i][j]/=A[i][i];
        for(int j=0; j<n; j++)
            if(i!=j)
                for(int k=i+1; k<=n; k++)
                    A[j][k]-=A[j][i]*A[i][k];
    }
    for(int i=0; i<n; i++)printf("%.2lfn",A[i][n]);
    return 0;
}

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!