P3389 【模板】高斯消元法

2020-09-10 16:33:19 阅读：352 来源： 互联网

高斯消元模板题

把第i列除了第i行外所有的系数变成0

#include <bits/stdc++.h>
#define inf 2333333333333333
#define N 110
#define p(a) putchar(a)
#define For(i,a,b) for(int i=a;i<=b;++i)
//by war
//2020.9.10
using namespace std;
int n;
double t,eps=1e-7;
double a[N][N];
void in(int &x){
    int y=1;char c=getchar();x=0;
    while(c<'0'||c>'9'){if(c=='-')y=-1;c=getchar();}
    while(c<='9'&&c>='0'){ x=(x<<1)+(x<<3)+c-'0';c=getchar();}
    x*=y;
}
void o(int x){
    if(x<0){p('-');x=-x;}
    if(x>9)o(x/10);
    p(x%10+'0');
}

void gs(){
    For(i,1,n){
        t=a[i][i];
        For(j,i+1,n) t=fmax(t,a[j][i]);
        if(fabs(t)<eps){
            puts("No Solution");
            return ;
        }
        For(j,1,n){
            if(i!=j){
                double temp=a[j][i]/a[i][i];
                For(k,i,n+1){
                    a[j][k]-=a[i][k]*temp;
                }
            }
        }
    }
    For(i,1,n) a[i][n+1]/=a[i][i];
    For(i,1,n) printf("%.2lf\n",a[i][n+1]);
}

signed main(){
    in(n);
    For(i,1,n)
        For(j,1,n+1)
            scanf("%lf",&a[i][j]);
    gs();
    return 0;
}

标签：10,gs,int,P3389,高斯消,模板,元法,define
来源： https://www.cnblogs.com/war1111/p/13646540.html

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

ICode9

P3389 【模板】高斯消元法