（算法竞赛进阶指南）计算系数（逆元+组合数）

题目链接：211. 计算系数 - AcWing题库

分析：这道题目是要我们求组合数，显然答案等于a^n*b^m*C(k,n)，前两个指数都比较容易求出，主要是求后面一个组合数会涉及到除法，因为需要取模，所以当过程中对mod取余后不一定能整除，所以我们在对mod取余的同时除以一个数，只能转化为乘以该数关于mod的逆元，逆元其实也比较好求，由费马定理得当a和p互质时有a^(p-1)%p=1，所以a关于p的逆元就是a^(p-2)，由于该题目中p为10007是一个质数，所以a一定与其互质，所以求逆元就可以直接用快速幂求

下面是代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
using namespace std;
typedef long long ll;
const int N=1003,mod=10007;
ll fact[N];
ll qpow(ll a,ll b)
{
	ll ans=1;
	while(b)
	{
		if(b&1) ans=ans*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return ans;
}
int main()
{
	fact[1]=1;
	for(int i=2;i<=1000;i++)
		fact[i]=fact[i-1]*i%mod;
	ll a,b,k,n,m;
	cin>>a>>b>>k>>n>>m;
	ll ans=qpow(a,n)*fact[k]%mod*qpow(fact[n],mod-2)%mod*qpow(fact[k-n],mod-2)%mod*qpow(b,m)%mod;
	printf("%lld",ans);
	return 0;
}