记录编号 |
336900 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[HZOI 2016]寻找saber |
最终得分 |
100 |
用户昵称 |
Go灬Fire |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
1.432 s |
提交时间 |
2016-11-03 19:27:17 |
内存使用 |
152.88 MiB |
显示代码纯文本
#include<stdio.h>
#include<algorithm>
using namespace std;
#define LL long long
#define Begin freopen("saber.in","r",stdin);freopen("saber.out","w",stdout);
#define End fclose(stdin);fclose(stdout);
const int maxn=19999999+10;
const LL mod=19999999;
LL jc[maxn];
inline LL qpow(LL x,LL h){
LL ans=1;
while(h){
if(h&1)ans*=x,ans%=mod;
x*=x;x%=mod;h>>=1;
}
return ans;
}
inline LL C(LL n,LL m){
LL tot=(jc[n]%mod*qpow(jc[m],mod-2)%mod*qpow(jc[n-m],mod-2));
tot%=mod;
return tot;
}
inline LL Lucas(LL n,LL m){
if(m==0)return 1;
return C(n%mod,m%mod)*Lucas(n/mod,m/mod)%mod;
}
void Init(void);
int main(){
Begin;
Init();
getchar();getchar();
//for(;;);
End;
return 0;
}
void Init(void){
LL n,m;scanf("%lld%lld",&n,&m);
jc[0]=1;
for(int i=1;i<=min(n+m+2,mod+2);i++)jc[i]=jc[i-1]*i,jc[i]%=mod;
LL ans=Lucas(n+m,m),now=(n+1)*(n+2)-m*(m-1);now%=mod;
ans*=now;ans%=mod;
ans*=qpow(n+1,mod-2);ans%=mod;
ans*=qpow(n+2,mod-2);ans%=mod;
if(ans<0)ans+=mod;
printf("%lld\n",ans);
}