洛谷 P1313 计算系数（杨辉三角+快速幂）

传送门

---

我们知道(by+ax)^k的展开式中，第x^n*y^m的系数为杨辉三角形中的第(k,m)项乘以a^n*b^m，因为数据太大，后者需要用快速幂解决。关于快速幂，可以看前文

Code:

#include<cstdio>
#include<cstdlib>
#define MOD 10007

int f[2010][2010];
int a,b,k,n,m;

int ksm(int x,int y)
{
    int ans=1;
    x%=10007;
    while(y>=1)
    {
        if(y&1) ans=ans*x%MOD;
        x=x*x%MOD;
        y/=2;
    }
    return ans;
}

int main()
{
    scanf("%d %d %d %d %d",&a,&b,&k,&n,&m);
    for(int i=1;i<=k;i++) f[i][0]=f[i][i]=1;
    for(int i=1;i<=k;i++)
        for(int j=1;j<i;j++)
            f[i][j]=(f[i-1][j-1]+f[i-1][j])%MOD;
    printf("%d",(f[k][m]*ksm(a,n)%MOD)*ksm(b,m)%MOD);
}