题解
大佬的完整题解 - 以及解题的关键步骤:特征法
图片来自大佬博客↓
之前一直想不通,为什么等式 ⌈ ( a + b ) n ⌉ = ( a + b ) n + ( a − b ) n \left \lceil (a+\sqrt b)^n\right \rceil=(a+\sqrt b)^n+(a-\sqrt b)^n ⌈(a+b)n⌉=(a+b)n+(a−b)n 成立,它不会超吗?现在有点想明白了,
假设 ( a + b ) n (a+\sqrt b)^n (a+b)n 的小数部分是小于0.5的, ( a − b ) n (a-\sqrt b)^n (a−b)n的小数部分是大于0.5的,小数部分相加大于1,
但这是不可能的,因为无论n是奇数还是偶数,公式 ( a + b ) n + ( a − b ) n (a+\sqrt b)^n+(a-\sqrt b)^n (a+b)n+(a−b)n 都是整数,
而如果为了凑整,使得小数部分相加等于2,这也是不可能的,因为 0 < ( a − b ) n < 1 0<(a-\sqrt b)^n<1 0<(a−b)n<1,
故,该记一记了 ⌈ ( a + b ) n ⌉ = ( a + b ) n + ( a − b ) n \left \lceil (a+\sqrt b)^n\right \rceil=(a+\sqrt b)^n+(a-\sqrt b)^n ⌈(a+b)n⌉=(a+b)n+(a−b)n
斐波那契数列通项公式推导及其类似构造共轭公式反求递推式
当有 f n = ( a + b ) n + ( a − b ) n f_n=(a+\sqrt b)^n+(a-\sqrt b )^n fn=(a+b)n+(a−b)n 时
存在递推公式: f n = 2 a ∗ f n − 1 + ( b − a 2 ) ∗ f n − 2 f_n=2a*f_{n-1}+(b-a^2)*f_{n-2} fn=2a∗fn−1+(b−a2)∗fn−2
所以,如果是求某个 f n f_n fn 确切的值时,要么找循环节,要么用矩阵快速幂
矩阵快速幂的转移公式:
[ f 2 f 1 ] = [ 2 a b − a 2 1 0 ] ⋅ [ f 1 f 0 ] \begin{bmatrix} f_2 \\ f_1\end{bmatrix}=\begin{bmatrix} 2a & b-a^2 \\ 1 & 0\end{bmatrix}·\begin{bmatrix} f_1\\ f_0 \end{bmatrix} [f2f1]=[2a1b−a20]⋅[f1f0]
[ f n f n − 1 ] = [ 2 a b − a 2 1 0 ] n − 1 ⋅ [ f 1 f 0 ] \begin{bmatrix} f_n \\ f_{n-1}\end{bmatrix}=\begin{bmatrix} 2a & b-a^2 \\ 1 & 0\end{bmatrix}^{n-1} ·\begin{bmatrix} f_1\\ f_0 \end{bmatrix} [fnfn−1]=[2a1b−a20]n−1⋅[f1f0]
再手动算出 f 1 f_1 f1和 f 0 f_0 f0,就可以了
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1e6 + 10;
const int INF = 0x3f3f3f3f;
int n, m, k;
int a, b;
namespace MAT {//矩阵快速幂板子
const int maxn = 2;
struct Mat {
ll mat[maxn][maxn];
void clear() {
memset(mat, 0, sizeof(mat));
}
Mat operator*(const Mat a) const {
Mat b;
b.clear();
for (int i = 0; i < maxn; ++i) {
for (int j = 0; j < maxn; ++j) {
for (int k = 0; k < maxn; ++k) {
b.mat[i][j] = (b.mat[i][j] + (mat[i][k] * a.mat[k][j]) % m + m) % m;
}
}
}
return b;
}
};
Mat pow(Mat m, ll k) {
Mat res;
res.clear();
for (int i = 0; i < maxn; ++i) {//单位矩阵
res.mat[i][i] = 1;
}
while (k) {
if (k & 1)res = res * m;
k >>= 1;
m = m * m;
}
return res;
}
Mat E;//转移矩阵
void init() {
E.clear();
/*
fn = a*f_n-1 + b*f_n-2
| f |=|a b|*|f_n-1|
|f_n-1| |1 0| |f_n-2|
需要一个转移矩阵 E = {a, b , 1, 0};
以及最开始的两个f1 f0
Mat f = {1, 0, 1, 0};// f1 f0
*/
E = {2 * a, b - a * a, 1, 0};
}
}
using namespace MAT;
int main() {
ios::sync_with_stdio(0);
while (cin >> a >> b >> n >> m) {
init();
Mat f = {2 * a, 0, 1, 0};// f1 f0
f = pow(E, n - 1) * f;
cout << f.mat[0][0] << endl;
}
return 0;
}