浅谈高斯消元

高斯消元

解一个n元的一次方程组，显然解这样的方程组我们需要得到n个方程。

大致思想

就像初中的时候解二元方程一样，先将所有方程某同一元前面的系数化为相同系数，抵消后只剩下$n-1$个方程和$n-1$个系数，当我们发现消到只剩一个元的时候显然可以直接求出，然后回代到两元的方程中，继续向多元回代直到解出$n$个元。
操作过程十分简单，但是不难发现存在非常严重的精度问题。

基本过程

把方程组的系数和值转换为一个$n*(n+1)$的増广矩阵（最后一列是值）。
当我们消元之后可以发现，矩阵的第$i$行前面有$i-1$个0，所以増广矩阵就变为一个阶梯矩阵。
从下往上回代。
如果在转换过程中无法找到系数不为0的元，就是无解的情况。
下面给出核心代码：

//消元过程
    for (int t=1;t<=n;t++){
        int h;
        for (int i=t;i<=n;i++) if (fabs(a[i][t])>1e-10){h=i;break;}
        swap(a[t],a[h]);
        for (int i=t+1;i<=n;i++){
            double x=a[i][t]/a[t][t];
            for (int j=t;j<=n+1;j++) a[i][j]-=a[t][j]*x;
        }
    }
//回代过程
    for (int i=n;i>=1;i--){
        ans[i]=a[i][n+1];
        for (int j=i+1;j<=n;j++) ans[i]-=ans[j]*a[i][j];
        ans[i]/=a[i][i];
    }