codevs1515 瞎搞+Lucas

xiaoxiao2021-02-28  134

题目分析

首先我们打出一张地图: 1 1 1 1 1 1 2 3 4 5 1 3 6 10 15 1 4 10 20 35 1 5 15 35 70 别告诉我你对这张图没有任何想法,如果你确实没有想法,说明你与组合数还不够亲近。。。事实上,把图旋转45度后可以得到一个杨辉三角(意会可证),因此我们发现第x行y列的数就是 Cxx+y 。 然后我们脑补可知,假如m>n,我们就先沿着长边走m步,然后再沿着短边走。 所以得到答案就是 m+Cnm+n+1 Cm+1m+n+1=ni=0Cmm+i=Cnm+n+1 ,欲知如何证明,请看我曾经一篇总结组合数性质的博客:走你╭(′▽`)╯。) 好,然后你会发现这题其实可以直接上公式求的。。。然而。。。你就WA了4个点。。。 你会迷茫。。。会怀疑人生,怀疑鲁迅,怀疑枣树。。。 于是,你轻轻地将long long改成了unsigned long long,然后,你得到了这个点的正确结果。。。 你会提交,会发现,WA的点更多了。。。 你会怀疑出题人,怀疑世界,怀疑鲁迅,怀疑枣树。。。 于是你会去搜题解,你发现别人都是用的Lucas,并且没有开unsigned long long。 这时候,你会写Lucas,并A掉这题。 (boshi体结束) 天哪这题数据炸了long long!!!因为Lucas的蜜汁原因会有4个点不会炸!!!什么鬼标程。

代码

#include<iostream> #include<cstdio> #include<climits> #include<cstring> #include<algorithm> using namespace std; #define LL long long LL mod=1000000007; LL ksm(LL x,LL y){ LL re=1; while(y){ if(y&1)re=re*x%mod; x=x*x%mod,y>>=1; } return re; } LL C(LL m,LL n){ LL s1=1,s2=1; for(LL i=m-n+1;i<=m;i++)s1=s1*i%mod; for(LL i=1;i<=n;i++)s2=s2*i%mod; return s1*ksm(s2,mod-2)%mod; } LL lucas(LL m,LL n){ if(!n)return 1; return lucas(m/mod,n/mod)*C(m%mod,n%mod)%mod; } int main() { LL n,m; scanf("%lld%lld",&n,&m); if(n>m)swap(n,m); printf("%lld",(m+lucas(n+m+1,n))%mod); return 0; }
转载请注明原文地址: https://www.6miu.com/read-62695.html

最新回复(0)