Lucas定理(板子+理解

最新推荐文章于 2022-05-05 12:41:19 发布
最新推荐文章于 2022-05-05 12:41:19 发布
阅读量485 收藏
点赞数
分类专栏： 排列组合
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_36553623/article/details/77248568
版权

Lucas的定理，定理描述是，如果

 

 

那么得到

 

   

这样然后分别求，采用逆元计算即可。


传送门

C(n,m) mod p

//china no.1
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <vector>
#include <iostream>
#include <string>
#include <map>
#include <stack>
#include <cstring>
#include <queue>
#include <list>
#include <stdio.h>
#include <set>
#include <algorithm>
#include <cstdlib>
#include <cmath>
#include <iomanip>
#include <cctype>
#include <sstream>
#include <functional>
#include <stdlib.h>
#include <time.h>
#include <bitset>
using namespace std;

#define pi acos(-1)
#define PI acos(-1)
#define endl '\n'
#define srand() srand(time(0));
#define me(x,y) memset(x,y,sizeof(x));
#define foreach(it,a) for(__typeof((a).begin()) it=(a).begin();it!=(a).end();it++)
#define close() ios::sync_with_stdio(0); cin.tie(0);
#define FOR(x,n,i) for(int i=x;i<=n;i++)
#define FOr(x,n,i) for(int i=x;i<n;i++)
#define W while
#define sgn(x) ((x) < 0 ? -1 : (x) > 0)
#define bug printf("***********\n");
#define db double
typedef long long LL;
const int INF=0x3f3f3f3f;
const LL LINF=0x3f3f3f3f3f3f3f3fLL;
const int dx[]={-1,0,1,0,1,-1,-1,1};
const int dy[]={0,1,0,-1,-1,1,-1,1};
const int maxn=1e5+10;
const int maxx=1e6+100;
const double EPS=1e-8;
const double eps=1e-8;
const int mod=1e9+7;
template<class T>inline T min(T a,T b,T c) { return min(min(a,b),c);}
template<class T>inline T max(T a,T b,T c) { return max(max(a,b),c);}
template<class T>inline T min(T a,T b,T c,T d) { return min(min(a,b),min(c,d));}
template<class T>inline T max(T a,T b,T c,T d) { return max(max(a,b),max(c,d));}
template <class T>
inline bool scan_d(T &ret){char c;int sgn;if (c = getchar(), c == EOF){return 0;}
while (c != '-' && (c < '0' || c > '9')){c = getchar();}sgn = (c == '-') ? -1 : 1;ret = (c == '-') ? 0 : (c - '0');
while (c = getchar(), c >= '0' && c <= '9'){ret = ret * 10 + (c - '0');}ret *= sgn;return 1;}

inline bool scan_lf(double &num){char in;double Dec=0.1;bool IsN=false,IsD=false;in=getchar();if(in==EOF) return false;
while(in!='-'&&in!='.'&&(in<'0'||in>'9'))in=getchar();if(in=='-'){IsN=true;num=0;}else if(in=='.'){IsD=true;num=0;}
else num=in-'0';if(!IsD){while(in=getchar(),in>='0'&&in<='9'){num*=10;num+=in-'0';}}
if(in!='.'){if(IsN) num=-num;return true;}else{while(in=getchar(),in>='0'&&in<='9'){num+=Dec*(in-'0');Dec*=0.1;}}
if(IsN) num=-num;return true;}

void Out(LL a){if(a < 0) { putchar('-'); a = -a; }if(a >= 10) Out(a / 10);putchar(a % 10 + '0');}
void print(LL a){ Out(a),puts("");}
//freopen( "in.txt" , "r" , stdin );
//freopen( "data.txt" , "w" , stdout );
//cerr << "run time is " << clock() << endl;

//C(m,n) m在上 n在下,p为素数
LL n,m,p;

LL quick_mod(LL a, LL b)
{
    LL ans = 1;
    a %= p;
    while(b)
    {
        if(b & 1)
        {
            ans = ans * a % p;
            b--;
        }
        b >>= 1;
        a = a * a % p;
    }
    return ans;
}

LL C(LL n, LL m)
{
    if(m > n) return 0;
    LL ans = 1;
    for(int i=1; i<=m; i++)
    {
        LL a = (n + i - m) % p;
        LL b = i % p;
        ans = ans * (a * quick_mod(b, p-2) % p) % p;
    }
    return ans;
}

LL Lucas(LL n, LL m)
{
    if(m == 0) return 1;
    return C(n % p, m % p) * Lucas(n / p, m / p) % p;
}

int main()
{
    int T;
    scan_d(T);
    while(T--)
    {
        scan_d(n),scan_d(m),scan_d(p);
        print(Lucas(n,m));
    }
    return 0;
}

和，并且可能为合数

//china no.1
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <vector>
#include <iostream>
#include <string>
#include <map>
#include <stack>
#include <cstring>
#include <queue>
#include <list>
#include <stdio.h>
#include <set>
#include <algorithm>
#include <cstdlib>
#include <cmath>
#include <iomanip>
#include <cctype>
#include <sstream>
#include <functional>
#include <stdlib.h>
#include <time.h>
#include <bitset>
using namespace std;

#define pi acos(-1)
#define PI acos(-1)
#define endl '\n'
#define srand() srand(time(0));
#define me(x,y) memset(x,y,sizeof(x));
#define foreach(it,a) for(__typeof((a).begin()) it=(a).begin();it!=(a).end();it++)
#define close() ios::sync_with_stdio(0); cin.tie(0);
#define FOR(x,n,i) for(int i=x;i<=n;i++)
#define FOr(x,n,i) for(int i=x;i<n;i++)
#define W while
#define sgn(x) ((x) < 0 ? -1 : (x) > 0)
#define bug printf("***********\n");
#define db double
typedef long long LL;
const int INF=0x3f3f3f3f;
const LL LINF=0x3f3f3f3f3f3f3f3fLL;
const int dx[]={-1,0,1,0,1,-1,-1,1};
const int dy[]={0,1,0,-1,-1,1,-1,1};
const int maxn=1e5+10;
const int maxx=1e6+100;
const double EPS=1e-8;
const double eps=1e-8;
const int mod=1e9+7;
template<class T>inline T min(T a,T b,T c) { return min(min(a,b),c);}
template<class T>inline T max(T a,T b,T c) { return max(max(a,b),c);}
template<class T>inline T min(T a,T b,T c,T d) { return min(min(a,b),min(c,d));}
template<class T>inline T max(T a,T b,T c,T d) { return max(max(a,b),max(c,d));}
template <class T>
inline bool scan_d(T &ret){char c;int sgn;if (c = getchar(), c == EOF){return 0;}
while (c != '-' && (c < '0' || c > '9')){c = getchar();}sgn = (c == '-') ? -1 : 1;ret = (c == '-') ? 0 : (c - '0');
while (c = getchar(), c >= '0' && c <= '9'){ret = ret * 10 + (c - '0');}ret *= sgn;return 1;}

inline bool scan_lf(double &num){char in;double Dec=0.1;bool IsN=false,IsD=false;in=getchar();if(in==EOF) return false;
while(in!='-'&&in!='.'&&(in<'0'||in>'9'))in=getchar();if(in=='-'){IsN=true;num=0;}else if(in=='.'){IsD=true;num=0;}
else num=in-'0';if(!IsD){while(in=getchar(),in>='0'&&in<='9'){num*=10;num+=in-'0';}}
if(in!='.'){if(IsN) num=-num;return true;}else{while(in=getchar(),in>='0'&&in<='9'){num+=Dec*(in-'0');Dec*=0.1;}}
if(IsN) num=-num;return true;}

void Out(LL a){if(a < 0) { putchar('-'); a = -a; }if(a >= 10) Out(a / 10);putchar(a % 10 + '0');}
void print(LL a){ Out(a),puts("");}
//freopen( "in.txt" , "r" , stdin );
//freopen( "data.txt" , "w" , stdout );
//cerr << "run time is " << clock() << endl;


const int N = 200005;

bool prime[N];
int p[N];
int cnt;

void isprime()
{
    cnt = 0;
    memset(prime,true,sizeof(prime));
    for(int i=2; i<N; i++)
    {
        if(prime[i])
        {
            p[cnt++] = i;
            for(int j=i+i; j<N; j+=i)
                prime[j] = false;
        }
    }
}

LL quick_mod(LL a,LL b,LL m)
{
    LL ans = 1;
    a %= m;
    while(b)
    {
        if(b & 1)
        {
            ans = ans * a % m;
            b--;
        }
        b >>= 1;
        a = a * a % m;
    }
    return ans;
}

LL Work(LL n,LL p)
{
    LL ans = 0;
    while(n)
    {
        ans += n / p;
        n /= p;
    }
    return ans;
}

LL Solve(LL n,LL m,LL P)
{
    LL ans = 1;
    for(int i=0; i<cnt && p[i]<=n; i++)
    {
        LL x = Work(n, p[i]);
        LL y = Work(n - m, p[i]);
        LL z = Work(m, p[i]);
        x -= (y + z);
        ans *= quick_mod(p[i],x,P);
        ans %= P;
    }
    return ans;
}

int main()
{
    int T;
    isprime();
    cin>>T;
    while(T--)
    {
        LL n,m,P;
        cin>>n>>m>>P;
        n += m - 2;
        m--;
        cout<<Solve(n,m,P)<<endl;
    }
    return 0;
}

题意：中有多少个奇数，其中。

 

分析：其实组合数判断奇偶性有一个优美的结论

          

            如果，那么为奇数，否则为偶数

 

            当然本题要判断的组合数很多，所以不能用上述结论，只能另辟蹊径。由于是判断奇偶性，那么就是判断

     是否为1，利用Lucas定理，先把和化为二进制，这样它们都是01序列了。我们又知道

     。这样中为0的地方对应的中的位置只有一种可能，那就是0。

 

      这样我们可以不用管中为0的地方，只考虑中为1的位置，可以看出，中为1的位置对应的中为0

      或1，其结果都是1，这样答案就是：1<<(二进制表示中1的个数)


//china no.1
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <vector>
#include <iostream>
#include <string>
#include <map>
#include <stack>
#include <cstring>
#include <queue>
#include <list>
#include <stdio.h>
#include <set>
#include <algorithm>
#include <cstdlib>
#include <cmath>
#include <iomanip>
#include <cctype>
#include <sstream>
#include <functional>
#include <stdlib.h>
#include <time.h>
#include <bitset>
using namespace std;

#define pi acos(-1)
#define PI acos(-1)
#define endl '\n'
#define srand() srand(time(0));
#define me(x,y) memset(x,y,sizeof(x));
#define foreach(it,a) for(__typeof((a).begin()) it=(a).begin();it!=(a).end();it++)
#define close() ios::sync_with_stdio(0); cin.tie(0);
#define FOR(x,n,i) for(int i=x;i<=n;i++)
#define FOr(x,n,i) for(int i=x;i<n;i++)
#define W while
#define sgn(x) ((x) < 0 ? -1 : (x) > 0)
#define bug printf("***********\n");
#define db double
typedef long long LL;
const int INF=0x3f3f3f3f;
const LL LINF=0x3f3f3f3f3f3f3f3fLL;
const int dx[]={-1,0,1,0,1,-1,-1,1};
const int dy[]={0,1,0,-1,-1,1,-1,1};
const int maxn=1e5+10;
const int maxx=1e6+100;
const double EPS=1e-8;
const double eps=1e-8;
const int mod=1e9+7;
template<class T>inline T min(T a,T b,T c) { return min(min(a,b),c);}
template<class T>inline T max(T a,T b,T c) { return max(max(a,b),c);}
template<class T>inline T min(T a,T b,T c,T d) { return min(min(a,b),min(c,d));}
template<class T>inline T max(T a,T b,T c,T d) { return max(max(a,b),max(c,d));}
template <class T>
inline bool scan_d(T &ret){char c;int sgn;if (c = getchar(), c == EOF){return 0;}
while (c != '-' && (c < '0' || c > '9')){c = getchar();}sgn = (c == '-') ? -1 : 1;ret = (c == '-') ? 0 : (c - '0');
while (c = getchar(), c >= '0' && c <= '9'){ret = ret * 10 + (c - '0');}ret *= sgn;return 1;}

inline bool scan_lf(double &num){char in;double Dec=0.1;bool IsN=false,IsD=false;in=getchar();if(in==EOF) return false;
while(in!='-'&&in!='.'&&(in<'0'||in>'9'))in=getchar();if(in=='-'){IsN=true;num=0;}else if(in=='.'){IsD=true;num=0;}
else num=in-'0';if(!IsD){while(in=getchar(),in>='0'&&in<='9'){num*=10;num+=in-'0';}}
if(in!='.'){if(IsN) num=-num;return true;}else{while(in=getchar(),in>='0'&&in<='9'){num+=Dec*(in-'0');Dec*=0.1;}}
if(IsN) num=-num;return true;}

void Out(LL a){if(a < 0) { putchar('-'); a = -a; }if(a >= 10) Out(a / 10);putchar(a % 10 + '0');}
void print(LL a){ Out(a),puts("");}
//freopen( "in.txt" , "r" , stdin );
//freopen( "data.txt" , "w" , stdout );
//cerr << "run time is " << clock() << endl;


int n;
int main()
{
    W(scan_d(n))
    {
        int t=1,cnt=0;
        W(n)
        {
            if(n&t)
                cnt++;
            n>>=1;
        }
        print(1<<cnt);
    }
}




twh233
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
卢卡斯定理---板子
three_trees的blog
09-25 192
先贴个代码，这个是p是素数的情况 #include<bits/stdc++.h> using namespace std; typedef long long ll; typedef unsigned long long ull; const int N = 2e5 + 5; ll fast_pow(ll a, ll b, ll p) { ll c = 1; while(b) {if(...
板子：卢卡斯定理
少主大人殿下的博客
03-23 737
卢卡斯算法 原理 用于求大组合数 c(n,m)%mod=C(n%mod,m%mod)*C(n/mod,m/mod)%mod 代码 递归版 LL fac[mod+105],inv[mod+105]; void getInv() { fac[0]=fac[1]=inv[1]=1; for(int i=2;i&lt;mod;i++) { ...
数论——卢卡斯（Lucas）定理（板子）
lymww的博客
05-05 647
文章目录Lucas定理：Lucas板子：Lucas例题ac代码 Lucas定理： CabC_a^bCab​ modmodmod ppp === Cb/pa/p∗Camod  pbmod  pC_{b/p}^{a/p}*C_{a\\\mod p}^{b\\\mod p}Cb/pa/p​∗Camodpbmodp​ modmodmod ppp Lucas板子： int qmi(int a, int k, int p) // 快速幂模板 { int res = 1 % p; while (k)
lucas板子
Bluelanzhan的博客
10-08 291
求c(n,m)mod p ，p是素数； n = qp + r m = sp + t c(n,m) = c(q,s)*c(r,t) (mod p) = c(q,s) * c(q%mod,s%mod); int c(int nn, int mm, int mod) { if(mm &gt; nn) return 0; return 1ll * f[nn] * inv[mm] %mod *...
Lucas定理板子 洛谷P3807
五彩蒟蒻DDD的博客
11-08 346
数论Lucas定理是用来求 C(m，n)%p的值,【p是素数】【辣鸡杭电针对我，在洛谷上能A在杭电上就要T】 HDU3037还涉及到隔板法，详情见此 洛谷直接给了公式。 1、O（n）预处理逆元#include<cstdio> #include<cstring> using namespace std; #define ll long long const int N=1e5+5; /*ll f
数论（组合数求法+Lucas定理+线性求逆元
08-01
数论
lucas定理扩展论文
08-19
讲述了lucas定理如何求C（n，m）%p^k
排列组合及其推论例题内含有lucas定理
08-15
排列组合及其推论例题内含有lucas定理
BakeraMatthews.zip_Lucas-Kanade_lucas
07-15
Lucas-kanade 算法经典文献 总结了该算法的理论和延伸 Lucas-Kanade 20 years on: A unifying framework
fulltext141012017_Facialimage_lucas_
09-29
an efficient feature extraction method based on the Kande-Lucas-Tomasi (KLT) using fast independentcomponent analysis (Fast ICA) & Anthropometric Model as the distance measure is proposed. Each ...
扩展卢卡斯 (板子)
King8611的博客
10-25 256
求c(m,n)%a 不知道原理，不过套板子就对了：   #include&lt;bits/stdc++.h&gt; using namespace std; typedef long long ll; ll exgcd(ll a,ll b,ll &amp;x,ll &amp;y) { if(!b){x=1;y=0;return a;} ll res=exgcd(b,a%...
Lucas–Kanade光流算法
JustRemind的专栏
04-15 4642
作者：gnuhpc@gmail.com  简介：在计算机视觉中，Lucas–Kanade光流算法是一种两帧差分的光流估计算法。它由Bruce D. Lucas 和 Takeo Kanade提出。 光流的概念:(Optical flow or optic flow) 它是一种运动模式，这种运动模式指的是一个物体、表面、边缘在一个视角下由一个观察者（比如眼睛、摄像头等）和背景之间形成
三项式与组合数(lucas板子）
十七的博客
04-20 568
1 、三项式展开 其实可以推广到多项式展开 2、排列组合 %%%大佬博客：https://blog.csdn.net/qq_34531807/article/details/79795261 1、当n,m都很小的时候直接用杨辉三角求 C(n,m)=C(n-1,m)+C(n-1,m-1) 2、利用乘法逆元 （初始化完之后如果需要多次计算组合相...
扩展Lucas板子
include / Hdhd
09-08 224
实在看不惯自己以前的代码风格，改了一下。。。 #include &lt;cstdio&gt; #include &lt;iostream&gt; using namespace std; inline int add(int a,int b,int p) {a+=b;return a&gt;=p?a-p:a;} inline int ch(int a,int b,int p) {return...
数论13——组合数学——lucas 算法（模板0.5 ）
T.J的博客
08-24 271
#include<bits/stdc++.h> #define ll long long using namespace std; const int maxn = 1e6; ll fac[maxn]; void init(ll p)// 求阶乘打表 { fac[0] = 1; for(ll i = 1;i <= p;i++)fac[i] = fac[i-1]*...
ACM_算法_Lucas定理
RaAlGhul的博客
06-24 4550
Lucas定理是用于求解C(n,m)%p的问题 这里小编用一张图： 这张图就很完整的说明了Lucas定理的内容，比较简单，也比较好理解，小编也就不多说了。 #include #include #include #include using namespace std; typedef long long LL; LL Power_mod(LL a, LL b, LL p) {
Lucas定理
leimingfu的博客
08-01 214
简介 Lucas定理用于解决数据较大时的组合数学问题。其中C(n,m)mod p中的模数p为素数。 公式 C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p 。 就是不断递归求C(n%p,m%p)%p ，每次更新n=n/p，m=m/p。C(n,m)%p=(n!/m!*(n-m)!)，这里的除法需要转换为乘法逆元来计算。 算法实现 逆元的计算这里用的是费马小定理算快速幂，也...
Lucas 定理
java开发指南博客 【转载】
08-08 199
转载请注明出处，谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents by---cxlove 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...
java1e6,大组合数取模之lucas定理模板,1<=n<=m<=1e9,1<p<=1e6，p必须为素数
weixin_33380613的博客
03-21 287
typedef long long ll;/**********************************大组合数取模之lucas定理模板,1<=n<=m<=1e9,1输入：C(n,m)%p 调用lucas(n,m,p)复杂度：min(m,p)*log(m)***********************************///ax + by = gcd(a,b)//传...
lucas定理c++
最新发布
09-22
Lucas定理是组合数学中的一种定理，用于计算组合数取模的结果。根据引用和引用的代码，可以看出Lucas定理的计算公式为C(n,m)%p = Lucas(n,m)%p = C(n%p, m%p) * lucas(n/p, m/p) % p。其中，C(n,m)表示组合数，lucas...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值