题目1 : 骨牌覆盖问题·一
时间限制:
10000ms
单点时限:
1000ms
内存限制:
256MB
-
62247088
样例输出
-
17748018
描述
骨牌,一种古老的玩具。今天我们要研究的是骨牌的覆盖问题:
我们有一个2xN的长条形棋盘,然后用1x2的骨牌去覆盖整个棋盘。对于这个棋盘,一共有多少种不同的覆盖方法呢?
举个例子,对于长度为1到3的棋盘,我们有下面几种覆盖方式:
输入
第1行:1个整数N。表示棋盘长度。1≤N≤100,000,000
输出
第1行:1个整数,表示覆盖方案数 MOD 19999997
们考虑在已经放置了部分骨牌(灰色)的情况下,下一步可以如何放置新的骨牌(蓝色):
那么通过对上面的观察,我们可以发现:
在任何一个放置方案最后,一定满足前面两种情况。而灰色的部分又正好对应了长度为N-1和N-2时的放置方案。由此,我们可以得到递推公式:
f[n] = f[n-1] + f[n-2];
这个公式是不是看上去很眼熟?没错,这正是我们的费波拉契数列。
f[0]=1,f[1]=1,f[2]=2,...
如何快速计算结果?
当N很小的时候,我们直接通过递推公式便可以计算。当N很大的时候,只要我们的电脑足够好,我们仍然可以直接通过递推公式来计算。
但是我们学算法的,总是这样直接枚举不是显得很Low么,所以我们要用一个好的算法来加速(装X)。
事实上,对于这种线性递推式,我们可以用矩阵乘法来求第n项。对于本题Fibonacci数列,我们希望找到一个2x2的矩阵M,使得(a, b) x M = (b, a+b),其中(a, b)和(b, a+b)都是1x2的矩阵。
显然,只需要取M = [0, 1; 1, 1]就可以了:
1. 先计算出所有的{a^1, a^2, a^4 ... a^(2^j)},因为该数列满足递推公式,时间复杂度为O(logN)
2. 将指数n二进制化,再利用公式将对应的a^j相乘计算出a^n,时间复杂度仍然为O(logN)
则总的时间复杂度为O(logN)
这种算法因为能够在很短时间内求出幂,我们称之为“快速幂”算法。
代码:
#include <iostream>
#include <cstdio>
using namespace std;
typedef long long ll;
const int MOD = 19999997;
ll a[2][2] = {{0, 1}, {1, 1}}, fa[2][2];
ll power_mod(ll b) {
ll ans[2] = {0, 1}, fans[2];
while(b>0) {
if(b%2 == 1) {
fans[0] = (ans[0]*a[0][0] + ans[1]*a[1][0])%MOD;
fans[1] = (ans[0]*a[0][1] + ans[1]*a[1][1])%MOD;
ans[0] = fans[0];
ans[1] = fans[1];
}
b /= 2;
fa[0][0] = (a[0][0]*a[0][0] + a[0][1]*a[1][0])%MOD;
fa[0][1] = (a[0][0]*a[0][1] + a[0][1]*a[1][1])%MOD;
fa[1][0] = (a[1][0]*a[0][0] + a[1][1]*a[1][0])%MOD;
fa[1][1] = (a[1][0]*a[0][1] + a[1][1]*a[1][1])%MOD;
a[0][0] = fa[0][0];
a[0][1] = fa[0][1];
a[1][0] = fa[1][0];
a[1][1] = fa[1][1];
}
return ans[1];
}
int main() {
ll n;
scanf("%lld", &n);
printf("%lld\n", power_mod(n));
return 0;
}