高斯消元专题训练

最新推荐文章于 2022-10-31 20:26:08 发布
最新推荐文章于 2022-10-31 20:26:08 发布
阅读量468 收藏 1
点赞数
文章标签： 矩阵 线性代数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_50815157/article/details/121106542
版权

文章目录

1、高斯消元求解开关问题

POJ 1222 EXTENDED LIGHTS OUT

题意：
给定一个 5 5 5行 6 6 6列的 01 01 01矩阵，每一个小格子代表一个开关， 0 0 0表示开关关闭， 1 1 1表示开关闭合，当按下一个开关后，其也会影响该开关上下左右的开关的状态。现在给定一个这样的01矩阵，输出一种方案，使得所有开关都变成闭合状态。

思路： 5 5 5行 6 6 6列，相当于一共是 30 30 30个未知数，可以根据题意列出 30 30 30个方程，然后高斯消元，因为都是 01 01 01，所以是解异或方程组。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <algorithm>
#define INF 0x3f3f3f3f

using namespace std;

typedef long long LL;

const int N = 35;

int n = 30, m;
int a[N][N], id[N][N];
int dx[] = {-1, 1, 0, 0}, dy[] = {0, 0, -1, 1};

void guass()
{
	int r, c;
	for(c = 0, r = 0; c < n; c ++)
	{
		int t = r; //找主元
		for(int i = r + 1; i < n; i ++)
			if(a[i][c])	{ t = i; break; }
		
		if(!a[t][c])	continue;//为0跳过

		for(int i = c; i <= n; 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; j ++)
					a[i][j] ^= a[r][j];

		r ++;
	}

	for(int i = r - 1; i >= 0; i --)//倒推求解
		for(int j = i + 1; j < n; j ++)
			a[i][n] ^= a[i][j] & a[j][n];
}

int main()
{
#ifdef LOCAL
	freopen("in.in", "r", stdin);
	freopen("out.out", "w", stdout);
#endif

	int t; cin >> t;
	for(int i = 0, k = 0; i < 5; i ++)
		for(int j = 0; j < 6; j ++)
			id[i][j] = k ++;
	
	int kase = 0;
	while(t --)
	{
		memset(a, 0, sizeof a);
		for(int i = 0; i < n; i ++)
		{
			cin >> a[i][n];
			a[i][i] = 1;
		}

		for(int i = 0; i < 5; i ++)
			for(int j = 0; j < 6; j ++)
				for(int k = 0; k < 4; k ++)
				{
					int ni = i + dx[k], nj = j + dy[k];
					if(ni < 0 || ni >= 5 || nj < 0 || nj >= 6)	continue;
					int u = id[i][j], v = id[ni][nj];
					a[u][v] = 1;
				}
		
		guass();

		printf("PUZZLE #%d\n", ++ kase);
		for(int i = 0; i < 5; i ++)
		{
			for(int j = 0; j < 6; j ++)	
				printf("%d ", a[id[i][j]][n]);
			putchar('\n');
		}
	}

	return 0;
}

POJ 1753 Flip Game

题意：
相当于给定一个 01 01 01矩阵，求全部变成0或则全部变成1的最小步数。

注意：
当求解方程时，出现了自由元的时候，要去枚举自由元的值。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <algorithm>
#define INF 0x3f3f3f3f

using namespace std;

typedef long long LL;

const int N = 20;

int n, a[N][N], b[N][N], cp[N][N];
int dx[] = {-1, 1, 0, 0}, dy[] = {0, 0, -1, 1};

int get(int x, int y)
{
	return 4 * x + y;
}

int guass(int a[][N])
{
	int r, c;
	for(r = 0, c = 0; c < n; c ++)
	{
		int t = r;
		for(int i = r + 1; i < n; i ++)
			if(a[i][c])	{ t = i; break; }

		if(!a[t][c])	continue;

		for(int i = c; i <= n; 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; j ++)
					a[i][j] ^= a[r][j];

		r ++;
	}

	if(r < n)
	{
		for(int i = r; i < n; i ++)
			if(a[i][n])	return INF;
	}

    //需要枚举变元的值
	int res = INF, t = n - r; //t是变元的个数
	
	for(int i = 0; i < 1 << t; i ++)
	{
	    int tmp = 0;
	    memcpy(cp, a, sizeof cp);
	    for(int j = 0; j < t; j ++)
	        if(i >> j & 1)  cp[n - j - 1][n] = 1, tmp ++;
	    
    	for(int j = r - 1; j >= 0; j --)
    	{
    		for(int k = j + 1; k < n; k ++)
    			cp[j][n] ^= cp[j][k] & cp[k][n];		
    		tmp += cp[j][n];
    	}
    	//cout << tmp << endl;
    	res = min(res, tmp);
	}

	return res;
}

int main()
{
#ifdef LOCAL
	freopen("in.in", "r", stdin);
	freopen("out.out", "w", stdout);
#endif

	n = 16;
	for(int i = 0; i < 4; i ++)
		for(int j = 0; j < 4; j ++)
		{
			char ch; cin >> ch;
			int id = get(i, j);
			if(ch == 'b')	a[id][n] = 1;
			else	b[id][n] = 1;

			a[id][id] = b[id][id] = 1;
		}

	for(int i = 0; i < 4; i ++)
		for(int j = 0; j < 4; j ++)
			for(int k = 0; k < 4; k ++)
			{
				int ni = i + dx[k], nj = j + dy[k];
				if(ni < 0 || ni >= 4 || nj < 0 || nj >= 4)	continue;
				int u = get(i, j), v = get(ni, nj);
				a[u][v] = b[u][v] = 1;
			}
	
	int ans = min(guass(a), guass(b));

	if(ans == INF)	puts("Impossible");
	else	cout << ans << endl;

	return 0;
}

P2962 [USACO09NOV]Lights G

在这里插入图片描述
在这里插入图片描述

本题写法与上面有些不同，因为自由元可能较多，上面的枚举方法肯定会TLE，所以采用dfs+剪枝枚举。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <algorithm>
#define INF 0x3f3f3f3f

using namespace std;

typedef long long LL;

const int N = 55;

int n, m, a[N][N], x[N];
int ans = INF;

void gauss()
{//这里的写法没有将所有无解或无穷解的行移到最下面，方便dfs枚举
//保证了第i行一定是第i个未知数的方程
	for(int c = 0; c < n; c ++)
	{
		int t = c;
		for(int i = c; i < n; i ++)
			if(a[i][c])	{ t = i; break; }
		
		if(!a[t][c])	continue;

		for(int i = c; i <= n; i ++)	swap(a[t][i], a[c][i]);
		for(int i = c + 1; i < n; i ++)
			if(a[i][c])
				for(int j = c; j <= n; j ++)
					a[i][j] ^= a[c][j];
	}
}

void dfs(int k, int tot)
{
	if(tot >= ans)	return ; //剪枝
	if(k < 0)	{ ans = tot; return ; }

	if(a[k][k]) //不是自由元，正常倒推求解
	{
		x[k] = a[k][n];
		for(int i = k + 1; i < n; i ++)
			x[k] ^= a[k][i] & x[i];
		dfs(k - 1, tot + x[k]);
	}
	else//是自由元，枚举取值
	{
		x[k] = 0; dfs(k - 1, tot);
		x[k] = 1; dfs(k - 1, tot + 1);
	}
}

int main()
{
#ifdef LOCAL
	freopen("in.in", "r", stdin);
	freopen("out.out", "w", stdout);
#endif

	cin >> n >> m;
	while(m --)
	{
		int u, v;
		cin >> u >> v;
		u --, v --;
		a[u][v] = a[v][u] = 1;
	}
	for(int i = 0; i < n; i ++)	a[i][i] = a[i][n] = 1;

	gauss();
	dfs(n - 1, 0);
	cout << ans << endl;

	return 0;
}

2、高斯消元求解线性方程组

P3389 【模板】高斯消元法

在这里插入图片描述

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <algorithm>
#define INF 0x3f3f3f3f

using namespace std;

typedef long long LL;

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

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

int guass()
{
	int r, c;
	for(r = 0, c = 0; c < n; c ++)
	{
		int t = r;//找主元
		for(int i = r + 1; i < n; i ++)
			if(fabs(a[i][c]) > fabs(a[t][c]))	t = i;

		if(fabs(a[t][c]) < eps)	continue;//为0跳过

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

		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[i][c] * a[r][j];
		r ++;
	}

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

	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];
}

int main()
{
#ifdef LOCAL
	freopen("in.in", "r", stdin);
	freopen("out.out", "w", stdout);
#endif

	cin >> n;
	for(int i = 0; i < n; i ++)
		for(int j = 0; j <= n; j ++)	cin >> a[i][j];
	
	int t = guass();

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

	return 0;
}

P2973 [USACO10HOL]Driving Out the Piggies G

在这里插入图片描述

做法:
设 f i f_i fi​表示到达i点的概率，则在i点爆炸的概率就是 p q ∗ f i \frac{p}{q}*f_i qp​∗fi​

f i f_i fi​的递推公式可以表示为 f i = ∑ j → i f j ∗ ( 1 − p q ) ∗ 1 d e g r e e j f_i=\sum_{j\to i}f_j*(1-\frac{p}{q})*\frac{1}{degree_j} fi​=∑j→i​fj​∗(1−qp​)∗degreej​1​，其中 d e g r e e j degree_j degreej​是 j j j的的度数。

但是本题无法直接 d p dp dp，本图是无向图，所以这个递推方程是有环的，有后效性。

我们可以将每一个 f i f_i fi​看成一个未知数，这样就有 n n n个未知数，根据递推公式对于每一个 i i i都可以列出一个方程，相当于有了 n n n个方程，所以我们可以直接高斯消元求解。

要注意一点的是：本题的起点是 1 1 1号点，在列方程的时候要考虑 1 1 1号点的特殊性。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <algorithm>
#define INF 0x3f3f3f3f

using namespace std;

typedef long long LL;

const int N = 310;
const double eps = 1e-9;

int n, m;
double p, q, a[N][N];
int degree[N], g[N][N];

void guass()
{
	int r, c;
	for(r = 0, c = 0; c < n; c ++)
	{
		int t = r;
		for(int i = r + 1; 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; 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[i][c] * a[r][j];
		r ++;
	}

	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];
}

int main()
{
#ifdef LOCAL
	freopen("in.in", "r", stdin);
	freopen("out.out", "w", stdout);
#endif

	cin >> n >> m >> p >> q;
	p = 1.0 * p / q;
	while(m --)
	{
		int u, v;
		cin >> u >> v;
		u --, v --;
		g[u][v] = g[v][u] = 1;
		degree[u] ++, degree[v] ++;
	}

	for(int i = 0; i < n; i ++)
	{
		a[i][i] = 1.0;
		for(int j = 0; j < n; j ++)	
			if(g[j][i])	a[i][j] = -(1.0 - p) / degree[j];
	}
	a[0][n] = 1.0;//因为原本就在1号点，初值设为1

	guass();

	for(int i = 0; i < n; i ++)
		printf("%.9lf\n", a[i][n] * p);

	return 0;
}
proudzhao
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
高斯消元（非常详细） 高斯消元（非常详细）
06-19
高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元...
高斯消元附解方程神器
04-09
先输入一个数n表示方程元数。 然后再输入一个 n*(n+1)的矩阵 具体见我的博客
高斯消元
Nismilesucc的学习札记
08-20 2403
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 代码 ...
【Luogu 3389】【模板】高斯消元法（用高斯消元求解线性方程组）
小哈里的博客
01-01 4631
1、线性方程组 2、高斯消元法 步骤1：通过一系列的加减消元，得到类似 kx=b 的式子，求得最后一个未知量的结果 步骤2：然后逐一回代求解整个 x 向量 以下列方程为例： 第一次加减消元，用第1式子消去后面所有的x得到： 第一次加减消元，用第2个式子消去后面所有的y得到： 这样就完成了高斯消元的步骤1，形成了一个倒三角形的形状，接下来主意回代即可。 3、用矩阵表示高斯消元 （1）消元过程： （2）无解： 当消元完毕后，发现有一行系数都为 0，但是常数项不为 0，此时无解 例如： （3）多解： 当
第十九篇，解析法求解五阶多项式
qq_42466012的博客
10-31 662
五阶多项式解析解
线性代数--高斯消元法
qq_51776588的博客
11-19 799
Widget Factory 模线性方程组
数值计算(九)——线性代数方程组求解(一)高斯消元法
DFS的博客
03-30 2426
线性方程组求解 许多物理问题的数学模型最终都要归结为一个求解线性代数方程组的问题，因此掌握快速准确的求解线性代数方程组是极具必要性的，线性代数方程组的求解一般分为两大类别，直接求解和迭代求解两种方式，直接法是指在没有舍入误差的情况下经过有限次运算可求得方程组的精确解的方法，迭代法则是采取逐步逼近的方式，即从一个初始近似解出发，按某种迭代格式，逐步地向前推进，使其近似解逐步接近精确解，直到满足精度要求为止。 这时候长得比较帅的读者就会有疑问了，既然有直接求解方式，那为何还有搞一个迭代求解，求出来的结果还不
gaosixiaoyuanfa.rar_高斯消元
09-23
这个名为“gaosixiaoyuanfa.rar_高斯消元”的压缩包文件包含了一份详细的高斯消元法课程资料，通常以PPT的形式呈现，可以帮助学习者深入理解并掌握这一算法。 高斯消元法的基本思想是通过一系列行初等变换，将系数...
高斯消元模版深度认识完整版
06-19
高斯消元模版
gauss_seq.c.zip_高斯消元
09-23
此外，高斯-约旦消元（Gauss-Jordan elimination）是高斯消元的一种扩展，它不仅将矩阵转化为上三角形，还能进一步将其化为单位矩阵，从而直接得到解。 总结来说，`gauss_seq.c`是一个用C语言实现的高斯消元法串行...
三种消元法（全主元、Gauss消去法、列主元）
04-23
三种消元法（全主元、Gauss消去法、列主元）三种消元法（全主元、Gauss消去法、列主元）三种消元法（全主元、Gauss消去法、列主元）三种消元法（全主元、Gauss消去法、列主元）
线性代数-高斯消元法
02-19
详细讲解高斯消元法，讲解线性方程组的一般解法，引入矩阵初等变换这一工具
线性方程组消去法求解 高斯消元 追赶法 数值计算作业
12-15
结合一个题目演示了如何自编代码实现高斯消元法和追赶法，各函数文件独立，便于移植，题目附有解答，题目来自西工大数值计算方法作业。
matlab实现高斯消元法
06-01
matlab实现高斯消元法，测试通过，保证能用
NOIP ACM 高斯消元练习+数据题解
04-01
NOIP ACM 高斯消元练习+数据题解
高斯消元法
zhendong
01-21 677
在数学模块，我们将学习矩阵相关知识，主要是矩阵乘法、高斯消元、线性空间。高斯消元也可以与其他知识点结合起来，比如dp。坏掉的机器人这道题目就是dp套高斯消元。 高斯消元的思路参考进阶指南。主要目的是将矩阵通过初等行变换最后转化为简化阶梯型矩阵。系数和常数构成n∗(n+1)n*(n+1)n∗(n+1)的增广矩阵。有三种特殊情况： 最后有n个主元。那么解唯一。 有k个主元，自由元n-k。解无限多。 存在0=d0=d0=d这种形式。无解。 下面举四个例题。 洛谷： P3389 【模板】高斯消元法 #inclu
高斯消元习题
cheng__yu_的博客
07-23 2190
高斯消元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[]
高斯消元学习总结+例题
JKdd123456的博客
07-29 4477
有具体的类型代码，先整理下来以后有空仔细分析 转大佬博客 整数类型高斯消元 返回值的情况 -2表示有浮点数解，但无整数解 -1表示无解 0表示唯一解 大于0表示无穷解，并返回自由变元的个数 其他说明 有equ个方程，var个变元。 增广矩阵行数为equ,分别为0到equ-1,列数为var+1,分别为0到var. #include &lt;algorithm&gt; #inclu...
高斯消元解线性方程组
wmy0217_的博客
05-22 6173
文章目录基础知识算法步骤举个例子模拟一下算法例题和代码 基础知识 高斯消元可以通过初等行列变化把 增广矩阵 转换成 阶梯型矩阵，进而求解 n 个线性方程组的解，其时间复杂为 O(n^3^) 初等行列变换：(对一个方程组进行以下三个操作不会影响方程的解) 把某一行乘一个非0的数（方程左右两边同乘） 交换某两行的位置 把某一行的若干倍加到另一行上 例如线性方程组为： a11x1 + a12x2 + a13x3 + … + a1nxn = b1 a21x1 + a22x2 + a23x3 + … + a2nx
matlab 高斯消元
最新发布
11-12
高斯消元是一种线性代数中常用的求解线性方程组的方法，Matlab中也提供了相应的函数来实现高斯消元。下面是一个简单的例子： 假设有如下线性方程组： ``` 2x + y - z = 8 -3x - y + 2z = -11 -2x + y + 2z = -3 ``` 可以将其表示为增广矩阵的形式： ``` [2 1 -1 8; -3 -1 2 -11; -2 1 2 -3] ``` 在Matlab中，可以使用`rref`函数来进行高斯消元，该函数返回一个阶梯形矩阵和一个列向量，其中阶梯形矩阵是增广矩阵的等价形式，列向量则是方程组的解。具体操作如下： ```matlab A = [2 1 -1; -3 -1 2; -2 1 2]; b = [8; -11; -3]; Ab = [A b]; rref(Ab) ``` 运行结果为： ``` ans = 1.0000 0 0 2.0000 0 1.0000 0 3.0000 0 0 1.0000 1.0000 ``` 其中，阶梯形矩阵为： ``` [1 0 0; 0 1 0; 0 0 1] ``` 列向量为： ``` [2; 3; 1] ``` 即方程组的解为`x=2, y=3, z=1`。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值