去年多校 我买了块表 之 hdu 4301 简单DP

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int INF=1<<29;
const int maxn=1002;
const int mod=100000007;
long long dp[maxn][maxn*2];///表示第i列被分成j块，且第i列未被分开的情况数
int f[maxn][maxn*2];///表示第i列被分成j块，且第i列被分开的情况数
/*
边界条件 f[1][2]=dp[1][1]=1;
那么状态转移方程为：
  dp[i][j]=dp[i-1][j]+dp[i-1][j-1]+f[i-1][j-1]+f[i-1][j]*2;
  f[i][j]=2*(dp[i-1][j-1]+f[i-1][j-1])+f[i-1][j];
  f[i][j]+=dp[i-1][j-2]+f[i-1][j-2];

*/
int t,n,k;
int main()
{
  //  freopen("//media/学习/ACM/input.txt","r",stdin);
    cin>>t;
    int i,j;
    f[1][2]=dp[1][1]=1;
    for(i=2;i<=1000;i++)
    {
        for(j=1;j<=2*i;j++)
        {
            dp[i][j]=(dp[i-1][j]+dp[i-1][j-1]+f[i-1][j-1]+f[i-1][j]*2)%mod;
            f[i][j]=(2*(dp[i-1][j-1]+f[i-1][j-1])+f[i-1][j])%mod;
            if(j>=2)
            f[i][j]=(f[i][j]+dp[i-1][j-2]+f[i-1][j-2])%mod;
        }
    }
    while(t--)
    {
        cin>>n>>k;
        cout<<(dp[n][k]+f[n][k])%mod<<endl;
    }
    return 0;
}