1137 计算系数（组合数）

给定一个多项式(ax + by)^k，请求出多项式展开后x^n y^m项的系数。

输入描述 Input Description
共一行，包含 5 个整数，分别为a，b，k，n，m，每两个整数之间用一个空格隔开。

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

样例输入 Sample Input
1 1 3 1 2

样例输出 Sample Output
3
组合数的预处理啊，应用~

#include<iostream>
#include<cstdio>
#include<vector>
#include<algorithm>
#include<cstring>
#include<cmath>
#define mod 10007
using namespace std;
int _p(int a,int b)
{
    int ans=1;
    while(b)
    {
        if(b&1)ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}
int fac[1005];
int inv[1005];
void init()
{
    fac[0]=1;
    for(int i=1;i<=1000;i++)
        fac[i]=fac[i-1]*i%mod;
    inv[1000]=_p(fac[1000],mod-2);
    for(int i=999;i>=0;i--)
        inv[i]=inv[i+1]*(i+1)%mod;
}
int get(int n,int m)
{
    return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int main()
{
    int a,b,k,n,m;
    init();
    cin>>a>>b>>k>>n>>m;
    cout<<get(k,n)*_p(a,n)%mod*_p(b,m)%mod;
    return 0;
}