计算系数

先导

根据高中的数学知识二项式定理，我们可以得到

((x+y)^n=C_{n}^0x^ny^0+C_{n-1}^1x^{n-1}y1+...+C_{n}^nx^0y^n)

我们可以这样理解，比如这样一个式子，((x+y)^5)

如果我们要求(x^3y^2)的系数，那么我们想象成有(()()()()())个括号连续相乘，为了得到(x^3y^2)

我们需要从(5)个括号中选出(3)个来作为(x)相乘的项，即(C_{5}^3)((x)^3)

另外还剩下两个括号，我们选出(2)个括号作为(y)相乘的项，即(C_{2}^2(y)^2)

最后把这两个分散的项相乘得到(C_{5}^3(x)^3C_{2}^2(y)^2)

推广

那么进行推广我们可以得到

((x+y+z)^n)如果要求(x^ay^bz^c)的系数

还是想象成有(n)个括号连续相乘

从(n)个括号中选出(a)个凑成(x^a)，即(C_{n}^a(x)^a)。

从剩下的(n-a)个括号选出(b)个凑出(y^b)，即(C_{n-a}^b(y)^b)。

最后从剩下的(n-a-b)个括号中选出(c)个凑出(z^c)，即(C_{n-a-b}^c(z)^c)

最后的最后把这些分散的项相乘就是最终(x^ay^bz^c)的系数

（注意，我这里之所以要把(x,y,z)外面带上括号，就是因为某些情况下，(x)前面的系数可能不为1，这个时候还要把(x)前面的系数也算上）

比如要求(x^3y^2)

((2x+3y)^5=C_{5}^3(2x)^3C_{2}^2(3y)^2)

有了上面的知识储备，下面的这道题一定就非常简单了

真题

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

输入格式

共一行，包含 5 个整数，分别为 (a，b，k，n，m)，每两个整数之间用一个空格隔开。

输出格式

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

数据范围

(0≤n,m≤k≤1000)
(n+m=k)
(0≤a,b≤10^6)

输入样例：

1 1 3 1 2 

输出样例：

3

题解

从连续的(k)个括号中选出(n)个作为凑成(x^n)

(C_{k}^n(ax)^n)

从剩下的(k-n)个括号中选出(m)个凑成(y^m)

(C_{k-n}^{m}(by)^m)

最后答案(C_{k}^n(ax)^n)(C_{k-n}^{m}(by)^m)

所以我们的代码只需要处理

(a^n)，(b^m)，还有组合数的计算就行

代码

#include <cmath>
#include <iostream>

using namespace std;
const int mod = 10007;
const int N = 1010;
int a, b, k, n, m;
int c[N][N];
void init() {
  for (int i = 0; i <= 1000; i++)
    for (int j = 0; j <= i; j++) {
      if (j == 0)
        c[i][j] = 1;
      else
        c[i][j] = (c[i - 1][j - 1] + c[i - 1][j]) % mod;
    }
}

int qmi(int a, int b) {
  int ans = 1;
  while (b) {
    if (b & 1) ans = (long long)ans * a % mod;
    a = (long long)a * a % mod;
    b >>= 1;
  }
  return ans % mod;
}
int main() {
  init();  //预处理出组合数
  cin >> a >> b >> k >> n >> m;
  cout << (long long)qmi(a, n) * qmi(b, m) * c[k][n] % mod;
  return 0;
}