关注
#include<algorithm>
#include<cstdio>
#include<queue>
#include<cmath>
#include<cstring>
#include<iostream>
#include<ctime>
#include<cstdlib>
#define ll long long
using namespace std;
const int mod=1e9+7;
const int maxn=10000005;
ll n,m,k;
ll ans;
ll q_pow(ll n,ll mi){
ll res=1,temp=n%mod;
while(mi){
if(mi&1) res=res*temp%mod;
temp=temp*temp%mod;
mi>>=1;
}
return res;
}
ll cal(ll n,ll m){ // Cm(n,m)=(n!/(n-m)!) * (m!)^(mod-2)) mod mod
if(m>n) return 0; // important
ll res=1;
for(int i=1;i<=m;i++){
ll t1=(n-m+i)%mod,t2=i%mod;
res=res*(t1*q_pow(t2,mod-2)%mod)%mod;
}
return res;
}
/*Lucas(n,m,mod)=Cm(n%mod,m%mod)* Lucas(n/mod,m/mod,mod)
Lucas(x,0,mod)=1;*/
ll lucas(ll t1,ll t2){
if(t2==0) return 1;
return cal(t1%mod,t2%mod)*lucas(t1/mod,t2/mod)%mod;
}
int main(){
// int t;
// scanf("%d",&t);
// while(t--){
// scanf("%lld",&n);
// cout<<lucas(n-1,n/2)<<endl;
// }
// scanf("%lld%lld",&n,&m);
cin>>n>>m>>k;
ll N0=lucas(n+m-1,n-1);
ll num=(n+m+1)/(k+1);
for(int i=1;i<=num;i++){
ans+=pow(-1,i+1)*lucas(n,i)*lucas((n+m-1)-(i*(k+1)),(n-1));
ans%=mod;
}
printf("%lld\n",N0-ans);//取模正确吗?
// cout<<lucas(25,5)-(lucas(6,1)*lucas(16,5))+(lucas(6,2)*lucas(7,5));
// printf("Time used= %.2f\n",(double)clock()/CLOCKS_PER_SEC);
// system("pause");
return 0;
}
查看原帖
点赞 9
相关推荐
牛客热帖
正在热议
# 拼多多求职进展汇总 #
235167次浏览 2036人参与
# 实习,投递多份简历没人回复怎么办 #
2437622次浏览 34716人参与
# 北方华创开奖 #
66421次浏览 549人参与
# 25届秋招总结 #
402460次浏览 4034人参与
# ai智能作图 #
23863次浏览 281人参与
# 地方国企笔面经互助 #
6721次浏览 16人参与
# 阿里云管培生offer #
60605次浏览 1755人参与
# 机械求职避坑tips #
22866次浏览 244人参与
# 25届机械人为了秋招做了哪些准备? #
25669次浏览 359人参与
# 软件开发投递记录 #
1479870次浏览 23936人参与
# 我的实习求职记录 #
6126653次浏览 83987人参与
# 歌尔求职进展汇总 #
42831次浏览 294人参与
# 机械人怎么评价今年的华为 #
157469次浏览 1349人参与
# 如果再来一次,你还会选择这个工作吗? #
112939次浏览 1124人参与
# 硬件兄弟们 甩出你的华为奖状 #
78234次浏览 626人参与
# 我在牛爱网找对象 #
74413次浏览 552人参与
# 经纬恒润求职进展汇总 #
99510次浏览 960人参与
# 联想求职进展汇总 #
202949次浏览 1816人参与
# 如果可以,你希望哪个公司来捞你 #
32190次浏览 189人参与
# 阿里求职进展汇总 #
72494次浏览 793人参与
# 牛客租房专区 #
5614次浏览 141人参与
# 许愿池 #
217662次浏览 2545人参与