给定一个多项式(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;
}