中国剩余定理与扩展 Lucas定理与扩展学习笔记

最新推荐文章于 2024-02-08 15:16:38 发布

Clove_unique

最新推荐文章于 2024-02-08 15:16:38 发布

阅读量1w

点赞数 27

分类专栏：学习笔记中国剩余定理组合数学

本文链接：https://blog.csdn.net/Clove_unique/article/details/54571216

版权

学习笔记同时被 3 个专栏收录

41 篇文章 2 订阅

订阅专栏

组合数学

35 篇文章 0 订阅

订阅专栏

中国剩余定理

11 篇文章 0 订阅

订阅专栏

中国剩余定理

问题

求同余方程组

⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ x \equiv c 1 (mod m 1) x \equiv c 2 (mod m 2) x \equiv c 3 (mod m 3) . . . x \equiv c k (mod m k)

$\left\{ \begin{array}{c} x\equiv c_1\pmod {m_1}\\ x\equiv c_2\pmod {m_2} \\ x\equiv c_3\pmod {m_3}\\ ...\\ x\equiv c_k\pmod {m_k} \end{array} \right.$
其中满足

(mi,mj)=1,1<=i≠j<=k ( m i , m j ) = 1 , 1 <= i ≠ j <= k $(m_i,m_j)=1,1<=i\neq j<=k$
x的最小正（非负）整数解

结论

令 $M=m_1*m_2*m_3*...*m_k$
则 $x=(\sum\limits_{i=1}^k c_i*{M\over m_i}*inv({M\over m_i},mi))\%M$

证明

a.在模M意义下x只有唯一解 ~~（有多解那还了得）~~
b.令 $n_i$ 满足 $n_i\% {M\over m_i}=0$ 且 $n_i\% m_i=1$ ，则 $N=\sum\limits_{i=1}^k c_i*n_i$ 为原问题的一个解
c.根据上面的式子容易得出 $n_i={M\over m_i}x,n_i=m_iy+1$ ，则 ${M\over m_i}x=m_iy+1$ ，即 ${M\over m_i}x-m_iy=1$
d.由于 $m_i$ 两两互质，所以 ${M\over m_i}$ 与 $m_i$ 也互质，令 $a={M\over m_i},b=m_i$ ，则 $ax-by=1$
e.可以发现我们已经将其化简成扩展欧几里得的基本形式 $ax+by=(a,b)$ ，其中也可以认为x为a在模b意义下的逆元

代码

codevs 3990

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define LL long long

LL k,l,r,n,M,x,y,Min,ans,m[15],c[15];

void exgcd(LL a,LL b,LL &x,LL &y)
{
    if (!b) x=1,y=0;
    else exgcd(b,a%b,y,x),y-=a/b*x;
}
int main()
{
    scanf("%d%lld%lld",&k,&l,&r);
    M=1LL;
    for (int i=1;i<=k;++i)
    {
        scanf("%lld%lld",&m[i],&c[i]);
        M*=m[i];
    }
    for (int i=1;i<=k;++i)
    {
        LL a=M/m[i],b=m[i];
        exgcd(a,b,x,y);
        x=(x%b+b)%b;
        if (!x) x+=b;
        n+=c[i]*a*x;
    }
    n%=M;
    if (!n) n+=M;

    if (r>=n)
        ans=(r-n)/M+1;
    if (l>=n) ans=ans-((l-n)/M+1);
    if ((l-n)%M==0) ++ans;
    if (ans)
    {
        if (l<=n) Min=n;
        else Min=n+((l-n)/M+1)*M;
    }
    printf("%lld\n%lld\n",ans,Min);
}

扩展中国剩余定理

问题

求同余方程组

⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ x \equiv c 1 (mod m 1) x \equiv c 2 (mod m 2) x \equiv c 3 (mod m 3) . . . x \equiv c k (mod m k)

$\left\{ \begin{array}{c} x\equiv c_1\pmod {m_1}\\ x\equiv c_2\pmod {m_2} \\ x\equiv c_3\pmod {m_3}\\ ...\\ x\equiv c_k\pmod {m_k} \end{array} \right.$
x的最小正（非负）整数解

结论

对于两个方程
$x\equiv c_1\pmod {m_1}$
$x\equiv c_2\pmod {m_2}$
将其合并成一个方程，有解条件为 $(m_1,m_2)|(c_2-c_1)$
$m={m_1m_2\over (m_1,m_2)}$
$c=(inv({m_1\over (m_1,m_2)},{m_2\over (m_1,m_2)})*{(c_2-c_1)\over (m_1,m_2)})\%{m_2\over (m_1,m_2)}*m_1+c_1$
最终得出一个式子 $x\equiv c\pmod {m}$
$x=c\%m$ 即为原问题的一个解

证明

a.将两个方程写成 $x=m_1k_1+c_1,x=m_2k_2+c_2$ ，两式联立得 $m_1k_1+c_1=m_2k_2+c_2$ ，移项 $m_1k_1=c_2-c_1+m_2k_2$ ，
b.根据贝祖定理，以上等式有解充要条件为 $(m_1,m_2)|(c_2-c_1)$
c.将等式两边同除 $(m_1,m_2)$ 得 ${m_1\over(m_1,m_2)}*k_1={(c_2-c_1)\over (m_1,m_2)}+{m_2\over (m_1,m_2)}*k_2$ ，即 ${m_1\over(m_1,m_2)}k_1\equiv {(c_2-c_1)\over (m_1,m_2)}\pmod {{m_2\over(m_1,m_2)}}$ ，进一步化简 $k_1\equiv inv({m_1\over(m_1,m_2)},{m_2\over (m_1,m_2)})*{(c_2-c_1)\over (m_1,m_2)}\pmod {{m_2\over(m_1,m_2)}}$ ， $k_1=inv({m_1\over(m_1,m_2)},{m_2\over (m_1,m_2)})*{(c_2-c_1)\over (m_1,m_2)}+y{{m_2\over (m_1,m_2)}}$
d.将其回代 $x=m_1k_1+c_1$ ，得 $x=inv({m_1\over(m_1,m_2)},{m_2\over (m_1,m_2)})*{(c_2-c_1)\over (m_1,m_2)}*m_1+y{{m_1m_2\over (m_1,m_2)}}+c_1$ ，即 $x\equiv inv({m_1\over(m_1,m_2)},{m_2\over (m_1,m_2)})*{(c_2-c_1)\over (m_1,m_2)}*m_1+c_1\pmod {{m_1m_2\over (m_1,m_2)}}$
e.至此我们又将其化简成了 $x\equiv c\pmod m$ 的形式，其中 $m={m_1m_2\over (m_1,m_2)}$ ， $c=(inv({m_1\over (m_1,m_2)},{m_2\over (m_1,m_2)})*{(c_2-c_1)\over (m_1,m_2)})\%{m_2\over (m_1,m_2)}*m_1+c_1$

代码

pku 2891

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define LL long long
#define N 1005

int k;
LL c[N],m[N],c1,c2,m1,m2,t;
bool flag;

LL gcd(LL a,LL b)
{
    if (!b) return a;
    else return gcd(b,a%b);
}
void exgcd(LL a,LL b,LL &x,LL &y)
{
    if (!b) x=1LL,y=0LL;
    else exgcd(b,a%b,y,x),y-=a/b*x;
}
LL inv(LL a,LL b)
{
    LL x=0LL,y=0LL;
    exgcd(a,b,x,y);
    x=(x%b+b)%b;
    if (!x) x+=b;
    return x;
}
int main()
{
    while (~scanf("%d",&k))
    {
        flag=true;
        for (int i=1;i<=k;++i)
            scanf("%I64d%I64d",&m[i],&c[i]);
        for (int i=2;i<=k;++i)
        {
            m1=m[i-1],m2=m[i],c1=c[i-1],c2=c[i];
            t=gcd(m1,m2);
            if ((c2-c1)%t!=0) {flag=false;break;}
            m[i]=m1*m2/t;
            c[i]=inv(m1/t,m2/t)*((c2-c1)/t)%(m2/t)*m1+c1;
            c[i]=(c[i]%m[i]+m[i])%m[i];
        }
        if (!flag) puts("-1");
        else printf("%I64d\n",c[k]);
    }
}

Lucas定理

问题

求 $C_n^m\%p$ ，其中p为质数

结论

令
$n=n_k*p^k+n_{k-1}*p^{k-1}+...+n_2*p^2+n_1*p+n_0$
$m=m_k*p^k+m_{k-1}*p^{k-1}+...+m_2*p^2+m_1*p+m_0$
则 $C_n^m=\prod\limits_{i=0}^k C_{n_i}^{m_i}$

证明

~~去问卢卡斯~~
听说要用二项式定理什么的
我怎么可能会这种东西
~~记住就行了~~

代码

zoj 3557

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define LL long long

LL n,m,Mod;

LL fast_pow(LL a,LL p)
{
    LL ans=1LL;
    for (;p;p>>=1,a=a*a%Mod)
        if (p&1)
            ans=ans*a%Mod;
    return ans;
}
LL inv(LL x)
{
    return fast_pow(x,Mod-2);
}
LL C(LL n,LL m)
{
    if (m>n) return 0LL;
    LL up=1LL,down=1LL;
    for (LL i=n-m+1;i<=n;++i) up=up*i%Mod;
    for (LL i=1;i<=m;++i) down=down*i%Mod;
    return up*inv(down)%Mod;
}
LL lucas(LL n,LL m)
{
    if (m>n) return 0LL;
    LL ans=1;
    for (;m;n/=Mod,m/=Mod)
        ans=ans*C(n%Mod,m%Mod)%Mod;
    return ans;
}
int main()
{
    while (~scanf("%lld%lld%lld",&n,&m,&Mod))
        printf("%lld\n",lucas(n-m+1,m));
}

扩展Lucas定理

问题

求 $C_n^m\%p$

结论

令
$p={p_1}^{k_1}*{p_2}^{k_2}*...*{p_q}^{k_q}$
列出同余方程组

⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ a n s \equiv c 1 (mod p 1 k 1) a n s \equiv c 2 (mod p 2 k 2) . . . a n s \equiv c q (mod p q k q)

$\left\{ \begin{array}{c} ans\equiv c_1\pmod {{p_1}^{k_1}}\\ ans\equiv c_2\pmod {{p_2}^{k_2}}\\ ...\\ ans\equiv c_q\pmod {{p_q}^{k_q}}\\ \end{array} \right.$
其中

c1...cq c 1 . . . c q $c_1...c_q$ 是对于每一个

Cmn%piki C n m % p i k i $C_n^m\%{p_i}^{k_i}$ 求出的答案
然后根据中国剩余定理合并
可见

piki p i k i ${p_i}^{k_i}$ 并不是质数，而是一个质数的幂的形式
对于如何求

Cmn%piki C n m % p i k i $C_n^m\% {p_i}^{k_i}$ 将在证明中给出

证明

a.由于同于方程组在模 $p={p_1}^{k_1}*{p_2}^{k_2}*...*{p_q}^{k_q}$ 意义下有唯一解，可以证明上面做法的正确性
b.由于 $p={p_1}^{k_1}*{p_2}^{k_2}*...*{p_q}^{k_q}$ 实质上是将p进行质因数分解，所以 ${p_i}^{k_i}$ 满足两两互质，可以直接用中国剩余定理合并
c.对于如何求 $C_n^m\%{p_i}^{k_i}$
根据 $C_n^m={n!\over m!(n-m)!}$ ，只要分别求出 $n!\% {p_i}^{k_i},m!\% {p_i}^{k_i},(n-m)!\% {p_i}^{k_i}$ 就可以通过逆元求出 $C_n^m\%{p_i}^{k_i}$
d.对于如何求 $n!\%{p_i}^{k_i}$
我们以 $n=19,p_i=3,k_i=2$ 为例
$n!=1*2*3*4*5*6*7*8*9*10*11*12*13*14*15*16*17*18*19$
$=(1*2*4*5*7*8*10*11*13*14*16*17*19)*3^6*(1*2*3*4*5*6)$
根据这个例子发现，求解n!可以分为3部分：第一部分是 $p_i$ 的幂的部分，也就是 $3^6$ 即 ${p_i}^{\lfloor{n\over p_i}\rfloor}$ ，可以直接求解；第二部分是一个新的阶乘，也就是6!即 ${\lfloor{n\over p_i}\rfloor}!$ ，可以递归下去求解；第三部分是除前两部分之外剩下的数
考虑第三部分如何求解
发现第三部分在模 ${p_i}^{k_i}$ 意义下是以 ${p_i}^{k_i}$ 为周期的，即 $(1*2*4*5*7*8)\equiv(10*11*13*14*16*17)\pmod {{p_i}^{k_i}}$ ，所以只求 ${p_i}^{k_i}$ 长度的即可；但是还剩下一个孤立的19，可以发现剩下孤立的数长度不会超过 ${p_i}^{k_i}$ ，只需要暴力求解即可
e.最后一个问题是对于求出的 $m!\%{p_i}^{k_i}$ 和 $(n-m)!\%{p_i}^{k_i}$ 有可能与 ${p_i}^{k_i}$ 不互质，无法求逆元
所以要将 $m!\%{p_i}^{k_i}$ 和 $(n-m)!\%{p_i}^{k_i}$ 中质因子 $p_i$ 先全部除去，求出逆元后再全部乘回去
计算n!中质因子p的个数x的公式为 $x=\lfloor{n\over p}\rfloor+\lfloor{n\over p^2}\rfloor+\lfloor{n\over p^3}\rfloor+...$
递推式也可以写为 $f(n)=f(\lfloor{n\over p}\rfloor)+\lfloor{n\over p}\rfloor$

代码

codeforces2015ICL,Finals,Div.1#J

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define LL long long

LL n,m,MOD,ans;

LL fast_pow(LL a,LL p,LL Mod)
{
    LL ans=1LL;
    for (;p;p>>=1,a=a*a%Mod)
        if (p&1)
            ans=ans*a%Mod;
    return ans;
}
void exgcd(LL a,LL b,LL &x,LL &y)
{
    if (!b) x=1LL,y=0LL;
    else exgcd(b,a%b,y,x),y-=a/b*x;
}
LL inv(LL A,LL Mod)
{
    if (!A) return 0LL;
    LL a=A,b=Mod,x=0LL,y=0LL;
    exgcd(a,b,x,y);
    x=((x%b)+b)%b;
    if (!x) x+=b;
    return x;
}
LL Mul(LL n,LL pi,LL pk)
{
    if (!n) return 1LL;
    LL ans=1LL;
    if (n/pk)
    {
        for (LL i=2;i<=pk;++i)
            if (i%pi) ans=ans*i%pk;
        ans=fast_pow(ans,n/pk,pk);
    }
    for (LL i=2;i<=n%pk;++i)
        if (i%pi) ans=ans*i%pk;
    return ans*Mul(n/pi,pi,pk)%pk;
}
LL C(LL n,LL m,LL Mod,LL pi,LL pk)
{
    if (m>n) return 0LL;
    LL a=Mul(n,pi,pk),b=Mul(m,pi,pk),c=Mul(n-m,pi,pk);
    LL k=0LL,ans;
    for (LL i=n;i;i/=pi) k+=i/pi;
    for (LL i=m;i;i/=pi) k-=i/pi;
    for (LL i=n-m;i;i/=pi) k-=i/pi;
    ans=a*inv(b,pk)%pk*inv(c,pk)%pk*fast_pow(pi,k,pk)%pk;
    return ans*(Mod/pk)%Mod*inv(Mod/pk,pk)%Mod;
}
int main()
{
    scanf("%I64d%I64d%I64d",&n,&m,&MOD);
    for (LL x=MOD,i=2;i<=MOD;++i)
        if (x%i==0)
        {
            LL pk=1LL;
            while (x%i==0) pk*=i,x/=i;
            ans=(ans+C(n,m,MOD,i,pk))%MOD;
        }
    printf("%I64d\n",ans);
}

Clove_unique

关注

27
点赞
踩
37

收藏

觉得还不错? 一键收藏
14
评论
中国剩余定理与扩展 Lucas定理与扩展学习笔记

中国剩余定理问题求同余方程组 ⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪x≡c1(modm1)x≡c2(modm2)x≡c3(modm3)...x≡ck(modmk){x≡c1(modm1)x≡c2(modm2)x≡c3(modm3)...x≡ck(modmk)\left\{ \begin{array}{c}x\equiv c_1\pmod {m_1}\\ x\equiv c_2\pmo...
复制链接

扫一扫