高斯消元板子题

Problem Link

#include <bits/stdc++.h>
using namespace std;
const int N = 15;
double b[N][N],a[N][N];
int n;
void guass(){
    for(int r=1,c=1;c<=n;c++,r++){
        int t=r;
        for(int i=r;i<=n;i++)
            if(fabs(a[i][c])>fabs(a[t][c]))t=i;
        for(int i=1;i<=n+1;i++)swap(a[t][i],a[r][i]);
        for(int i=n+1;i>=c;i--)a[r][i]/=a[r][c];
        for(int i=r+1;i<=n;i++)
            for(int j=n+1;j>=c;j--)a[i][j]-=a[r][j]*a[i][c];
    }
    for(int c=n;c>1;c--)
        for(int r=c-1;r>=1;r--){
            a[r][n+1]-=a[r][c]*a[c][n+1];
            a[r][c]=0;
        }
        
}
int main()
{
    cin>>n;
    for(int i=0;i<n+1;i++)
        for(int j=1;j<=n;j++)cin>>b[i][j];
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++){
            a[i][j]=2*(b[i][j]-b[0][j]);
            a[i][n+1]+=-b[0][j]*b[0][j]+b[i][j]*b[i][j];
        }
    guass();
    for(int i=1;i<=n;i++)printf("%.3lf ",a[i][n+1]);
}