[NOIP模拟][动态规划]Table

题目描述：
题目背景： SOURCE：NOIP2016-RZZ-2 T2
给定一个 n×m 的矩阵，行列均从 1 开始标号。一个矩阵被认为是稳定的，当且仅当对于任意的 2≤i≤n，第 i 行的数的和不小于第 i−1 行的数的和，且最后一行的数的和小于等于 m ，并且要求矩阵中所有的元素都是非负的。求所有 n×m 的稳定矩阵的方案数，答案对 $10^9$ 取模。
输入格式:
第一行一个整数 T ，表示数据组数。
每组数据一行两个整数 n，m 。
输出格式:
输出 T 行，每行一个整数，表示方案数。
样例输入：
3
1 1
2 2
2 3
样例输出：
2
25
273
数据规模与约定：
对于 30% 的数据，n,m≤3。
对于 60% 的数据，n,m≤50。
对于 100% 的数据，1≤n,m≤2000；1≤T≤10。
题目分析：
动态规划。不过我的方法和标算不太一样。
我的方法：先预处理出一个二维数组p[i][j],代表i个格子时，和为j时的方案数。转移是p[i][j]=p[i-1][j]+p[i][j-1]。可以这样想对于p[i][j]，第一个格子可以放从j到0,对于剩下的i-1格子，方案数不就是p[i-1][0]+p[i-1][1]+·······+p[i-1][j-1]+p[i-1][j];而同时除去p[i-1][j]的前面部分，就是同理的p[i][j-1];于是得到了p[i][j]的转移方程。然后进行答案的计算，二维数组f[i][j]代表i行时和为j的方案数（列数固定为m）。于是有转移方程f[i][j]=f[i][j-1]+f[i-1][j]*p[m][j]。这是显然的，因为总的方案数f[i][j]是包括f[i][j-1]的，当和增大时，增加的方案数为m个格子和为j时的方案数（即p[m][j]）乘上前面i -1行的总的方案数（f[i-1][j]）。于是得到正确答案。
标算的方法就自己看代码了······懒·····下面分别附了我的方法的代码和标算。
附代码：
我的代码：

#include<iostream>
#include<cstring>
#include<string>
#include<cstdlib>
#include<cstdio>
#include<ctime>
#include<queue>
#include<iomanip>
#include<cmath>
#include<cctype>
#include<set>
#include<map>
#include<algorithm>
using namespace std;

const int mod=1e9;
int t,n,m;
long long p[2010][2010],f[2010][2010],sum[2010][2010];

int readint()
{
    char ch;int f=1,i=0;
    for(ch=getchar();(ch<'0'||ch>'9')&&ch!='-';ch=getchar());
    if(ch=='-')
    {
        ch=getchar();
        f=-1;
    }
    for(;ch>='0'&&ch<='9';ch=getchar())
        i=(i<<3)+(i<<1)+ch-'0';
    return i*f;
}

void pre()
{
    for(int i=1;i<=2000;++i)
    {
        p[i][0]=1;//无论多少个格子，和为0时，永远只有1种放法。
        p[1][i]=1;//一个格子，也永远只有1种放法
    }
    for(int i=1;i<=2000;++i)
        for(int j=1;j<=2000;++j)
            p[i][j]=(p[i][j-1]+p[i-1][j])%mod;
    return ;
}

int main()
{
    //freopen("table.in","r",stdin);
    //freopen("table.out","w",stdout);

    pre();
    t=readint();
    while(t--)
    {
        n=readint();m=readint();
        memset(f,0,sizeof(f));
        for(int i=0;i<=m;i++)
            f[1][i]=(f[1][i-1]+p[m][i])%mod;//第一行就是累加，其实是可以把f[0][j]看作1
        for(int i=2;i<=n;i++)
            for(int j=0;j<=m;j++)
                f[i][j]=(f[i][j-1]+f[i-1][j]*p[m][j])%mod;
        printf("%I64d\n",f[n][m]);
    }

    return 0;
}

标算：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<string>
#include<cstring>
#include<queue>
using namespace std;

const int nod = 1000000000;

int n,m,t;
int c[4005][4005],s[2005][2005];

void pre()
{
    for(int i=0;i<=4000;i++) c[i][0] = c[i][i] = 1;
    for(int i=2;i<=4000;i++)
      for(int j=1;j<i;j++)
        c[i][j] = (c[i-1][j] + c[i-1][j-1])%nod;
}

void solve()
{
    memset(s,0,sizeof(s));
    cin >> n >> m;
    for(int i=0;i<=m;i++)   s[0][i] = 1;
    for(int i=1;i<=n;i++)
      for(int j=0;j<=m;j++)
      {
        long long now = s[i-1][j];
        now = now*c[j+m-1][m-1]%nod;
        if(j)   s[i][j] = s[i][j-1] + now;
        else    s[i][j] = now;
        s[i][j] %= nod;
      }
      cout << s[n][m] << endl;
}

int main()
{
    pre();
    cin >> t;
    while(t--)  solve();
}