AcWing 885. 求组合数 I(dp)

Acwing 885
给定n组询问，每组询问给定两个整数a，b，请你输出C^b_a mod (10⁹+7)的值。

输入格式
第一行包含整数n。

接下来n行，每行包含一组a和b。

输出格式
共n行，每行输出一个询问的解。

数据范围
1≤n≤10000,
1≤b≤a≤2000
输入样例：

3
3 1
5 3
2 2

输出样例：

3
10
1

这道题我刚开始想到用阶乘去解决。
C^b_a=a!/(b!*!(a-b))
然后小数是可以算出来的，但是mod的问题没解决好，时间复杂度也比较高，这个类似于筛法。
用了一个简单dp
我们高中数学学过C^b_a=C^b_a-1+C^b-1_a-1
这个刚好可以用来去dp
我们怎么去dp呢
代码如下

#include<iostream>

using namespace std;

const int N=2020;
const int mod=1e9+7;

int c[N][N];

int main(void)
{
    for(int i=0;i<N;i++)
    for(int j=0;j<=i;j++)
    {
        if(!j)
        c[i][j]=1;
        else
        c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
    }
    int n;
    scanf("%d",&n);
    while(n--)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        printf("%d\n",c[a][b]);
    }
}