今天的testB调错告诉我一件事，做了算法调整之后，千万注意数组空间是否在一些无关紧要的歩奏越界，不然天知道会出些什么鬼魅的错误。
在DP的时候，只要压维，单调性的理解就显得尤为重要，具体的表现是循环时变量的改变方向一类的细节。
testA复习了一下数位DP，但据说正解是数学推导。naive。
testC周长除以二
今天比较水，难得
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
long long dp[52][2002],c[2008][52];
long long t,k,n,v,s;
int main()
{
freopen("testB.in","r",stdin);
freopen("testB.out","w",stdout);
cin>>t;
memset(dp,0,sizeof(dp));
dp[0][0]=1;
for (int z=1;z<1003;z++)
{
for (int i=50;i>=0;i--)
{
for (int j=z;j<=1003;j++)
{
dp[i+1][j]+=dp[i][j-z];
dp[i+1][j]%=1000000007;
}
}
}
memset(c,0,sizeof(c));
for (int i=1;i<1003;i++)
{
c[i][0]=1;
c[i][i]=1;
}
for (int i=2;i<1003;i++)
{
for (int j=1;j<min(i,50);j++)
{
c[i][j]=c[i-1][j]+c[i-1][j-1];
c[i][j]%=1000000007;
}
}
for (int i=0;i<t;i++)
{
cin>>n>>k;
if (k>50)
{
cout<<0<<endl;
continue;
}
s=0;
for (int j=k;j<=n;j++)
{
v=(dp[k][j]%1000000007)*(c[n-j+k][k]%1000000007)%1000000007;
s=(v+s)%1000000007;
}
for (int z=2;z<=k;z++)
{
s=s*z%1000000007;
}
cout<<s<<endl;
}
return 0;
}