luogu1313【2011提高】计算系数（求组合数取模）

拖了很久的组合数取模求法。。。先更一个，逆元+快速幂。
何为逆元？又有人称之为数论倒数。那a的数论倒数当然就是a^-1啦，然而数论中都是整数，所以严格地定义为 a关于p的逆元为满足$ax\equiv1(mod\ p)$的最小正整数x。那逆元有什么用呢？四则运算的取余运算基本是相同的，但是除法是不满足的，即a/b%p!=(a%p)/(b%p)%p，那怎么办呢。。。a/b=a*b^-1！！！也就是说，a/b%p=a * b^-1 %p =(a%p) *(b^-1%p) %p。注意这里的b^-1可是b的数论倒数哟。那a关于p的逆元怎么求呢？如果(a,p)=1并且p为素数，则根据费马小定理，有a关于p的逆元为$a^{p-2}$。具体证明请见百度。这样我们就可以用逆元+快速幂来计算组合数取模了！

#include <bits/stdc++.h>
using namespace std;
#define N 1010
#define inf 0x3f3f3f3f
#define ll long long
#define mod 10007
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int a,b,n,m,kk,ans=1;
inline int ksm(int x,int k){
    int res=1,base=x;
    for(;k;k>>=1,base=base*base%mod)
        if(k&1) res=res*base%mod;
    return res;
}
int main(){
//  freopen("a.in","r",stdin);
    a=read()%mod;b=read()%mod;kk=read();n=read();m=read();
    for(int i=1;i<=n;++i) ans=ans*i%mod;
    ans=ksm(ans,mod-2);
    for(int i=kk-n+1;i<=kk;++i) ans=ans*i%mod;
    ans=(ll)ans*ksm(a,n)%mod*ksm(b,m)%mod;
    printf("%d\n",ans);
    return 0;
}