格子刷油漆（动态规划（递推））

X国的一段古城墙的顶端可以看成 2*N个格子组成的矩形(如图所示)
现需要把这些格子刷上保护漆。

你可以从任意一个格子刷起，刷完一格，可以移动到和它相邻的格子（对角相邻也算数），但不能移动到较远的格子（因为油漆未干不能踩！）

比如：a d b c e f 就是合格的刷漆顺序。
c e f d a b 是另一种合适的方案。

当已知 N 时，求总的方案数。当N较大时，结果会迅速增大，请把结果对 1000000007 (十亿零七) 取模。
输入数据为一个正整数（不大于1000）
输出数据为一个正整数。

例如：
用户输入：
2
程序应该输出：
24

再例如：
用户输入：
3
程序应该输出：
96

再例如：
用户输入：
22
程序应该输出：
359635897
思路：

---

#include<iostream>
#define M 1000000007
using namespace std; 
//从边缘某一点出发到临近点 
long long fb(int n){
if(n==1)return 1;
return fb(n-1)*2%M; 

}
//从边缘某一点出发的所有情况 
long long fa(int n){
    if(n==1)return 1;
    if(n==2)return 6;                   
    return (fb(n)+2*fa(n-1)+4*fa(n-2))%M;
}

//从中间某一点i出发的所有状况。 
long long fk(int i,int n){
    return (fb(i)*fa(n-i)*2%M + fb(n-i+1)*fa(i-1)*2%M)*2%M;
}

long long f(int n){
    if(n==1)return 2;

    //四个边缘角出发的情况 
    long long sum=0;
    sum=fa(n)*4%M;

    for(int i=2;i<n;i++){
        sum=(sum+fk(i,n))%M;
    }
    return sum;
}

int main(){
    for(int i=1;i<=30;i++){
        cout<<i<<": "<<f(i)<<endl; 
    }

}

实际运行，输入数据一大，很慢，所以我们优化为记忆话搜索的方式，避免重复计算。

#include<iostream>
#define M 1000000007
using namespace std; 

long long fa[1000];
long long fb[1000];

//从边缘某一点出发到临近点 
//long long fb(int n){
//if(n==1)return 1;
//return fb(n-1)*2%M; 
//  
//}
////从边缘某一点出发的所有情况 
//long long fa(int n){
//  if(n==1)return 1;
//  if(n==2)return 6;                   
//  return (fb(n)+2*fa(n-1)+4*fa(n-2))%M;
//}
//从中间某一点i出发的所有状况。 
long long fk(int i,int n){
//  return (fb(i)*fa(n-i)*2%M + fb(n-i+1)*fa(i-1)*2%M)*2%M;
return (fb[i]*fa[n-i]*2%M+fb[n-i+1]*fa[i-1]*2%M)*2%M; 
}

long long f(int n){
    if(n==1)return 2;

    //四个边缘角出发的情况 
    long long sum=0;
    sum=fa[n]*4%M;

    for(int i=2;i<n;i++){
        sum=(sum+fk(i,n))%M;
    }
    return sum;
}

int main(){
    int i;
    cin>>i;
    fb[1]=1;
    for(int j=2;j<1000;j++){
        fb[j]=fb[j-1]*2%M;
    }
    fa[1]=1;
    fa[2]=6;
    for(int j=3;j<1000;j++){
        fa[j]=(2*fa[j-1]+4*fa[j-2]+fb[j])%M;
    }

    cout<<f(i)<<endl; 
return 0;
}