组合数打表

当n,m都很小的时候可以利用杨辉三角直接求(<=3000)
C(n,m)=C(n-1,m)+C(n-1,m-1);

const int CN = 55;
long long c[CN][CN] = {};
void cinit()
{
    for (int i = 0; i < CN; i++)
    {
        c[i][0] = c[i][i] = 1;
        for (int j = 1; j < i; j++)
        {
            c[i][j] = c[i - 1][j] + c[i - 1][j - 1];
        }
    }
}

当n固定时对C(n,i)打表,n可以很大

typedef long long ll;
const int M = 1e9 + 7;
const ll MAXN = 1e6 + 50;
ll inv[MAXN + 200];
ll C[MAXN];
void getinv()
{
    inv[1] = 1;
    for (int i = 2; i <= MAXN; i++)
    {
        inv[i] = inv[M % i] * (M - M / i) % M;
    }
}
void getfac(ll n)
{
    C[0] = 1;
    for (int i = 1; i <= n; i++)
    {
        C[i] = C[i - 1] * (n - i + 1) % M * inv[i] % M;
    }
}
int main()
{
    getinv();
    ll n, m;
    cin >> n;
    getfac(n);
    while (cin >> m)
    {
        cout << C[m] << endl;
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值