题245.数论-洛谷P1313 [NOIP2011 提高组] 计算系数

---

题245.数论-洛谷P1313 [NOIP2011 提高组] 计算系数

---

一、题目

二、题解

首先推导公式我们可以知道我们要求的系数其实就是a的n次方乘以b的m次方再乘上m,k组合数。对于a的n次方和b的m次方，我们可以使用快速幂直接得出。

之后，看k,m组合数，它可作如下推导：

则我们把除以(k-m)!分解，其实就是每次去乘i逆元(即倒数)，又因为取模是素数，所以可以费马小定理a^-1%mod=a^mod-2，得到结果。
则代码如下：

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

const int mod=10007;

ll fastPow(ll base,ll pow)
{
    ll res=1;
    base=base%mod;
    while(pow>0)
    {
        if(pow&1)
        {
            res=res*base%mod;
        }
        pow>>=1;
        base=(base*base)%mod;
    }
    return res;
}

int main()
{
    int a,b,k,n,m;
    cin>>a>>b>>k>>n>>m;
    int res;
    res=fastPow(a,n)*fastPow(b,m)%mod;
    for(int i=1,j=0;i<=n,j<=m-1;i++,j++)
    {
        res=res*(k-j)%mod;//计算k(k-1)...(k-m+1)=k(k-1)...(k-(m-1))
        res=res*fastPow(i,mod-2)%mod;//计算(1/i)%mod，用费马小定理
    }
    cout<<res<<endl;
}

---