hdoj 4549 M斐波那契数列 【矩阵快速幂 + 费马小定理 + 快速幂】

题目

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.