MnZn刚学矩阵-114514秒，50pts求助

评测记录

#include<cstdio>
#define L 110
int n;
double a[L][L];
int main(){
    scanf("%d",&n);
    for(register int i=1;i<=n;i++)
        for(register int j=1;j<=n+1;j++)scanf("%lf",&a[i][j]);
    for(register int j=1;j<=n;j++){
        if(!a[j][j]){
            if(a[j][n+1])printf("-1");
            else printf("0");
            return 0;
        }
        for(register int i=1;i<=n;i++)
            if(i!=j){
                double tmp=a[i][j]/a[j][j];
                for(register int k=1;k<=n+1;k++)a[i][k]-=tmp*a[j][k];
            }
    }
    for(register int i=1;i<=n;i++)printf("x%d=%.2lf\n",i,a[i][n+1]/a[i][i]);
    return 0;
}