高斯消元题集

最新推荐文章于 2021-11-19 19:24:00 发布
最新推荐文章于 2021-11-19 19:24:00 发布
阅读量748 收藏
点赞数
分类专栏： Gauss-Jordan Elimination 文章标签： Gauss-Jordan Elimina
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Mrx_Nh/article/details/75196333
版权

高斯消元题集

高斯消元本质上不过是线性代数中对矩阵进行的初等变换而已，结题的关键是如何构造出增广矩阵，再看方程的类型，比如整数解、浮点数解、模线性方程、或者求自由变元的数目等等再分别用不同的方式解决即可。


        关于高斯消元的原理，可以参考我的另一篇博客：高斯消元原理


HDU: 3976 4818

POJ: 1222 1830 2065 2947 3185


HDU 3976 

Electric resistance

#include <bits/stdc++.h>

using namespace std;
const int MAXN = 284;
const double eps = 1e-9;

double a[MAXN][MAXN];
double x[MAXN];
int free_x[MAXN];

int Gauss(int equ, int var) {
  int col, max_r;
  for(int k = 0, col = 0; k < equ && col < var; ++k, ++col) {
    max_r = k;
    for(int i = k + 1; i < equ; ++i) {
      if(fabs(a[i][col]) > fabs(a[max_r][col])) max_r = i;
    }
    if(fabs(a[max_r][col]) < eps) return 0; // No solution
    if(max_r != k) {
      for(int j = col; j < var; ++j) {
        swap(a[k][j], a[max_r][j]);
      }
      swap(x[k], x[max_r]);
    }
    x[k] /= a[k][col];
    for(int j = col + 1; j < var; ++j) a[k][j] /= a[k][col];
    a[k][col] = 1;
    for(int i = 0; i < equ; ++i) {
      if(i != k) {
        x[i] -= x[k] * a[i][col];
        for(int j = col + 1; j < var; ++j) a[i][j] -= a[k][j] * a[i][col];
        a[i][col] = 0;
      }
    }
  }
  return 1;
}

int main() {
  int T, cas = 1;
  scanf("%d", &T);
  while(T--) {
    int n, m;
    scanf("%d %d", &n, &m);
    int u, v, w;
    memset(a, 0, sizeof(a));
    while(m--) {
      scanf("%d %d %d", &u, &v, &w);
      a[u-1][v-1] += 1.0 / w;
      a[u-1][u-1] += -1.0 / w;
      a[v-1][u-1] += 1.0 / w;
      a[v-1][v-1] += -1.0 / w;  
    }
    for(int i = 0; i < n; ++i) x[i] = 0;
    x[0] = 1;
    for(int i = 0; i < n; ++i) {
      a[n-1][i] = 0;
    }
    a[n-1][0] = 1;
    Gauss(n, n);
    printf("Case #%d: %.2f\n", cas++, x[n-1]);
  }
  return 0;
}



HDU 4818

RP problem

#include <bits/stdc++.h>

using namespace std;
const int MAXN = 284;
const double eps = 1e-6;

double a[MAXN][MAXN], x[MAXN];
int equ, var;

int Gauss() {
  int col, max_r;
  for(int k = 0, col = 0; k < equ && col < var; ++k, ++col) {
    max_r = k;
    for(int i = k + 1; i < equ; ++i) {
      if(fabs(a[i][col]) > fabs(a[max_r][col])) {
        max_r = i;
      }
    }
    if(fabs(a[max_r][col]) < eps) return 0; // No solution
    if(max_r != k) {
      for(int j = col; j < var; ++j) {
        swap(a[k][j], a[max_r][j]);
      }
      swap(x[k], x[max_r]);
    }
    x[k] /= a[k][col];
    for(int j = col + 1; j < var; ++j) a[k][j] /= a[k][col];
    a[k][col] = 1;
    for(int i = 0; i < equ; ++i) {
      if(i != k) {
        x[i] -= x[k] * a[i][col];
        for(int j = col + 1; j < var; ++j) a[i][j] -= a[k][j] * a[i][col];
        a[i][col] = 0;
      }
    }
  }
  return 1;
}

vector<int> linker[MAXN];
int g[MAXN][MAXN], deg[MAXN], add[MAXN];
inline void init(int n) {
  for(int i = 0 ; i < n; ++i) {
    linker[i].clear();
    deg[i] = x[i] = 0;
  }
  memset(g, 0, sizeof(g));
  memset(a, 0, sizeof(a));
}

int main() {
  int T;
  scanf("%d", &T);
  while(T--) {
    int n, m;
    scanf("%d %d", &n, &m);
    equ = var = n;
    init(n);
    int u, v;
    while(m--) {
      scanf("%d %d", &u, &v);
      if(u == v) continue;
      g[u][v] = 1;
    }
    for(int i = 0; i < n; ++i) {
      for(int j = 0; j < n; ++j) {
        if(j != i && g[i][j]) {
          deg[i]++;
          linker[j].push_back(i);
        }
      } 
    }
    for(int i = 0; i < n; ++i) {
      a[i][i] = -1;
      for(int j = 0; j < linker[i].size(); ++j) {
        int v = linker[i][j];
        if(i == v) continue;
        a[i][v] = 1.0 / deg[v];
      }
    }
    for(int i = 0; i < n; ++i) {
      a[n-1][i] = 1;
    }
    x[n-1] = 1;

    for(int k = 0; k < n-1; ++k) {
      if(g[n-1][k] == 0) {
        for(int i = 0; i < n-1; ++i) {
          if(g[n-1][i]) a[i][var] = 1.0 / (deg[n-1] + 1);
          else a[i][var] = 0;
        }                                               
        a[k][var] = 1.0 / (deg[n-1] + 1);              
        a[n-1][var] = 1;
        add[var] = k;
        var++;
      }
    }
    if(!Gauss()) {
      puts("INF");
      continue;
    }
    double now = x[n-1];
    int ans  = -1;
    for(int i = n; i < var; ++i) {
      if(x[n-1] / a[n-1][i] > now) {
        ans = add[i];
        now = x[n-1] / a[n-1][i];
      }
    }
    printf("1 %d\n", ans);
  }
  return 0;
}


POJ 1222

EXTENDED LIGHTS OUT

  

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>

using namespace std;
const int MAXN = 128;

int a[MAXN][MAXN];
int x[MAXN];
bool free_x[MAXN];

int gcd(int a, int b) {
  return b == 0 ? a : gcd(b, a % b);
}

int lcm(int a, int b) {
  return a / gcd(a, b) * b;
}

int Gauss(int equ, int var) {
  for(int i = 0; i <= var; ++i) {
    x[i] = 0;
    free_x[i] = true;
  }
  int max_r;
  for(int col = 0, row = 0; col < var && row < equ; ++col, ++row) {
    max_r = row;
    for(int i = row + 1; i < equ; ++i) {
      if(abs(a[max_r][col]) < abs(a[i][col])) {
        max_r = i;
      }
    }
    if(max_r != row) {
      for(int j = row; j < var + 1; ++j) {
        swap(a[row][j], a[max_r][j]);
      }
    }
    if(a[row][col] == 0) {
      row--;
      continue;
    }
    for(int i = row + 1; i < equ; ++i) {
      if(a[i][col] != 0) {
        for(int j = col; j < var + 1; ++j) {
          a[i][j] ^= a[row][j];
        }
      }
    }
  }
  for(int i = var - 1; i >= 0; --i) {
    x[i] = a[i][var];
    for(int j = i + 1; j < var; ++j) {
      x[i] ^= (a[i][j] & x[j]);
    }
  }
  return 0;
}

int main() {
  int T, cas = 1;
  scanf("%d", &T);
  while(T--) {
    memset(a, 0, sizeof(a));
    for(int i = 0; i < 30; ++i) {
      scanf("%d", &a[i][30]);
    }  
    for(int i = 0; i < 5; ++i) {
      for(int j = 0; j < 6; ++j) {
        int t = i * 6 + j;
        a[t][t] = 1;
        if(i > 0) a[(i - 1) * 6 + j][t] = 1;
        if(i < 4) a[(i + 1) * 6 + j][t] = 1;
        if(j > 0) a[i * 6 + j - 1][t] = 1;
        if(j < 5) a[i * 6 + j + 1][t] = 1;
      }
    }
    Gauss(30, 30);
    printf("PUZZLE #%d\n", cas++);
    for(int i = 0; i < 30; ++i) {
      printf("%d", x[i]);
      if((i + 1) % 6 == 0) printf("\n");
      else printf(" ");
    }
  }
  return 0;
}


POJ  1830 

开关问题

  


#include <iostream>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int MAXN = 284;
const double eps = 1e-6;

int a[MAXN][MAXN];
int x[MAXN];
int free_x[MAXN];

int Gauss(int equ, int var) {
  int k, max_r, col, ta, tb;
  for(int i = 0; i <= var; ++i) {
    x[i] = 0;
    free_x[i] = 1;
  }
  for(k = 0, col = 0; k < equ && col < var; ++k, ++col) {
    max_r = k;
    for(int i = k + 1; i < equ; ++i) {
      if(abs(a[i][col]) > abs(a[max_r][col])) max_r = i;
    }
    if(max_r != k) {
      for(int i = k; i < var + 1; ++i) swap(a[max_r][i], a[k][i]);
    }
    if(a[k][col] == 0) {
      k--;
      continue;
    }
    for(int i = k + 1; i < equ; ++i) {
      if(a[i][col] != 0) {
        for(int j = col; j < var + 1; ++j) {
          a[i][j] ^= a[k][j];
        }
      }
    }
  }
  for(int i = k; i < equ; ++i) {
    if(a[i][col] != 0) return -1;
  }
  return var - k;
}

int st[MAXN], en[MAXN];

int main() {
  int T;
  scanf("%d", &T);
  while(T--) {
    int n;
    scanf("%d", &n);
    for(int i = 0; i < n; ++i) scanf("%d", st + i);
    for(int i = 0; i < n; ++i) scanf("%d", en + i);
    memset(a, 0, sizeof(a));
    int u, v;
    while(scanf("%d %d", &u, &v), u || v) {
      a[v-1][u-1] = 1;
    }
    for(int i = 0; i < n; ++i) a[i][i] = 1;
    for(int i = 0; i < n; ++i) a[i][n] = st[i] ^ en[i];
    int ans = Gauss(n, n);
    if(ans == -1) {
      puts("Oh,it's impossible~!!");
    } else {
      printf("%d\n", 1 << ans);
    }
  }
  return 0;
}


POJ 2065 

SETI

  

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <map>

using namespace std;
const int MAXN = 108;
const int inf = 0x7fffffff;
int n, mod;

int a[MAXN][MAXN];
int x[MAXN];
bool free_x[MAXN];

char str[MAXN];

int quickPow(int a, int p, int mod) {
  int ans = 1;
  a %= mod;
  while(p) {
    if(p & 1) {
      ans = (ans * a % mod) % mod;
    }
    a = a * a % mod;
    p >>= 1;
  }
  return ans;
}

int gcd(int a, int b) {
  return b == 0 ? a : gcd(b, a % b);
}

int lcm(int a, int b) {
  return a / gcd(a, b) * b;
}

int Gauss(int equ, int var, int mod) {
  int k, max_r;
  int col;
  int ta, tb;
  int LCM;
  int tmp;
  int free_x_num;
  int free_index;

  for(int i = 0; i <= var; ++i) {
    x[i] = 0;
    free_x[i] = true;
  }
  col = 0;
  for(k = 0; k < equ && col < var; ++k, ++col) {
    max_r = k;
    for(int i = k + 1; i < equ; ++i) {
      if(abs(a[i][col]) < abs(a[max_r][col])) max_r = i;
    }
    if(max_r != k) {
      for(int j = k; j < var + 1; ++j) swap(a[k][j], a[max_r][j]);
    }
    if(a[k][col] == 0) {
      k--;
      continue;
    }
    for(int i = k + 1; i < equ; ++i) {
      if(a[i][col] != 0) {
        LCM = lcm(abs(a[i][col]), abs(a[k][col]));
        ta = LCM / abs(a[i][col]);
        tb = LCM / abs(a[k][col]);
        if(a[i][col] * a[k][col] < 0) tb = -tb;
        for(int j = col; j < var + 1; ++j) {
          a[i][j] = ((a[i][j] * ta - a[k][j] * tb) % mod + mod) % mod; 
        }
      }
    }
  }
  for(int i = var - 1; i >= 0; --i) {
    tmp = a[i][var];
    for(int j = i + 1; j < var; ++j) {
      if(a[i][j] != 0) tmp -= a[i][j] * x[j];
      tmp = (tmp % mod + mod) % mod;
    }
    while(tmp % a[i][i] != 0) tmp += mod;
    x[i] = (tmp / a[i][i]) % mod;
  }
  return 0;
}


int main() {
  int T;
  scanf("%d", &T);
  while(T--) {
    memset(a, 0, sizeof(a));
    scanf("%d %s", &mod, &str);
    n = strlen(str);
    for(int i = 0; i < n; ++i) {
      if(str[i] != '*') {
        a[i][n] = str[i] - 'a' + 1;
      } else {
        a[i][n] = 0;
      }
      for(int j = 0; j < n; ++j) {
        a[i][j] = quickPow(i + 1, j, mod);
      }
    }
    Gauss(n, n, mod);
    for(int i = 0; i < n; ++i) {
      if(i == n - 1) printf("%d\n", x[i]);
      else printf("%d ", x[i]);
    }
  }
  return 0;
}

POJ 2947 

Widget Factory

  

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <map>

using namespace std;
const int MAXN = 328;
const int mod = 7;
const int inf = 0x7fffffff;

int a[MAXN][MAXN];
int x[MAXN];
bool free_x[MAXN];

map<string, int> mp;

int gcd(int a, int b) {
  return b == 0 ? a : gcd(b, a % b);
}

int lcm(int a, int b) {
  return a / gcd(a, b) * b;
}

int Gauss(int equ, int var, int mod) {
  int i, j, k;
  int max_r, col;
  int ta, tb;
  int LCM, tmp;
  int free_x_num;
  int free_index;
  for(int i = 0; i <= var; ++i) {
    x[i] = 0;
    free_x[i] = true;
  }
  col = 0;
  for(k = 0; k < equ && col < var; k++, col++) {
    max_r = k;
    for(i = k + 1; i < equ; ++i) {
      if(abs(a[i][col]) > abs(a[max_r][col])) max_r = i;
    }
    if(max_r != k) {
      for(j = k; j < var + 1; ++j) swap(a[k][j], a[max_r][j]);
    }
    if(a[k][col] == 0) {
      k--;
      continue;
    }
    for(i = k + 1; i < equ; ++i) {
      if(a[i][col] != 0) {
        LCM = lcm(abs(a[i][col]), abs(a[k][col]));
        ta = LCM / abs(a[i][col]);
        tb = LCM / abs(a[k][col]);
        if(a[i][col] * a[k][col] < 0) tb = -tb;
        for(j = col; j < var + 1; ++j) {
          a[i][j] = ((a[i][j] * ta - a[k][j] * tb) % mod + mod) % mod;
        }
      }
    }
  }
  
  for(i = k; i < equ; ++i) {
    if(a[i][col] != 0) return -1;  // No Solution
  }
  if(k < var) return var - k;      // Infinity solutions
  for(int i = var - 1; i >= 0; --i) {
    tmp = a[i][var];
    for(int j = i + 1; j < var; ++j) {
      if(a[i][j] != 0) tmp -= a[i][j] * x[j];
      tmp = (tmp % mod + mod) % mod;
    }
    while(tmp % a[i][i] != 0) tmp += mod;
    x[i] = (tmp / a[i][i]) % mod;
  }
  return 0;
 }

int main() {
  mp["MON"] = 1;
  mp["TUE"] = 2;
  mp["WED"] = 3;
  mp["THU"] = 4;
  mp["FRI"] = 5;
  mp["SAT"] = 6;
  mp["SUN"] = 7;
  
  int n, m;
  char d1[20], d2[20];
  while(scanf("%d %d", &n, &m), n || m) {
    memset(a, 0, sizeof(a));
    int k;
    for(int i = 0; i < m; ++i) {
      scanf("%d%s%s", &k, &d1, &d2);
      a[i][n] = mp[d2] - mp[d1] + 1;
      // cerr << a[i][n] << endl;
      int id;
      while(k--) {
        scanf("%d", &id);
        a[i][id - 1]++;
        a[i][id - 1] %= 7; 
      }
    }
    // test starts
    /*
    cerr << " n = " << n << " m = " << m << endl;
    for(int i = 0; i < m; ++i) {
      for(int j = 0; j < n + 1; ++j) {
        cerr << a[i][j] << " ";
      }
      cout << endl;
    }
    */
    // test ends
    int ans = Gauss(m, n, mod);
    // test starts
    /*
    cerr << " ans = " << ans << endl;
    for(int i = 0; i < n; ++i) {
      cerr << x[i] << " ";
    }
    */
    // test ends
    if(ans == 0) {
      for(int i = 0; i < n; ++i) {
        if(x[i] <= 2) x[i] += 7;
      }
      for(int i = 0; i < n; ++i) {
        if(i == n - 1) printf("%d\n", x[i]);
        else printf("%d ", x[i]);
      }
    }
    else if(ans == -1) {
      puts("Inconsistent data.");
    } else { 
      puts("Multiple solutions.");  
    }  
  }
  return 0;
}

POJ 3185 

The Water Bowls

  

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>

using namespace std;
const int INF = 0x7fffffff;
const int MAXN = 27;

int a[MAXN][MAXN];
int x[MAXN];

//高斯消元法解模线性方程组
bool free_x[MAXN];//标记是否是不确定的变元
int gcd(int a, int b) {
  return b ? gcd(b, a % b) : a;
}

int lcm(int a, int b) {
  return a / gcd(a, b) * b;
}

// 高斯消元法解方程组
//-2表示有浮点数解,但无整数解
//-1表示无解
//0表示唯一解
//大于0表示无穷解，并返回自由变元的个数
//有equ个方程，var个变元
//增广矩阵行数为equ,分别为0到equ-1,列数为var+1,分别为0到var.
int Gauss(int equ, int var, int mod) {
  int i, j, k;
  int max_r;//当前这列绝对值最大的行.
  int col;//当前处理的列
  int ta, tb;
  int LCM;
  int temp;
  int free_x_num;
  int free_index;
  for (int i = 0; i <= var; i++) {
    x[i] = 0;
    free_x[i] = true;
  }
  //转换为阶梯阵.
  col = 0;//当前处理的列
  for (k = 0; k < equ && col < var; k++, col++) {
    // 枚举当前处理的行.
    // 找到该col列元素绝对值最大的那行与第k行交换(为了在除法时减小误差)
    max_r = k;
    for (i = k + 1; i < equ; i++) {
      if (abs(a[i][col]) > abs(a[max_r][col])) max_r = i;
    }
    if (max_r != k) {// 与第k行交换.
      for (j = k; j < var + 1; j++) swap(a[k][j], a[max_r][j]);
    }
    if (a[k][col] == 0) {// 说明该col列第k行以下全是0了,则处理当前行的下一列.
      k--;
      continue;
    }
    for (i = k + 1; i < equ; i++) {// 枚举要删去的行.
      if (a[i][col] != 0) {
        LCM = lcm(abs(a[i][col]), abs(a[k][col]));
        ta = LCM / abs(a[i][col]);
        tb = LCM / abs(a[k][col]);
        if (a[i][col] * a[k][col] < 0)tb = -tb;//异号的情况是相加
        for (j = col; j < var + 1; j++) {
          a[i][j] = ((a[i][j] * ta - a[k][j] * tb) % mod + mod) % mod;
        }
      }
    }
  }
  //无解的情况
  for (i = k; i < equ; i++) {
    if (a[i][col] != 0) return -1;
  }
  //找出和最小的解
  int ans = INF;
  for (int state = 0; state < (1 << (var - k)); state++) {
    int cnt = 0;
    for (int i = 0; i < var - k; i++)
      if (state & (1 << i))x[i + k] = 1, cnt++;
      else x[i + k] = 0;
    for (int i = k - 1; i >= 0; i--) {
      temp = a[i][var];
      for (j = i + 1; j < var; j++) {
        if (a[i][j] != 0) temp -= a[i][j] * x[j];
        temp = (temp % mod + mod) % mod;
      }
      while (temp % a[i][i] != 0) temp += mod;
      x[i] = (temp / a[i][i]) % mod;
      cnt += x[i];
    }
    ans = min(ans, cnt);
  }
  return ans;
}

int main() {
  int n = 20;
  memset(a, 0, sizeof(a));
  for (int i = 0; i < 20; ++i) {
    scanf("%d", &a[i][n]);
  }
  for (int i = 0; i < n; ++i) {
    a[i][i] = 1;
    if (i > 0) a[i - 1][i] = 1;
    if (i < n - 1) a[i + 1][i] = 1;
  }
  int ans = Gauss(n, n, 2);
  printf("%d\n", ans);
  return 0;
}


Southan97
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
高斯消元（非常详细） 高斯消元（非常详细）
06-19
高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元...
高斯消元附解方程神器
04-09
先输入一个数n表示方程元数。 然后再输入一个 n*(n+1)的矩阵 具体见我的博客
高斯消元
Nismilesucc的学习札记
08-20 2385
1.高斯消元 在O(n^3)下求解一个包含n个方程组n个未知数的多元线性方程组 n元线性方程组解的三种情况：无穷多组解、唯一解、无解 例题1 解线性方程组 输入 3 1.00 2.00 -1.00 -6.00 2.00 1.00 -3.00 -9.00 -1.00 -1.00 2.00 7.00 输出 1.00 -2.00 3.00 代码 例题2 解异或线性方程组 输入 3 1 1 0 1 0 1 1 0 1 0 0 1 输出 1 0 0 代码 ...
高斯消元专题训练
HeartWine的博客
11-04 457
文章目录POJ 1222 EXTENDED LIGHTS OUT POJ 1222 EXTENDED LIGHTS OUT 题意： 给定一个555行666列的010101矩阵，每一个小格子代表一个开关，000表示开关关闭，111表示开关闭合，当按下一个开关后，其也会影响该开关上下左右的开关的状态。现在给定一个这样的01矩阵，输出一种方案，使得所有开关都变成闭合状态。 思路：555行666列，相当于一共是303030个未知数，可以根据题意列出303030个方程，然后高斯消元，因为都是010101，所以是解异或
线性代数--高斯消元法
qq_51776588的博客
11-19 795
Widget Factory 模线性方程组
BZOJ 3503: [Cqoi2014]和谐矩阵（高斯消元）
少主大人殿下的博客
03-22 157
题目 给定一个矩形，要求你去填上0/1，其中每个格子周围4个格子+自己的1的个数位偶数，求一个方案数，不能全为0，有SPJ。 分析 暴力算法是一个裸的高斯消元，把每个元素都列出来，而且可以过。 进阶版的就是知道第一排元素之后，可以递推出后面的所有元素，如果m+1排元素为0，那么就是可行解。 把所有自由元弄成1可以防止出现全0的情况。 递推的时候可以用状压。 找可行解的方法是先算出...
高斯消元法
热门推荐
幸福在路上
01-10 1万+
高斯消元法(Gauss Elimination) 分析 & 题解 & 模板——czyuan原创 高斯消元法，是线性代数中的一个算法，可用来求解线性方程组，并可以求出矩阵的秩，以及求出可逆方阵的逆矩阵。 高斯消元法的原理是： 若用初等行变换将增广矩阵 化为 ，则AX = B与CX = D是同解方程组。 所以我们可以用初等行变换把增广矩阵转换为行阶梯阵，然后回代求出方程的解。
gaosixiaoyuanfa.rar_高斯消元
09-23
高斯消元法课件，将的很详细，对求解最小二乘法等算法有较大的帮助。
高斯消元模版深度认识完整版
06-19
高斯消元模版
高斯消元线性方程
09-26
这个例子中，`GaussianElimination`方法会执行上述的高斯消元步骤，而`BackSubstitution`方法则负责回代求解。当然，实际编程时还需要处理更复杂的边界条件和异常处理，以确保程序的健壮性。 总之，高斯消元法是...
线性方程组消去法求解 高斯消元 追赶法 数值计算作业
12-15
结合一个题目演示了如何自编代码实现高斯消元法和追赶法，各函数文件独立，便于移植，题目附有解答，题目来自西工大数值计算方法作业。
三种消元法（全主元、Gauss消去法、列主元）
04-23
三种消元法（全主元、Gauss消去法、列主元）三种消元法（全主元、Gauss消去法、列主元）三种消元法（全主元、Gauss消去法、列主元）三种消元法（全主元、Gauss消去法、列主元）
高斯消元学习总结+例题
JKdd123456的博客
07-29 4472
有具体的类型代码，先整理下来以后有空仔细分析 转大佬博客 整数类型高斯消元 返回值的情况 -2表示有浮点数解，但无整数解 -1表示无解 0表示唯一解 大于0表示无穷解，并返回自由变元的个数 其他说明 有equ个方程，var个变元。 增广矩阵行数为equ,分别为0到equ-1,列数为var+1,分别为0到var. #include &lt;algorithm&gt; #inclu...
高斯消元习题
cheng__yu_的博客
07-23 2176
高斯消元P3389 【模板】高斯消元法 （浮点数模板题） P3389 【模板】高斯消元法 （浮点数模板题） 链接：https://www.luogu.com.cn/problem/P3389 #include <bits/stdc++.h> #define ll long long using namespace std; const int maxn=110; int n; double a[maxn][maxn]; double x[maxn]; int guess(double a[]
【Luogu 3389】【模板】高斯消元法（用高斯消元求解线性方程组）
小哈里的博客
01-01 4625
1、线性方程组 2、高斯消元法 步骤1：通过一系列的加减消元，得到类似 kx=b 的式子，求得最后一个未知量的结果 步骤2：然后逐一回代求解整个 x 向量 以下列方程为例： 第一次加减消元，用第1式子消去后面所有的x得到： 第一次加减消元，用第2个式子消去后面所有的y得到： 这样就完成了高斯消元的步骤1，形成了一个倒三角形的形状，接下来主意回代即可。 3、用矩阵表示高斯消元 （1）消元过程： （2）无解： 当消元完毕后，发现有一行系数都为 0，但是常数项不为 0，此时无解 例如： （3）多解： 当
高斯消元解线性方程组
wmy0217_的博客
05-22 6169
文章目录基础知识算法步骤举个例子模拟一下算法例题和代码 基础知识 高斯消元可以通过初等行列变化把 增广矩阵 转换成 阶梯型矩阵，进而求解 n 个线性方程组的解，其时间复杂为 O(n^3^) 初等行列变换：(对一个方程组进行以下三个操作不会影响方程的解) 把某一行乘一个非0的数（方程左右两边同乘） 交换某两行的位置 把某一行的若干倍加到另一行上 例如线性方程组为： a11x1 + a12x2 + a13x3 + … + a1nxn = b1 a21x1 + a22x2 + a23x3 + … + a2nx
高斯消元(写(shui)题必备)
Coco_T的博客
10-07 1273
前言： 在一次学校hu测中， 遇到一道正解不用高斯消元，但是部分分需要的中档题 用舒老师的话说，只要是会高斯消元和树形dp 乱搞一下那道题就可以水到70 所以还是学习一下这个很有用算法：高斯消元简介 数学上，高斯消元法(或译:高斯消去法)，是线性代数中的一个算法， 可用来为线性方程组求解，求出矩阵的秩，以及求出可逆方阵的逆矩阵。 当用于一个矩阵时，高斯消元法会产生出一个“行梯阵式”。
[BZOJ4820] 硬币游戏 - 高斯消元
whzzt的博客
04-13 1588
不会捉题选手的瞎口胡
bzoj2155(高斯消元，线性基)
a1035719430的博客
09-20 363
Description Input 第一行包含两个整数N和 M， 表示该无向图中点的数目与边的数目。 接下来M 行描述 M 条边，每行三个整数Si，Ti ，Di，表示 Si 与Ti之间存在 一条权值为 Di的无向边。 图中可能有重边或自环。 Output 仅包含一个整数，表示最大的XOR和（十进制结果）,注意输出后加换行回车。 Sample Input 5 7 1 2 2 1 3 2 2 4 ...
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、付费专栏及课程。

余额充值