数论-中国剩余定理、高斯消元

中国剩余定理（先挖坑）

高斯消元 $O（n^3）$

求解线性方程组

初等行列变换

1. 把某一行乘以一个非0的数（等式俩边同时乘以一个非零数）
2. 交换某2行（方程交换位置）
3. 把某行的若干倍加到另一行上去（把一个方程的若干倍加到另一个方程上去）

变换成上三角形式后

• 完美阶梯型----------唯一解
• 0=非0 ----------------无解
• 0=0--------------------无穷多组解

算法步骤

• 枚举每一列c

1. 找到当前列绝对值最大的一行
2. 将该行换到最上面
3. 将该行第一个数变成1
4. 将下面所有行的当前列消成0

模板

#include <iostream>
#include <algorithm>
#include <cmath>

using namespace std;

const int N = 110;
const double eps = 1e-6;

int n;
double a[N][N];


int gauss()
{
    int c, r;
    for (c = 0, r = 0; c < n; c ++ )
    {
        int t = r;
        for (int i = r; i < n; i ++ )
            if (fabs(a[i][c]) > fabs(a[t][c]))
                t = i;

        if (fabs(a[t][c]) < eps) continue;

        for (int i = c; i < n + 1; i ++ ) swap(a[t][i], a[r][i]);
        for (int i = n; i >= c; i -- ) a[r][i] /= a[r][c];

        for (int i = r + 1; i < n; i ++ )
            if (fabs(a[i][c]) > eps)
                for (int j = n; j >= c; j -- )
                    a[i][j] -= a[r][j] * a[i][c];

        r ++ ;
    }

    if (r < n)
    {
        for (int i = r; i < n; i ++ )
            if (fabs(a[i][n]) > eps)
                return 2;
        return 1;
    }

    for (int i = n - 1; i >= 0; i -- )
        for (int j = i + 1; j < n; j ++ )
            a[i][n] -= a[j][n] * a[i][j];

    return 0;
}

int main()
{
    cin >> n;
    for (int i = 0; i < n; i ++ )
        for (int j = 0; j < n + 1; j ++ )
            cin >> a[i][j];

    int t = gauss();

    if (t == 0)
    {
        for (int i = 0; i < n; i ++ ) printf("%.2lf\n", a[i][n]);
    }
    else if (t == 1) puts("Infinite group solutions");
    else puts("No solution");

    return 0;
}