HDU 3037 (大数&LUCAS定理)
2016-08-04 22:23
302 查看
题目大意:求在n棵树上摘不超过m颗豆子的方案,结果对p取模。
解题思路:
题目可以转换成 x1+x2+……+xn=m 有多少组解,m在题中可以取0~m。
利用插板法可以得出x1+x2+……+xn=m解的个数为C(n+m-1,m);
则题目解的个数可以转换成求 sum=C(n+m-1,0)+C(n+m-1,1)+C(n+m-1,2)……+C(n+m-1,m)
利用公式C(n,r)=C(n-1,r)+C(n-1,r-1) == > sum=C(n+m,m)。
现在就是要求C(n+m,m)%p。
因为n,m很大,这里可以直接套用Lucas定理的模板即可。
Lucas(n,m,p)=C(n%p,m%p,p)*Lucas(n/p,m/p,p); ///这里可以采用对n分段递归求解,
Lucas(x,0,p)=1;
AC代码:(LUCAS函数可直接使用)
解题思路:
题目可以转换成 x1+x2+……+xn=m 有多少组解,m在题中可以取0~m。
利用插板法可以得出x1+x2+……+xn=m解的个数为C(n+m-1,m);
则题目解的个数可以转换成求 sum=C(n+m-1,0)+C(n+m-1,1)+C(n+m-1,2)……+C(n+m-1,m)
利用公式C(n,r)=C(n-1,r)+C(n-1,r-1) == > sum=C(n+m,m)。
现在就是要求C(n+m,m)%p。
因为n,m很大,这里可以直接套用Lucas定理的模板即可。
Lucas(n,m,p)=C(n%p,m%p,p)*Lucas(n/p,m/p,p); ///这里可以采用对n分段递归求解,
Lucas(x,0,p)=1;
AC代码:(LUCAS函数可直接使用)
#include <iostream> #include <algorithm> #include <cstring> #include <cmath> #include <complex> #include <cstdio> using namespace std; #define ll long long #define mod 110119 //ll n,m,r,ca=0,cnt,a,b,rig,up; #define N 111000 ll mod_pow(ll a,ll n,ll p) { ll ret=1,A=a; while(n) { if (n & 1) ret=(ret*A)%p; A=(A*A)%p; n>>=1; } return ret; } ll factorial ; void init(ll p) { factorial[0] = 1; for(int i = 1;i <= p;i++) factorial[i] = factorial[i-1]*i%p; } ll Lucas(ll a,ll k,ll p) //求C(n,m)%p p最大为10^5。a,b可以很大! a个数中挑k个的组合数 { ll re = 1; while(a && k) { ll aa = a%p;ll bb = k%p; if(aa < bb) return 0; //这个是最后的改动! re = re*factorial[aa]*mod_pow(factorial[bb]*factorial[aa-bb]%p,p-2,p)%p;//这儿的求逆不可先处理 a /= p; k /= p; } return re; } int main() { int t; ll a,b,p; cin >> t; while(t--) { scanf("%I64d%I64d%I64d",&a,&b,&p); init(p); ll ans=Lucas(a+b,b,p); printf("%I64d\n",ans); } return 0; }
相关文章推荐
- HDU 3037 Saving Beans 大组合数 lucas定理
- HDU 3037 Saving Beans(Lucas定理模板题)
- hdu 3037 Saving Beans(Lucas定理)
- hdu 3037 Lucas定理
- HDU 3037 Saving Beans (Lucas定理)
- HDU 3037 Saving Beans(Lucas定理的直接应用)
- 【HDU 3037】大数组合取模之Lucas定理+扩展欧几里得求逆元与不定方程一类问题
- [HDU 3461] Saving Beans & 组合计数Lucas定理模板
- hdu 3037 lucas定理
- hdu 3037 Saving Beans(lucas定理模板)
- hdu - 4349 - Xiao Ming's Hope - 大大的Lucas定理 && 小小的乘法逆元
- hdu 3037 saving beans (lucas定理)
- 【HDU 3037】大数组合取模之Lucas定理
- HDU 4349 组合数的奇数个数-杨辉三角&Lucas定理
- [ACM] hdu 3037 Saving Beans (Lucas定理,组合数取模)
- [Lucas定理] hdu 3037 Saving Beans
- hdu 3037 lucas定理
- hdu 3037 Saving Beans 【大组合数取模-Lucas定理+逆元取模】
- HDU 3037 Saving Beans (Lucas定理)
- hdu 4349 Xiao Ming's Hope,lucas定理