【数学】高斯消元

Laoj1387 贾老二学数学

#include<iostream>
#include<cstdio>
using namespace std;
int n;
double a[104][104],ans[104];
void gauss(){
    for(int i=1;i<=n;i++){
        if(!a[i][i]){
            for(int j=i+1;j<=n;j++)
                if(a[j][i]) {swap(a[i],a[j]);break;}
        }
        for(int j=i+1;j<=n;j++){
            double x=a[j][i]/a[i][i];
            for(int k=i;k<=n+1;k++)a[j][k]-=a[i][k]*x;
        }
    }
    for(int i=n;i;i--){
        ans[i]=a[i][n+1]/a[i][i];
        for(int j=1;j<=i-1;j++) a[j][n+1]-=a[j][i]*ans[i];
    }
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n+1;j++)
            scanf("%lf",&a[i][j]);
    gauss();        
    for(int i=1;i<=n;i++) printf("%d ",int(ans[i]+0.5));
}