计算系数

给定一个多项式(ax+by)^k，请求出多项式展开后xⁿy^m项的系数。
输入描述:
共一行，包含5个整数，分别为a，b，k，n，m，每两个整数之间用一个空格隔开。

输出描述:
输出共1行，包含一个整数，表示所求的系数，这个系数可能很大，输出对10007取模后的结果。
示例1
输入
1 1 3 1 2
输出
3

备注:
对于30%的数据，有0≤k≤10；
对于50%的数据，有a=1，b=1；
对于100%的数据，有0≤k≤1,000，0≤n,m≤k，且n+m=k，0≤a，b≤1,000,000。

#include<bits/stdc++.h>
using namespace std;
const int MOD=10007;
long long qpow(long long a,long long b)
{
    long long res=1;
    while(b)
    {
        if(b&1)res*=a;
        res%=MOD;
        b/=2;
        a=(a*a)%MOD;
    }
    return res%MOD;
}
long long C(long long n,long long m)
{
    long long c[1005];
    c[0]=1;
    for(int i=1;i<=n;i++)
    {
        c[i] = ((c[i - 1] * (m - i + 1))%MOD) * qpow(i, MOD - 2);
        c[i] %= MOD;
    }
    return c[n]%MOD;
}
int main()
{
    int a,b,k,n,m;
    cin>>a>>b>>k>>n>>m;
    cout<<((C(m,k)*qpow(a,n))%MOD*qpow(b,m))%MOD<<endl;
    return 0;
}