const int M=1e3;
const int N=2e3;
int c[N][N];
void init(){
int i,j;
c[1][1]=1;
for(i=0;i<=M;i++) c[i][0]=1;
for(i=2;i<=M;i++)
for(j=1;j<=i;j++)
c[i][j]=(c[i-1][j]+c[i-1][j-1]);
}
n较小
long long comb(long long n, int k) { auto res = n; for (int i = 2; i <= k; i++) res = res * --n / i; return res %mod; }