填格子3*N的方框使用2*1的矩形进行填充

考虑每个位置的前一个状态 可以发现有

我们分别给他们编号 假设 现在填充到了i+1行，我们可以发现从i行可以通过填充转到i+1行的状态

第i行第j列表示 可以 从上一个转态 j 可以到达这个状态的j ， 然后又了这个转态转化在使用矩阵

{00000001}为初始状态

#include <iostream>
#include <algorithm>
#include <string.h>
#include <vector>
#include <cstdio>
using namespace std;
typedef long long LL;
const LL mod=;
struct Matir{
   LL va[][];
   Matir(){
       memset(va,,sizeof(va));
   }
   Matir operator *(const Matir rhs){
      Matir ans;
         for(int i=; i<; i++)
         for(int j=; j<; j++){
                 for(int k=; k<; k++){
                   ans.va[i][j]=(ans.va[i][j]+(va[i][k]*rhs.va[k][j])%mod )%mod;
                 }
         }
         return ans;
   }
}model;
Matir solve(int N){
     Matir ans;
     for(int i=; i<; i++)ans.va[i][i]=;
     Matir A=model;
     while(N){
        if(N&)ans=ans*A;
        N>>=;
        A=A*A;
     }
     return ans;
}
int main()
{
    for(int i=; i<; i++)
        model.va[i][-i-]=;
        model.va[][]=; model.va[][]=; model.va[][]=;model.va[][]=;
      int N;
      while(scanf("%d",&N)==){
          if(N&){
             printf("0\n");continue;
          }
          Matir ans=solve(N);
          cout<<ans.va[][]<<endl;
      }
    return ;
}