19.11.30联赛 Problem E. （组合数打表）

Problem E. 小埋与集合

组合数打表

1.求组合数

ll C(ll n,ll m) 
{
   ll ret=1;
   for(ll i=1; i<=m; ++i)
      ret=ret*(n-i+1)/i;
   return ret;
}

2.快速幂求

typedef long long int ll;
const int mod=1e9+7;
ll qpow(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1)ans=ans*a%mod;
        b>>=1;
        a=a*a%mod;
    }
    return ans;
}
ll jc(ll n){
    ll ans=1;
    for(ll i=1;i<=n;i++)ans=ans*i%mod;
    return ans;
}
ll C(ll n,ll m){
    ll ans=1;
    ans=ans*jc(n)%mod;
    ans=ans*qpow(jc(m),mod-2)%mod;
    ans=ans*qpow(jc(n-m),mod-2)%mod;
    return ans;
}

3.适用于N<=3000

//c[i][j]表示从i个中选j个的选法。
long long C[N][N];

void get_C(int maxn)
{
    C[0][0] = 1;
    for(int i=1;i<=maxn;i++)
    {
        C[i][0] = 1;
        for(int j=1;j<=i;j++)
            C[i][j] = C[i-1][j]+C[i-1][j-1];
        //C[i][j] = (C[i-1][j]+C[i-1][j-1])%MOD;
    }
}

4.线性打表

const ll maxn=1e5+50;
ll a[maxn],b[maxn];
//a[i]是i的阶乘，b[i]是阶乘的除法逆元，两者用于求组合数
ll qpow(ll n,ll k,ll mod){
    ll res=1;
    n=n%mod;
    while(k){
        if(k&1)res=res*n%mod;
        n=n*n%mod;
        k>>=1;
    }
    return res;
}
void init(){
    a[0]=b[0]=1;
	for(int i=1;i<maxn;++i){
        a[i]=(a[i-1]*i)%mod;
	}
	b[maxn-1]=qpow(a[maxn-1],mod-2,mod);		
	//maxn-1必须小于mod
	for(int i=maxn-1;i>0;i--)
        b[i-1]=b[i]*i%mod;
}
ll C(ll n,ll m){
    if(n<m)return 0;
    return  a[n]*b[m]%mod*b[n-m]%mod;
}
 
//Lucas(n,m)      //C(N,M)
ll Lucas(ll n, ll m)
{	
//卢卡斯定理，求大组合数，且mod必须是质数，
//否则需要先对mod拆分然后用中国剩余定理取合并
    if(m ==0) return 1;
    else return (C(n%mod, m%mod)*Lucas(n/mod, m/mod))%mod;
}

题目：

Problem E. 小埋与集合
Input file: standard input
Output file: standard output
Time limit: 1 seconds
Memory limit: 512 megabytes
小埋是一个大学霸同时体育也是达人，在外面她是一位无可挑剔的美少女。聪明的小埋也喜欢数学问题，这天她在思考一道集合问题：
给定一个包含 n 个正整数元素的集合，这个集合满足如下要求
𝑎1 < 𝑎2 < 𝑎3 < ⋯ < 𝑎𝑛，
并且保证给定一个正整数 m
1 ≤ 𝑎1, 𝑎2, 𝑎3, … , 𝑎𝑛 ≤ 𝑛 + 𝑚
现在她想知道这 n 个元素的集合有多少种组成方式？答案可
能过于巨大我们需要将答案模 998244353。
Input
输入包含多组测试数据。
第一行一个 T (1 ≤ 𝑇 ≤ 106
) 代表测试数据组数。
第二行两个正整数 n，m 代表 n 个元素的集合以及上面提到的 m
值(1 ≤ 𝑛 ≤ 2019 0 ≤ 𝑚 ≤ 2019 )。
Output
对于每组样例输出一行一个答案代表结果模 998244353。
Sample input and output
1
3 0
Sample Input Sample Output
1
提示
对于第一个测试数据我们仅仅能组成{1，2，3}所以答案为 1。

AC代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
using namespace std;
long long c[4050][4050];
int main()
{
    c[0][0]=1;
    for(int i=1;i<4050;i++)
    {
        c[i][0]=1;
        for(int j=1;j<=i;j++)
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%998244353;
    }///组合数打表
    int t,n,m;
    cin>>t;
    while(t--)
    {
        scanf("%d%d",&n,&m);
        printf("%lld\n",c[n+m][n]);
    }
    return 0;
}