重蹈覆辙

题目描述

用L图形（大小为3，也就是去掉一个角的2x2的正方形）和1x2的矩形，覆盖2xn的矩形，问有多少种方案。 覆盖要求不重不漏，整体翻转和旋转均算作不同的方案。 用于覆盖的图形可以旋转，比如可以把L旋转为Г，把1x2的矩形旋转成为2x1的矩形等。 输出方案数模10007的结果。 对于100%的数据，1 <= n <= 10100000 对于30%的数据，n <= 105 对于60%的数据，n <= 1018 对于80%的数据，n <= 101000
输入描述:
输入一行一个整数n。
输出描述:
输出方案数模10007的结果。

示例1

输入

5

输出

24

示例2

输入

9

输出

569

示例3

输入

12

输出

6105

示例4

输入

13

输出

3458

示例5

输入

1

输出

1

示例6

输入

2

输出

2

示例7

输入

3

输出

5

示例8

输入

4

输出

11

思路

用 L 形和条形来覆盖 2 × ? 的矩形，求方案数

从左到右拼，每次都填充横坐标最小的空格

那么每时每刻，我们当前的状态一定只有两种
——最后一列完全填满的
——最后一列空一格的

记录总共填到第 ? 列，且最后一列空 ? 格的方案数为 ?（?,?） ，那么我们有 ?( ?,0 )= ? (? − 1,0 )+ ? (? − 1,1) + ?( ? − 2,0 )
? (?,1 )= 2? (? − 2,0) + ? (? − 1,1)
或者也可以找规律发现 ? (?,0 )= 2? (? − 1,0) + ?( ? − 3,0)

但是 ? 特别的大

解决方法一：十进制（矩阵）快速幂，
$A^n=({a\lfloor \frac{n}{10}\rfloor )}^{10}{A}^{nmod10}$

解决方法二：找循环节，发现循环节是 10006……于是直接暴力就好了

我用的解法二

代码

#include<bits/stdc++.h>
#define endline putchar('\n')
using namespace std;
template <typename T>
void read(T &x)
{
    x=0;
    char c=getchar();
    int sgn=1;
    while(c<'0'||c>'9')
 {
  if(c=='-')
   sgn=-1;
  c=getchar();
   }
    while(c>='0'&&c<='9')
  x=x*10+c-'0',c=getchar();
    x*=sgn;
}
template <typename T>
void out(T x)
{
    if(x<0)
 {
  putchar('-');
  x=-x;
 }
    if(x>=10)
  out(x/10);
    putchar(x%10+'0');
}
long long a[100005];
const int p=10007;
int main()
{
    a[0]=a[1]=1;
    a[2]=2;
    for(int i=3;i<=100000;i++)
        a[i]=(2*a[i-1]+a[i-3])%p;
    string str;
    cin>>str;
    int n=str.size();
    long long x=0;
    for(int i=0;i<n;i++)
        x=(x*10+str[i]-'0')%10006;
    cout<<a[x]<<endl;
    return 0;
}

来源：nkw