刷题记录（NC16596 [NOIP2011]计算系数（排列组合））

最新推荐文章于 2024-07-25 11:38:39 发布

rootting

最新推荐文章于 2024-07-25 11:38:39 发布

阅读量51

点赞数

分类专栏：刷题记录文章标签：算法

本文链接：https://blog.csdn.net/m0_60531106/article/details/128335267

版权

刷题记录专栏收录该内容

103 篇文章 0 订阅

订阅专栏

NC16596 [NOIP2011]计算系数

题目链接

关键点：

1、对于 $(a*x + b*y)^{k}$ 里面 $a^{n} * b^{m}$ 的系数为（因为不会打排列组合的符号，就用文字描述）k的n的组合 * $a^{n}$ * $b^{m}$

2、对于求一个数的次方，可以采用快速幂

ll kuai(ll x, ll y)
{
    ll ans = 1;
    ll base = x;
    while (y)
    {
        if (y&1) ans = (ans%MOD * base%MOD)%MOD;
        base = (base%MOD*base%MOD)%MOD;
        y>>=1;
    }
    return ans;
}

3、对于组合数，可以利用杨辉三角快速求出

    for (int i=0; i<=k; i++)
        f[i][0] = 1;
    for (int i=1; i<=k; i++)
    {
        for (int j=1; j<=i; j++)
            f[i][j] = (f[i-1][j] + f[i-1][j-1])%MOD;
    }

完整代码：

# include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MOD = 10007;
ll a, b, k, n, m;
ll f[2000][2000];
ll kuai(ll x, ll y)
{
    ll ans = 1;
    ll base = x;
    while (y)
    {
        if (y&1) ans = (ans%MOD * base%MOD)%MOD;
        base = (base%MOD*base%MOD)%MOD;
        y>>=1;
    }
    return ans;
}
int main()
{
    scanf("%lld%lld%lld%lld%lld", &a, &b, &k, &n, &m);
    for (int i=0; i<=k; i++)
        f[i][0] = 1;
    for (int i=1; i<=k; i++)
    {
        for (int j=1; j<=i; j++)
            f[i][j] = (f[i-1][j] + f[i-1][j-1])%MOD;
    }
    ll ans = f[k][n]%MOD * kuai(a, n)%MOD * kuai(b, m)%MOD;
    cout<<ans%MOD<<endl;
    
    
    
    return 0;
}