我们知道(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);
}