首先看k的范围,k小于1000,所以可以用杨辉三角来做,也可以用组合数,但杨辉三角好写。
x^n *y^m的系数为杨辉三角的f[k][m] * a^m * b^m;用快速幂来做,不要忘记取模a,b;
#include<iostream>
#include<cstdio>
using namespace std;
const int P=10007;
int a,b,k,n,m,ans;
int f[1100][1100];
int g(int x,int y)
{
int sum=1;x=x%P;
while(y)
{
if(y%2) sum=sum*x%P;
x=x*x%P;
y/=2;
}
return sum%P;
}
int main()
{
scanf("%d%d%d%d%d",&a,&b,&k,&n,&m);
for(int i=1;i<=k+2;i++) f[i][i]=f[i][0]=1;
for(int i=2;i<=k;i++)
for(int j=1;j<=i;j++)
f[i][j]=(f[i-1][j-1]%P+f[i-1][j]%P)%P;
ans=f[k][m]%P*g(a,n)%P*g(b,m)%P;
cout<<ans%P;
return 0;
}