全网最最最最最详细的c++算法讲解（三）矩阵乘法进阶与高斯消元法

建议各位在看本篇博客之前，先看一下我的博客的算法讲解（二），了解一些关于矩阵乘法和矩阵的相关知识

这里先来一道矩阵乘法好题，前几天看到的

设n元向量x=[x1,x2,…,xn]，y=[y1,y2,…,yn]，若A=yTx，求A^k      (yT表示矩阵y的转置）

n的范围是n<=100000

一看到，这个题，肯定很多人会想：”这不是个sb题吗？直接做不就好了？”

但其实不然，如果我们直接进行矩阵快速幂，复杂度是n^2*logk的

显然跑不过。

那么我们这时候就可以，动动脑筋了。 好像在前一篇文章里，我提到过这么一个结论

只要当A矩阵为x*y，B矩阵的乘法y*z时，乘出的C矩阵是一个x*z的矩阵！！！

那么，我们要是可以将x和yT相乘 就好了哎

展开式子，根据矩阵乘法的结合律，我们就会惊奇的发现：

Ak=(yTx)(yTx)…(yTx)
    =yT(xyT)(xyT)…(xyT)x

    =^k-1A

woc 骚操作呀！时间复杂度直接O(n^2 log k)-->O(log k+n^2)

社会社会

一道华丽的分割线------------------------------------

好啦，进入正题。

原谅我的懒.....

对于一个Ax=b的线性方程组

若b=0，则称Ax=b为齐次线性方程组

若b≠0，则称Ax=b为非齐次线性方程组

在解方程之前，我们先介绍一下矩阵的初等行变换

矩阵的三种初等行变换
(1)将A的第i行与第j行对换
(2)将A的第i行乘以非零常数k
(3)将A第i行的k倍加到第j行

并且并且并且------矩阵的初等行变换不改变线性方程组的解！

对矩阵A做一次初等行变换，等价于，在A的左边乘一个初等矩阵

这是矩阵行变换中常用的三个初等矩阵

Eij                  对换第i行与第j行
Ei(k)              第i行乘以k
Eij(k)             第i行的k倍加到第j行

在高斯消元中，我们普遍是将矩阵化为最简阶梯型矩阵：

阶梯形矩阵

若矩阵的每个非零行的上方没有零行，且各个非零行自左向右的第一个非零元素a1j1,a2j2,…,arjr所在列的编号满足j1

若阶梯形矩阵中个非零行自左向右第一个非零元素均为1，他们所在列的其余元素均为0，则称该阶梯形矩阵为最简阶梯形矩阵。

任意非零矩阵A只经初等行变换可化为阶梯形矩阵。
该阶梯形矩阵中主元列的个数称为矩阵A的秩，记作rank(A)

那么我们该如何求矩阵的秩？           化为阶梯形矩阵！

如何求方程组的解？                        化为最简洁阶梯型矩阵！

如何只通过行变换将任意矩阵化为阶梯形矩阵？   Gauss消元法！！！

大概步骤就是

i：1~n      k：当前已经选中的行数（矩阵的秩）
①挑选未被选中的某一行，满足这一行的第i列不为0
(若所有未被选择的行中，第i列都为0，则第i列不为主元列)
②将这一行换到第k+1行

③将这一行的若干倍加到其它未被选中的行上，使得它们的第i列变成0

最后求解即可。

这里注意

齐次方程组有解的条件：

若rank(A)=n，则只有零解
若rank(A)=r且非零解有n-r个自由变量

基础解系有n-r个解

非齐次方程组有解的条件：

若rank(A)≠rank([A,b])，则非齐次方程组无解
若rank(A)=rank([A,b])=n，则非齐次方程组有唯一解解

若rank(A)=rank([A,b])

我们需要在求解之前判断解的情况！！！！

上代码

bool gsxy()
{
	int k=1;
	for (int i=1;i<=n;i++)
	{
		int now = k; 
		while (now<=n && !a[now][i]) now++;
		if (now == n+1) continue;
		for (int j=1;j<=n+1;j++) swap(a[now][j],a[k][j]);
		for (int j=1;j<=n;j++)
		{
			if (j!=k){
				double t = a[j][i]/a[k][i];
				for (int p=1;p<=n+1;p++) a[j][p]=a[j][p]-a[k][p]*t;
			}
		}
		k++;
	}
	if (a[k][n+1]) return false;
	if (k-1

那么对于一些对解取余的题

我就在除法上，就需要逆元了

void gaosixiaoyuan()
{
	int k=1;
	for (int i=1;i<=n;i++)
	{
		int now=k;
		while (now<=n && !a[now][i]) now++;
		if (now==n+1) continue;
		for (int j=1;j<=n+1;j++) 
		  swap(a[now][j],a[k][j]);
		int inv=power(a[k][i],mod-2);
		for (int j=1;j<=n;j++)
		  if (j!=k)
		  {
		  	int t=(long long)*a[j][i]*inv%mod;
		  	for (int p=1;p<=n+1;p++) 
			  a[j][p]=(a[j][p]-(long long)*t*a[k][p]%mod+mod)%mod;
		  }
		k++;
	}
	if (a[k][n+1]) return false;
	if (k-1

差别不大 主要是逆元

啦啦啦 高斯消元的讲解就到这了~

希望大家能懂 嗯 不懂可以加我qq 752742355 有事情随意交流

如果您们能从我的博客中得到任何一点收获，都是我莫大的荣幸

共勉！

from y_immortal