记录编号 |
174891 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[NOIP 2011]计算系数 |
最终得分 |
100 |
用户昵称 |
<蒟蒻>我要喝豆奶 |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
0.003 s |
提交时间 |
2015-08-03 17:48:25 |
内存使用 |
0.32 MiB |
显示代码纯文本
#include<cstdio>
#include<iostream>
#define p 10007
using namespace std;
typedef long long LL;
inline int in(){
char c=getchar();
int x=0;
while(c<'0'||c>'9')c=getchar();
for(;c>='0'&&c<='9';c=getchar())x=x*10+c-'0';
return x;
}
LL n,m;
LL a,b;
LL k;
LL frm(LL a,LL b){
LL ans=1;
while(b){
if(b&1){
ans=ans*a%p;
}
b>>=1,a=(a*a)%p;
}
return ans;
}
LL C(LL n,LL m){
if(m>n)
return 0;
LL ans=1;
for(int i=1; i<=m; i++) {
LL a=(n+i-m)%p;
LL b=i%p;
ans=ans*(a*frm(b, p-2)%p)%p;
}
return ans;
}
LL Lucas(LL n,LL m){
if(m==0)
return 1;
return C(n%p,m%p)%p*(Lucas(n/p,m/p)%p);
}
int main(){
freopen("factor.in","r",stdin);
freopen("factor.out","w",stdout);
a=in(),b=in(),k=in(),n=in(),m=in(),printf("%lld",((Lucas(k,n)%p)*(frm(a,n)%p)*(frm(b,m)%p))%p);
return 0;
}