高斯消元

最新推荐文章于 2022-09-07 17:45:26 发布
最新推荐文章于 2022-09-07 17:45:26 发布
阅读量1.1k 收藏 2
点赞数 1
分类专栏： 算法 文章标签： 数学 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_42638946/article/details/115736209
版权

高斯消元

1. 高斯消元原理

原理

  • 高斯消元是用于求解线性方程组的方法。给定如下方程组：

{ a 11 × x 1 + a 12 × x 2 + . . . . . . + a 1 n × x n = b 1 a 21 × x 1 + a 22 × x 2 + . . . . . . + a 2 n × x n = b 2 . . . . . . a n 1 × x 1 + a n 2 × x 2 + . . . . . . + a n n × x n = b n \begin{cases} a_{11} \times x_1 + a_{12} \times x_2 + ...... + a_{1n} \times x_n = b_1 \\ a_{21} \times x_1 + a_{22} \times x_2 + ...... + a_{2n} \times x_n = b_2 \\ ...... \\ a_{n1} \times x_1 + a_{n2} \times x_2 + ...... + a_{nn} \times x_n = b_n \\ \end{cases} ⎩⎪⎪⎪⎨⎪⎪⎪⎧​a11​×x1​+a12​×x2​+......+a1n​×xn​=b1​a21​×x1​+a22​×x2​+......+a2n​×xn​=b2​......an1​×x1​+an2​×x2​+......+ann​×xn​=bn​​

  • 我们是通过初等变换对上述方程组进行求解的，初等变换存在三种：

    （1）把某一行乘以一个非零的数；

    （2）交换某两行；

    （3）把某行的若干倍加到另一行上去；

  • 总的来说高斯消元分为两大步骤：第一步将系数矩阵通过初等变化变为上三角矩阵；回代求解。

  • 第一步：转为上三角矩阵：枚举每一列c

    （1）找到绝对值最大的一行；

    （2）将该行换到最上面；

    （3）将该行第1个数变为1；

    （4）将下面所有行的第c列消成0；

    类似于下图，但是缺少（2）（3）步。

在这里插入图片描述

  • 第二步：回代求解。

    从最后一个方程进行回代可以求出所有的答案。

  • 结果存在三种情况（做完第一步之后）：

    （1）唯一解：满秩，即系数矩阵没有全0行；

    （2）无解：系数矩阵存在全0行，且该行对应的b值不为0；

    （3）无穷多组解：系数矩阵存在全0行，且0行对应的b值都为0；

2. AcWing上的高斯消元题目

AcWing 883. 高斯消元解线性方程组

问题描述

分析

  • 根据上述分析实现代码即可。

代码

  • C++
#include <iostream>
#include <cmath>

using namespace std;

const int N = 110;
const double eps = 1e-6;  // 浮点数小于eps认为是0

int n;
double a[N][N];  // 存储增广矩阵

void print() {
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n + 1; j++)
            printf("%10.2lf", a[i][j]);
        puts("");
    }
    puts("");
}

// 返回0：存在唯一解; 返回1：存在无穷多组解; 返回2：无解
int gauss() {
    
    int r, c;  // 行、列
    for (r = 0, c = 0; c < n; c++) {
        // 第一步：(1) 找到绝对值最大的一行
        int t = r;  // 当前考察的为第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;  // 说明后面的行都为0
        
        // 第一步：(2) 将该行换到最上面
        for (int i = c; i < n + 1; i++) swap(a[t][i], a[r][i]);
        // 第一步：(3) 将第r行第c列对应的数变为1
        for (int i = n; i >= c; i--) a[r][i] /= a[r][c];
        // 第一步：(4) 将下面所有行的第c列消成0
        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) {  // 说明系数矩阵存在0行
        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++)  // j+1行到最后一行
            a[i][n] -= a[j][n] * a[i][j];
    
    return 0;
}

int main() {
    
    scanf("%d", &n);
    for (int i = 0; i < n; i++)
        for (int j = 0; j < n + 1; j++)
            scanf("%lf", &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;
}

AcWing 884. 高斯消元解异或线性方程组

问题描述

分析

  • 异或又可以被看成不进位的加法。我们将常规高斯消元中的加减法换成异或运算，乘法换成与运算即可。

  • 这里的步骤是和高斯消元法一样的，只不过需要将运算换成异或运算。

  • 第一步：转为上三角矩阵：枚举每一列c

    （1）找到非零行；

    （2）将该行换到最上面；

    （3）将下面所有行的第c列消成0；

  • 第二步：回代求解。

    从最后一个方程进行回代可以求出所有的答案。

  • 结果存在三种情况：

    （1）唯一解：做完第一步之后，满秩，即系数矩阵没有全0行；

    （2）无解：系数矩阵存在全0行，且该行对应的b值不为0；

    （3）无穷多组解：系数矩阵存在全0行，且0行对应的b值都为0；

代码

  • C++
#include <iostream>

using namespace std;

const int N = 110;

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

// 返回0：存在唯一解; 返回1：存在无穷多组解; 返回2：无解
int gauss() {
    
    int r, c;
    for (r = c = 0; c < n; c++) {
        int t = r;
        for (int i = r; i < n; i++)
            if (a[i][c]) {
                t = i;
                break;
            }
        
        if (!a[t][c]) continue;
        
        for (int i = c; i < n + 1; i++) swap(a[t][i], a[r][i]);
        
        for (int i = r + 1; i < n; i++)
            if (a[i][c]) 
                for (int j = c; j < n + 1; j++)
                    a[i][j] ^= a[r][j];
        
        r++;
    }
    
    if (r < n) {
        for (int i = r; i < n; i++)
            if (a[i][n])
                return 2;
        return 1;
    }
    
    for (int i = n - 1; i >= 0; i--)
        for (int j = i + 1; j < n; j++)
            a[i][n] ^= a[i][j] & a[j][n];
    
    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("%d\n", a[i][n]);
    } else if (t == 1) {
        puts("Multiple sets of solutions");
    } else puts("No solution");
    
    return 0;
}

AcWing 207. 球形空间产生器

问题描述

分析

  • 假设球心的坐标为： ( x 1 , x 2 , . . . , x n ) (x_1, x_2, ..., x_n) (x1​,x2​,...,xn​)，假设球的半径为R；给定的n+1个点的坐标为： ( a i , 1 , a i , 2 , . . . , a i , n ) (a_{i,1}, a_{i,2}, ..., a_{i,n}) (ai,1​,ai,2​,...,ai,n​)。根据圆上的点到圆心的距离相等，可得到方程组：

{ ( a 0 , 1 − x 1 ) 2 + ( a 0 , 2 − x 2 ) 2 + . . . + ( a 0 , n − x n ) 2 = R 2 ( 1 ) ( a 1 , 1 − x 1 ) 2 + ( a 1 , 2 − x 2 ) 2 + . . . + ( a 1 , n − x n ) 2 = R 2 ( 2 ) . . . . . . ( a n , 1 − x 1 ) 2 + ( a n , 2 − x 2 ) 2 + . . . + ( a n , n − x n ) 2 = R 2 ( n + 1 ) \begin{cases} (a_{0,1} - x_1)^2 + (a_{0,2} - x_2)^2 + ... + (a_{0,n} - x_n)^2 = R^2 \quad (1) \\ (a_{1,1} - x_1)^2 + (a_{1,2} - x_2)^2 + ... + (a_{1,n} - x_n)^2 = R^2 \quad (2) \\ ...... \\ (a_{n,1} - x_1)^2 + (a_{n,2} - x_2)^2 + ... + (a_{n,n} - x_n)^2 = R^2 \quad (n+1) \end{cases} ⎩⎪⎪⎪⎨⎪⎪⎪⎧​(a0,1​−x1​)2+(a0,2​−x2​)2+...+(a0,n​−xn​)2=R2(1)(a1,1​−x1​)2+(a1,2​−x2​)2+...+(a1,n​−xn​)2=R2(2)......(an,1​−x1​)2+(an,2​−x2​)2+...+(an,n​−xn​)2=R2(n+1)​

  • 理论上，上述方程组中有n+1个未知数，有n+1个方程，我们可以求出这些未知数，下面考虑如何求出这些未知数。

  • 考虑让上述方程组从第(2)项开始都减去第(1)项，现在考虑(2)-(1)中的第一项相减的结果：

( a 1 , 1 − x 1 ) 2 − ( a 0 , 1 − x 1 ) 2 = ( a 1 , 1 2 − a 0 , 1 2 ) − 2 ( a 1 , 1 − a 0 , 1 ) × x 1 (a_{1,1} - x_1)^2 - (a_{0,1} - x_1)^2 = (a_{1,1}^2 - a_{0,1}^2) - 2(a_{1,1} - a_{0,1}) \times x_1 (a1,1​−x1​)2−(a0,1​−x1​)2=(a1,12​−a0,12​)−2(a1,1​−a0,1​)×x1​

所以(2)-(1)的结果如下：
2 ( a 1 , 1 − a 0 , 1 ) × x 1 + . . . + 2 ( a 1 , n − a 0 , n ) × x n = ( a 1 , 1 2 − a 0 , 1 2 ) + . . + ( a 1 , n 2 − a 0 , n 2 ) 2(a_{1,1} - a_{0,1}) \times x_1 + ... + 2(a_{1,n} - a_{0,n}) \times x_n = (a_{1,1}^2 - a_{0,1}^2) + .. + (a_{1,n}^2 - a_{0,n}^2) 2(a1,1​−a0,1​)×x1​+...+2(a1,n​−a0,n​)×xn​=(a1,12​−a0,12​)+..+(a1,n2​−a0,n2​)
因此(i)-(1)的结果为：
2 ( a i , 1 − a 0 , 1 ) × x 1 + . . . + 2 ( a i , n − a 0 , n ) × x n = ( a i , 1 2 − a 0 , 1 2 ) + . . + ( a i , n 2 − a 0 , n 2 ) 2(a_{i,1} - a_{0,1}) \times x_1 + ... + 2(a_{i,n} - a_{0,n}) \times x_n = (a_{i,1}^2 - a_{0,1}^2) + .. + (a_{i,n}^2 - a_{0,n}^2) 2(ai,1​−a0,1​)×x1​+...+2(ai,n​−a0,n​)×xn​=(ai,12​−a0,12​)+..+(ai,n2​−a0,n2​)

  • 因为上述i可以取值范围是1~n，因此我们得到n个线性方程，一共n个未知数，可以使用高斯消元法求解。这里使用数组b替换上述系数，则有：

b i , 1 × x 1 + . . . + b i , n × x n = b i , n + 1 b_{i,1} \times x_1 + ... + b_{i,n} \times x_n = b_{i, n+1} bi,1​×x1​+...+bi,n​×xn​=bi,n+1​

代码

  • C++
#include <iostream>
#include <cmath>

using namespace std;

const int N = 15;

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

void gauss() {
    
    // 转化成上三角矩阵
    for (int r = 1, c = 1; c <= n; c++, r++) {
        // 找主元
        int t = r;
        for (int i = r + 1; i <= n; i++)
            if (fabs(b[i][c]) > fabs(b[t][c]))
                t = i;

        // 交换第t行和第r行
        for (int i = c; i <= n + 1; i++) swap(b[t][i], b[r][i]);
        // 归一化第r行
        for (int i = n + 1; i >= c; i--) b[r][i] /= b[r][c];
        // 消
        for (int i = c + 1; i <= n; i++)
            for (int j = n + 1; j >= c; j--)
                b[i][j] -= b[i][c] * b[r][j];
    }
    
    // 转换成对角矩阵(从系数矩阵最后一列开始消, 使用b[i][i]消上面的数据)
    for (int i = n; i > 1; i--)  // 枚举列
        for (int j = i - 1; j; j--) {  // 当前要把b[j][i]这个元素消成0
            b[j][n + 1] -= b[i][n + 1] * b[j][i];
            b[j][i] = 0;
        }
}

int main() {
    
    scanf("%d", &n);
    for (int i = 0; i < n + 1; i++)
        for (int j = 1; j <= n; j++)
            scanf("%lf", &a[i][j]);
    
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= n; j++) {
            b[i][j] = 2 * (a[i][j] - a[0][j]);
            b[i][n + 1] += a[i][j] * a[i][j] - a[0][j] * a[0][j];
        }
    
    gauss();
    
    for (int i = 1; i <= n; i++) printf("%.3lf ", b[i][n + 1]);
    
    return 0;
}

AcWing 208. 开关问题

问题描述

分析

  • 这里以一个例子讲解这个题目，有三个开关，分别为1、2、3，规则是：如果按1，则1、2的状态会改变；如果按2，则1、2、3的状态会改变；如果按3，则3的状态会改变。初始状态是0、0、0，转变为1、1、1存在多少种方案？

  • 我们使用 x 1 、 x 2 、 x 3 x_1、x_2、x_3 x1​、x2​、x3​ 分别表示每个开关有没有发生变化，没变为0，变了为1。根据上述规则，按1、2会影响1的状态，因此x1 ^ x2 = 1，同理我们可以得到如下异或方程组：

{ x 1 ⊕ x 2 = 1 第 一 个 开 关 x 1 ⊕ x 2 = 1 第 二 个 开 关 x 2 ⊕ x 3 = 1 第 三 个 开 关 \begin{cases} x_1 \oplus x_2 = 1 \quad 第一个开关 \\ x_1 \oplus x_2 = 1 \quad 第二个开关 \\ x_2 \oplus x_3 = 1 \quad 第三个开关 \end{cases} ⎩⎪⎨⎪⎧​x1​⊕x2​=1第一个开关x1​⊕x2​=1第二个开关x2​⊕x3​=1第三个开关​

  • 因此此时我们就将问题转化成为了AcWing 884. 高斯消元解异或线性方程组。对增广矩阵进行变换，变为上三角矩阵，最后如果有k个自由元（即k个行全0），最终的结果就为 2 k 2^k 2k。

  • 上述例题存在一个自由元，因此答案是2，两个答案分别是 ( x 1 = 0 , x 2 = 1 , x 3 = 0 ) 、 ( x 1 = 1 , x 2 = 0 , x 3 = 1 ) (x_1=0, x_2=1, x_3=0)、(x_1=1, x_2=0, x_3=1) (x1​=0,x2​=1,x3​=0)、(x1​=1,x2​=0,x3​=1)。

代码

  • C++
#include <iostream>
#include <cstring>

using namespace std;

const int N = 35;

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

int gauss() {
    
    int r, c;
    for (r = 1, c = 1; c <= n; c++) {
        // 找主元
        int t = r;
        for (int i = r; i <= n; i++)
            if (a[i][c])
                t = i;
        if (a[t][c] == 0) continue;
        // 交换
        for (int i = c; i <= n + 1; i++) swap(a[r][i], a[t][i]);
        // 消
        for (int i = r + 1; i <= n; i++)
            for (int j = n + 1; j >= c; j--)
                a[i][j] ^= a[i][c] & a[r][j];
        r++;
    }
    
    int res = 1;
    if (r < n + 1) {
        for (int i = r; i <= n; i++) {
            if (a[i][n + 1]) return -1;  // 出现了 0 == !0，无解
            res *= 2;
        }
    }
    return res;
}

int main() {
    
    int T;
    scanf("%d", &T);
    while (T--) {
        scanf("%d", &n);
        memset(a, 0, sizeof a);
        for (int i = 1; i <= n; i++) scanf("%d", &a[i][n + 1]);
        for (int i = 1; i <= n; i++) {
            int x;
            scanf("%d", &x);
            a[i][n + 1] ^= x;  // a[i][n + 1]存储的是第i个灯的状态是否变化
            a[i][i] = 1;  // 默认开关可以控制自己
        }
        
        int x, y;
        while (scanf("%d%d", &x, &y), x || y) a[y][x] = 1;  // 影响的是第y个方程
        
        int t = gauss();
        if (t == -1) puts("Oh,it's impossible~!!");
        else printf("%d\n", t);
    }
    
    return 0;
}
你好世界wxx
关注
  • 1
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
高斯相除
01-26 2190
The distribution of Complex Gaussian is (1)Nc(x|a,A)∝exp⁡[−(x−a)HA−1(x−a)]" role="presentation" style="position: relative;">Nc(x|a,A)∝exp[−(x−a)HA−1(x−a)](
利用高斯消元法，列主元法求解线性方程组
11-27
高斯顺序消元法程序简单，操作便捷，结果比较精确，但是每次运算的时候必须保证对角线上元素不为0，否则将无法计算，并且当主元过小的时候（病态情况）会产生较大误差，影响算法的稳定性。通过对比，列主元法通过判断大小并行行交换的操作能够有效地克服高斯顺序消元法的缺陷。
高斯消元（Gaussian elimination）
努力中的老周的专栏
11-24 2万+
我们将分两个部分来说明高斯消去，第一部分是数学方面，第二个部分是如何用 C++ 来实现。 第一部分 数学知识 什么是高斯消去 高斯消元法（或译：高斯消去法），是线性代数规划中的一个算法，可用来为线性方程组求解。 高斯消元的目标是将矩阵变为上三角矩阵，例如，对应的矩阵为就是一个上三角矩阵（upper triangular matrix）。 对于一个线性方程组，我们可以写出两个矩阵（系数矩阵 coefficient matrix 和增广矩阵 Augmented Matrix），如下图所示。 .
高斯消元法讲解
辉小歌的博客
05-29 1万+
高斯消元的本质就是化简成一个阶梯式的行列式。 首先线性方程组的解有以下三种情况： 无解 有无穷多个解 有唯一解 高斯消元的步骤分为以下四步： 枚举每一行找到当前行(包括当前行)下面的，当前列的绝对值最大的一个数。 将其绝对值最大的一行移到上面 将改行的第一个数变成1 将下面所有的行的当前列都清成0 例子： 注意：这里刚开始是第一行和第一列，以此类推就会化简成一个阶梯行列式。 https://www.acwing.com/problem/content/description/885/ 按照上面.
高斯消元法解线性方程组数学原理及Python手工实现
杨老师的小跟班(*∩_∩*)
09-07 2404
高斯消元法解线性方程组数学原理及Python手工实现 先逐次消去未知数，将方程化为同解的上三角形方程组，此过程称为消元。然后按方程相反顺序求解上三角形方程组，得到原方程组的解，此过程称为回代。这种方法称为高斯消元法。
数论之高斯消元解线性方程，组合数的四种求法
lihua777的博客
02-17 1517
高斯消元法 例题： 高斯消元步骤： ！！！枚举每一列col 找到绝对值最大的那一行 将该行换到最上面 将该行的第一个数变成1 将下面所有行的第col列的数消成0（通过原系数-原系数*上一个方程为什么要变成1的原因） 核心思想： 将多组方程构造为阶梯型（三角形）：即如下图所示，那么就可以由下层依次推出上层的解 多种情况：（可以看上面的人话版本） 无解:若在最后化成的上三角形矩阵中,正对角线中某个元素为0,但其所在行的最后一列元素不为0时,此时矩阵无解 有无数解:若在最后化..
高斯消元（非常详细） 高斯消元（非常详细）
最新发布
06-19
高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元...
gaosixiaoyuanfa.rar_高斯消元
09-23
这个名为“gaosixiaoyuanfa.rar_高斯消元”的压缩包文件包含了一份详细的高斯消元法课程资料，通常以PPT的形式呈现，可以帮助学习者深入理解并掌握这一算法。 高斯消元法的基本思想是通过一系列行初等变换，将系数...
gauss_seq.c.zip_高斯消元
09-23
此外，高斯-约旦消元（Gauss-Jordan elimination）是高斯消元的一种扩展，它不仅将矩阵转化为上三角形，还能进一步将其化为单位矩阵，从而直接得到解。 总结来说，`gauss_seq.c`是一个用C语言实现的高斯消元法串行...
全主元高斯消元
05-31
在VC6.0环境下，可以使用标准C++库来创建、读取和输出矩阵，编写函数实现全主元高斯消元算法，并在主函数中调用这些函数，最后打印出解。通过调试和测试，可以优化代码并提高其效率。 总之，全主元高斯消元法是解决...
高斯消元附解方程神器
04-09
先输入一个数n表示方程元数。 然后再输入一个 n*(n+1)的矩阵 具体见我的博客
高斯消元法详解
热门推荐
ityanger的技术栈
05-05 10万+
1.概述 高斯消元法(Gaussian elimination)是求解线性方阵组的一种算法，它也可用来求矩阵的秩，以及求可逆方阵的逆矩阵。它通过逐步消除未知数来将原始线性系统转化为另一个更简单的等价的系统。它的实质是通过初等行变化(Elementary row operations)，将线性方程组的增广矩阵转化为行阶梯矩阵(row echelon form). 2.原理 已知线性方程组: ...
【数据结构与算法】高斯消元法
zzy_NIC的博客
11-02 6792
高斯消元法 数学上，高斯消元法（Gaussian Elimination或译：高斯消去法），是线性代数规划中的一个算法，可用来为线性方程组求解。在许多应用中，我们需要解一个包含n个方程的n元联立方程组： 其中，n是一个大数。高斯消去法的思路是把n个线性方程构成的n元联立方程组变成一个等价方程组（也就是说，它的解和原来的方程组相同），该方程组有着一个上三角形的系数矩阵，这种矩阵的主对角线下方元素全部为0。 用矩阵的符号，我们可以把它写成Ax=b=>A’x=b′{\boldsymbol{A}x=\bo
高斯消元法
Dextrad_ihacker的博客
11-26 1019
高斯消元法，是线性代数中的一个算法，可用来求解线性方程组，并可以求出矩阵的秩，以及求出可逆方阵的逆矩阵。 高斯消元法的原理是： 若用初等行变换将增广矩阵 化为 ，则AX = B与CX = D是同解方程组。 所以我们可以用初等行变换把增广矩阵转换为行阶梯阵，然后回代求出方程的解。以上是线性代数课的回顾，下面来说说高斯消元法在编程中的应用。首先，先介绍程序中高斯消元法的步骤： (我们设方程组中方
Gauss消元法matlab实现（原理、代码、案例及详解）
m0_56603583的博客
05-01 2万+
第一部分：问题分析 （1）实验题目：高斯消元算法 具体实验要求：要求学生运用高斯列主元消元法计算出线性方程组Ax=b的近似解。用matlab编写高斯列主元消元法的代码，要求代码实现用户输入了矩阵行列数、稀疏矩阵A、行列式b之后，程序能够输出迭代的近似解。 实验目的：让同学们进一步掌握高斯列主元消元法的原理以及迭代过程，并且通过matlab编程培养实际的上机操作能力和代码能力。 第二部分：数学原理 列主元素消去法是为控制舍入误差而提出来的一种算法，列主元素消去法计算基本上能控制舍入误差的影响，其..
高斯消元快速入门
Pengwill's Blog
08-15 3万+
高斯消元快速入门 一、基本描述 学习一个算法/技能，首先要知道它是干什么的，那么高斯消元是干啥的呢？ 高斯消元主要用来求解线性方程组，也可以求解矩阵的秩，矩阵的逆。在ACM中是一个有力的数学武器. 它的时间复杂度是n^3，主要与方程组的个数，未知数的个数有关。 那么什么是线性方程组呢? 简而言之就是有多个未知数，并且每个未知数的次数均为一次，这样多个未知数组成的方程组为线性...
线段树
weixin_42638946的博客
04-08 1万+
线段树 1. 线段树原理 原理 线段树存在五个操作 （1）pushup：由子节点计算父节点的信息； （2）pushdown：把当前父节点的修改信息下传到子节点，也被称为懒标记（延迟标记）；这个操作比较复杂，一般不涉及到区间修改则不用写。 （3）build：将一段区间初始化成线段树； （4）modify：修改操作，分为两类：① 单点修改（需要使用pushup），② 区间修改（需要使用pushdown）； （5）query：查询一段区间的值。 线段树的原理：线段树除了最后一层之外，是一颗满二叉树，假
matlab 高斯消元
11-12
高斯消元是一种线性代数中常用的求解线性方程组的方法，Matlab中也提供了相应的函数来实现高斯消元。下面是一个简单的例子： 假设有如下线性方程组： ``` 2x + y - z = 8 -3x - y + 2z = -11 -2x + y + 2z = -3 ``...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值