拖了很久的组合数取模求法。。。先更一个,逆元+快速幂。
何为逆元?又有人称之为数论倒数。那a的数论倒数当然就是a^-1啦,然而数论中都是整数,所以严格地定义为 a关于p的逆元为满足
ax≡1(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的逆元为
ap−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;
}