线性代数的基础在这儿呢: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;
}
好吧,很简单吧,而且有一个横着构造了的了,你可以竖着也可以横着,喜欢怎么样就怎么样。
习题到此结束,以后(不知道时间)可能会继续讲解矩阵的高斯消元,敬请关注。