Lucas定理 A、B是非负整数,p是质数。AB写成p进制:A=a[n]a[n-1]...a[0],B=b[n]b[n-1]...b[0]。 则组合数C(A,B)与C(a[n],b[n])*C(a[n-1],b[n-1])*...*C(a[0],b[0])modp同 即:Lucas(n,m,p...
分类:
编程语言 时间:
2014-11-02 10:42:32
阅读次数:
482
题目大意:发上来就过不了审核了……总之大意就是求C(n,m) mod 10007 m,n∈[1,2*10^8]
卢卡斯定理:C(n,m)=C(n%p,m%p)*C(n/p,m/p) mod p 要求p是质数
其中n%p可能会小于m%p 这种情况下直接返回0即可
证明去问卢卡斯 我不知道
#include
#include
#include
#include
#define p 10007
...
分类:
其他好文 时间:
2014-10-27 15:41:33
阅读次数:
198
解题思路:
直接求C(n+m , m) % p , 因为n , m ,p都很大,所以要用Lucas定理来解决大组合数取模的问题。
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#de...
分类:
其他好文 时间:
2014-10-17 11:54:26
阅读次数:
185
输入p n 求杨辉三角的n+1不能被p整除的数有多少个
Lucas定理:
A、B是非负整数,p是质数。AB写成p进制:A=a[n]a[n-1]...a[0],B=b[n]b[n-1]...b[0]。
则组合数C(A,B)与C(a[n],b[n])*C(a[n-1],b[n-1])*...*C(a[0],b[0]) mod p同余
即:Lucas(n,m,p)=c(...
分类:
其他好文 时间:
2014-10-11 15:25:35
阅读次数:
226
喵喵的神?数
Time Limit: 1 Sec Memory Limit: 128 MB
Description
喵喵对组合数比较感兴趣,并且对计算组合数非常在行。同时为了追求有后宫的素质的生活,喵喵每天都要研究
质数。
我们先来复习一下什么叫做组合数。对于正整数P、T
...
分类:
其他好文 时间:
2014-10-06 17:32:50
阅读次数:
221
插空法 大组合数取余
#include
#include
using namespace std;
typedef long long LL;
//求整数x和y,使得ax+by=d, 且|x|+|y|最小。其中d=gcd(a,b)
void gcd(LL a, LL b, LL& d, LL& x, LL& y)
{
if(!b)
{
d = a;
x = 1;
y = ...
分类:
其他好文 时间:
2014-09-29 19:37:31
阅读次数:
218
Lucas定理这里有详细的证明。其实就是针对n, m很大时,要求组合数C(n, m) % p, 一般来说如果p 2 #include 3 #include 4 using namespace std; 5 6 #define N 100010 7 8 long long mod_pow(i...
分类:
其他好文 时间:
2014-09-28 22:10:55
阅读次数:
192
直接lucas降到10w以内搞组合数
#include
#include
typedef __int64 LL;
LL f[110010];
LL pow(LL a, LL b, LL c)
{
LL ans = 1;
while(b)
{
if(b&1)
ans = (ans*a) % c;
b >>= 1;
a = (a*a) % c;
}
return an...
分类:
其他好文 时间:
2014-09-28 21:49:55
阅读次数:
216
/*求在n棵树上摘不超过m颗豆子的方案,结果对p取模。
求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;
将n,m分解变小之后问题又转换成了求C(a/b)%p。
而C(a,b) =a! / ( b! *...
分类:
其他好文 时间:
2014-09-26 22:58:09
阅读次数:
176
考虑加多一颗树,这样的话当加的树放了k(0#include #include #include #include using namespace std;typedef long long ll;ll n,m,p;ll POW(ll x,ll n,ll p){ ll res=1; whi...
分类:
其他好文 时间:
2014-09-23 17:04:34
阅读次数:
264