前几天的题了
以前没做过公式推矩阵快速幂的题,讲道理,我学东西真的是 惰性太强,学的很慢,还好,相对来说比较扎实,但是应该追求又快又好啊
然后场上见了这个题,完全没有思路,后来想到了 斐波那契数列 的矩阵快速幂,1小时后推出来了
公式推矩阵快速幂理解:
首先能用到矩阵快速幂的,是 一个数的 n 次方 对 mod 取模的题目,矩阵的话 就是这个数转化成了 矩阵而已
所以,要推矩阵快速幂 要推出一个 固定的矩阵,然后还要有一个初始矩阵,也就是初值,然后跟那个固定矩阵相乘多少次 可以得到答案,
也就到了问题的关键,就是要由一个公式 推出 第 n-1 项和第 n 项之间的关系(这个关系更体现在 共性的点之间的关系),也就是 一个递推的关系式
然后通过矩阵的运算,保证原始矩阵的第一行 都是关键元素,(这样说可能不好理解,下面会由这个题来解释一下)
******相如这个题******
固定矩阵:
1 1 0 0 0 0 0
2 0 0 0 0 0 0
1 0 1 0 0 0 0
4 0 4 1 0 0 0
6 0 6 3 1 0 0
4 0 4 3 2 1 0
1 0 1 1 1 1 1
初始矩阵:
f2 f1 16 8 4 2 1
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
我们可以由原式: F ( n ) = n ^ 4 + 2 * F( n - 2 ) + 1 * F( n - 1 )
推出 : F ( n ) = 1 * F( n - 1 ) + 2 * F( n - 2 ) + 1 * (n-1)^4 + 4 * (n-1)^3 + 6 * (n-1)^2 + 4*(n-1)^1 + 1 * 1 ;
其实 固定矩阵是由原始矩阵推出来了,知道原始矩阵的第一行就好了,
一般我们 就用 原始矩阵 (0,0)位置的元素作为答案,然后保证 原始矩阵对固定矩阵运算后,第一行的各元素 意义和开始的时候一样就行(保证下次运算可以的到新的答案)
----------------------这个看不懂 可以先自己手推一遍 fibonaci 的矩阵快速幂
#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
typedef long long ll;
const ll mod = 2147493647;
int id;
ll f1, f2;
struct Matrix {
ll a[10][10];
void setMatrix(ll v) {
for(int i = 0; i < 7; ++i)
for(int j = 0; j < 7; ++j)
a[i][j] = v;
}
Matrix operator * (const Matrix& m2) {
Matrix t;
for(int i = 0; i < 7; ++i) {
for(int j = 0; j < 7; ++j) {
t.a[i][j] = 0;
for(int k = 0; k < 7; ++k) {
t.a[i][j] += (a[i][k] * m2.a[k][j]) % mod;
t.a[i][j] %= mod;
}
}
}
return t;
}
}A, t, E;
void init() {
E.setMatrix(0);
for(int i = 0; i < 7; ++i) E.a[i][i] = 1;
A.setMatrix(0); t.setMatrix(0);
t.a[0][0] = 1; t.a[0][1] = 1;
t.a[1][0] = 2;
t.a[2][0] = 1; t.a[2][2] = 1;
t.a[3][0] = 4; t.a[3][2] = 4; t.a[3][3] = 1;
t.a[4][0] = 6; t.a[4][2] = 6; t.a[4][3] = 3; t.a[4][4] = 1;
t.a[5][0] = 4; t.a[5][2] = 4; t.a[5][3] = 3; t.a[5][4] = 2; t.a[5][5] = 1;
t.a[6][0] = 1; t.a[6][2] = 1; t.a[6][3] = 1; t.a[6][4] = 1; t.a[6][5] = 1; t.a[6][6] = 1;
scanf("%d %lld %lld", &id, &f1, &f2);
if(id >= 3) { A.a[0][0] = f2%mod; A.a[0][1] = f1%mod; A.a[0][2] = 16; A.a[0][3] = 8; A.a[0][4] = 4; A.a[0][5] = 2; A.a[0][6] = 1; }
}
Matrix fast_pow(Matrix t, int n) {
Matrix ans = E;
while(n) {
if(n&1) ans = ans * t;
t = t * t;
n /= 2;
}
return ans;
}
void solve() {
if(id == 1) { cout << f1 << endl; return; }
if(id == 2) { cout << f2 << endl; return; }
t = fast_pow(t, id-2);
A = A * t;
cout << A.a[0][0] << endl;
}
int main() {
int T;
scanf("%d", &T);
while(T--) {
init();
solve();
}
return 0;
}