题目
M斐波那契数列
Problem Description
M斐波那契数列F[n]是一种整数数列,它的定义如下:
F[0] = a
F[1] = b
F[n] = F[n-1] * F[n-2] ( n > 1 )
现在给出a, b, n,你能求出F[n]的值吗?
Input
输入包含多组测试数据;
每组数据占一行,包含3个整数a, b, n( 0 <= a, b, n <= 10^9 )
Output
对每组测试数据请输出一个整数F[n],由于F[n]可能很大,你只需输出F[n]对1000000007取模后的值即可,每组数据输出一行。
Sample Input
0 1 0 6 10 2
Sample Output
0 60
大体思路
刚开始的斐波那契数列是一个递归式,因增长速度过快,一般求到50左右就差不多了,随着渐渐学习,数据量的增加,可以通过矩阵快速幂求解,并且时间复杂度很低。看越多越发现,矩阵快速幂在这上面的应用有很多,各种式子的变形也有很多。我感觉这篇还是不错的。根据递推公式构造系数矩阵用于快速幂
分析:
首先根据题意我们可以得到F[0] = a 、 F[1] = b 、F[2] = a*b 、 F[3] = a * b^2 、 F[4} = a^2 * b^3、F[5] = a^3 * b^5… 最后推出式子为a ^ (f[i-2]) * b ^ (f[i-1]) % mod,拆开后a^f[i-2] % MOD * b^f[i-1] % mod。
如果数据量不大,这样直接矩阵快速幂就可以过了,但是现在数据量太大,我们可以通过和模运算相关的费马小定理来处理。费马小定理推导公式a^n = a ^ (n % (m-1)) (%m)。成立条件a是整数,m是质数且a不是m的倍数。
这有一篇关于快速幂取模运算的模板,这样就可以解决之前做过的一道幂运算取模越界的题了。快速幂取模运算 有的题也有可能会用到欧拉函数处理。
总结一下就是:
思路:
先用矩阵快速幂求解,再用费马小定理对求解出来的数值优化,最后快速幂求解最终答案。
代码还没来得及敲,注意费马小定理优化时模的是这个数的欧拉函数mod-1.