Quad Tiling
Time Limit: 1000MS | Memory Limit: 65536K | |
Total Submissions: 1801 | Accepted: 698 |
Description
Tired of the Tri Tiling game finally, Michael turns to a more challengeable game, Quad Tiling:
In how many ways can you tile a 4 × N (1 ≤ N ≤ 109) rectangle with 2 × 1 dominoes? For the answer would be very big, output the answer modulo M (0 < M ≤ 105).
Input
Input consists of several test cases followed by a line containing double 0. Each test case consists of two integers, N and M, respectively.
Output
For each test case, output the answer modules M.
Sample Input
1 10000 3 10000 5 10000 0 0
Sample Output
1 11 95
Source
POJ Monthly--2007.10.06, Dagger
感慨:这道题我郁闷了很久..因为网上没有很好的解题报告...而且标程写得..我看不懂= =所以一直没有理解..现在静静地想一想,其实一切都顺理成章嘛..一次AC,庆祝下,虽然很爽,但郁闷的是,刚blog出了bug,写好的解题报告没了,郁闷= =以后一定记得发之前保存个..
题目大意:给出一个4*n的矩阵,问用1*2的骨牌将它完全覆盖有多少种方案。输出方案数mod m。
分析:传统方法因为n最大有10^9之多,所以会超时。若定义f[i,k]表示前i-1行已经完全覆盖,第i行状态为k的方案数的话,注意到f[i]只和f[i-1]有关,而且关系是线性的,且f的i域范围很大,但k域范围很小(0~15),再者每相邻两行的状态转移方程完全相同,所以可以用矩阵来表示这种转移方法。
构造16*16的矩阵G,使得f[i]*G=f[i+1]。这里,将f[i]看成是一个1*16的矩阵。
只要对矩阵乘法理解透彻,就不难用dfs来构造出矩阵G。然后,只要将第一行可能的状态org乘上Gn次,就能得到f[n+1],f[n+1,0]即答案。
由于矩阵乘法满足结合律,所以可以用快速幂求出G^N,再乘以org。至此,问题解决。
codes: