hdu4301状态较多的dp

给出2*n的巧克力，输入n，k，输出这2*n的巧克力分成k块有多少种分法，n《=1000，(想到n^2刚好。。。）

dp[i][j][0]代表前1到i列分成j块并且最后一列是分开

dp[i][j][1]代表前1到i列分成j快并且最后一列是断开

#include <iostream>
#include <cstdio>
using namespace std;
typedef long long ll;
const int maxn=1e3;
const ll mod=1e8+7;
ll dp[maxn+100][2*maxn+100][2];
void init();
int main()
{
    int t,n,k;
    init();
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&k);
        printf("%lld\n",(dp[n][k][0]+dp[n][k][1])%mod);
    }
    return 0;
}
void init()
{
    for(int i=1;i<=1000;i++)
        dp[i][1][1]=1;
    dp[1][2][0]=1;
    for(int i=2;i<=1000;i++)
    {
        for(int j=2;j<=i*2;j++)
        {
            dp[i][j][0]=(dp[i][j][0]+dp[i-1][j-2][0]+dp[i-1][j-2][1])%mod;
            dp[i][j][1]=(dp[i][j][1]+dp[i-1][j-1][0]+dp[i-1][j-1][1])%mod;
            dp[i][j][0]=(dp[i][j][0]+2*(dp[i-1][j-1][0]+dp[i-1][j-1][1]))%mod;
            dp[i][j][1]=(dp[i][j][1]+2*dp[i-1][j][0])%mod;
            dp[i][j][1]=(dp[i][j][1]+dp[i-1][j][1])%mod;
            dp[i][j][0]=(dp[i][j][0]+dp[i-1][j][0])%mod;
        }
    }
}