线性代数练习课程（hjyowl精讲练习题）

线性代数的基础在这儿呢：https://blog.csdn.net/hjyowl/article/details/140039350?spm=1001.2014.3001.5502

还记得上次的文章里面布置的习题吗？没看的同学赶紧看一遍：

P1939 矩阵加速（数列） - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

这道题呢，和上次讲解的很像，首先还是设置状态，如下：

P.S 由于直接写看着太乱了，放到代码板里
我们需要让矩阵    {f[i-2]} {f[i-1]}
             A * {f[i-1]}={f[i]}
                 {f[i]}   {f[i+1}
A[1][1]*f[i-2]+A[1][2]*f[i-1]+A[1][3]*f[i]=f[i-1]
容易得知只有A[1][2]=1,其他的都为0
A[2][1]*f[i-2]+A[2][2]*f[i-1]+A[2][3]*f[i]=f[i]
也容易知道A[2][3]=1,其他都为0
A[3][1]*f[i-2]+A[3][2]*f[i-1]+A[3][3]*f[i]=f[i+1]
根据递推式f[i]=f[i-1]+f[i-3],得知只有A[3][2]为0,其他的都是1.
然后矩阵是{1}
         {1}
         {1}
然后输出答案的第二行，因为他是f[n],不过，由于我们初始值多算了一项，所以说n要减一.

是不是很简单？有没有一种感觉很熟悉了的感觉？是不是觉得很套路？

我们来看一下今天讲解的另外一道题，非01初始矩阵的（我去的名字）：P1349 广义斐波那契数列 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

这道题呢，就是f[i]=p*f[i-1]+q*f[i-2],不过也很简单，如下： 

这次，我们横着来构造（因为竖着构造看着不好看）
{f[i-1] f[i]} * A = {f[i] f[i+1]}
A的第一列分别等于0,1第二列等于q,p(很容易想到)
所以说矩阵就是：
(0,q)
(1,p)
然后就可以了
初始者就是x1和x2,但是由于我们多算了一项，n需要-1

这是代码：

#include <bits/stdc++.h>
using namespace std;
long long mod = 1e9 + 7;
int p,q,x1,x2,n;
const long long N = 1010;
void mul_qmi1(long long f[3],long long a[3][3]){
	long long c[3];
	memset(c,0,sizeof c);
	for (long long i = 1; i <= 2; i ++ ){
		for (long long j = 1; j <= 2; j ++ ){
			c[i] = (c[i] + (long long)f[j] * a[j][i]) % mod;
		}
	}
	memcpy(f,c,sizeof c);
}
void mul_qmi2(long long a[3][3]){
	long long c[3][3];
	memset(c,0,sizeof c);
	for (long long i = 1; i <= 2; i ++ ){
		for (long long j = 1; j <= 2; j ++ ){
			for (long long k = 1; k <= 2; k ++ ){
				c[i][j] = (c[i][j] + (long long)a[i][k] * a[k][j]) % mod;
			}
		}
	}
	memcpy(a,c,sizeof c);
}
void qmi(){
	cin >> p >> q >> x1 >> x2 >> n >> mod;
	if (n == 1){
		cout << x1 % mod;
	}
	else if (n == 2){
		cout << x2 % mod;
	}
	else{
		long long f[3];
		long long a[3][3];
		f[1] = x1 % mod;
		f[2] = x2 % mod;
		a[1][1] = 0;
		a[1][2] = q % mod;
		a[2][1] = 1;
		a[2][2] = p % mod;
		n -= 1;
		while (n){
			if (n & 1){
				mul_qmi1(f,a);
			}
			mul_qmi2(a);
			n >>= 1;
		}
		cout << f[1] << endl;
	}
}
int main(){
	qmi();
	return 0;
}

好吧，很简单吧，而且有一个横着构造了的了，你可以竖着也可以横着，喜欢怎么样就怎么样。

习题到此结束，以后（不知道时间）可能会继续讲解矩阵的高斯消元，敬请关注。