杨辉三角+快速幂

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll qpow(ll a,ll b,ll p){ //求a^b%p
    ll ans=1;
    while(b){
        if(b&1)ans=ans*a%p;
        a=a*a%p;
        b>>=1;
    }
    return ans;
}
int a,b,k,n,m;
ll C[1002][1002];
int main(){
    cin>>a>>b>>k>>n>>m;
    for(int i=0;i<1002;++i)
        for(int j=0;j<=i;++j)
            C[i][j]=1;
    for(int i=2;i<=k;++i)
        for(int j=1;j<i;++j)
            C[i][j]=(C[i-1][j]+C[i-1][j-1])%10007;
    cout<<(C[k][n]*qpow(a,n,10007)*qpow(b,m,10007))%10007;
    return 0;
}