CodeVS3040 中国余数定理 1

最新推荐文章于 2022-03-15 18:31:17 发布

*ACoder*

最新推荐文章于 2022-03-15 18:31:17 发布

阅读量387

点赞数

分类专栏： # 中国剩余定理

本文链接：https://blog.csdn.net/FSAHFGSADHSAKNDAS/article/details/60960492

版权

中国剩余定理专栏收录该内容

6 篇文章 0 订阅

订阅专栏

链接

　　http://codevs.cn/problem/3040/

题解

　　中国剩余定理。
　　对于方程组
　　

⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ x = a 1 (m o d m 1) x = a 2 (m o d m 2) . . . x = a n (m o d m n)

$\begin{cases}x=a_1(mod\ m_1)\\x=a_2(mod\ m_2)\\...\\x=a_n(mod\ m_n)\end{cases}$
　　解的构造方法是
　　

M = m 1 m 2 m 3 . . . m k M i = M m i a n s = \sum i = 1 n a i M i M - 1 i

$M=m_1m_2m_3...m_k\\M_i=\frac{M}{m_i}\\ans=\sum_{i=1}^na_iM_iM_i^{-1}$
　　其中

M−1i $M_i^{-1}$ 是

Mi $M_i$ 在模

mi $m_i$ 意义下的逆元
　　

ans mod M $ans\ mod\ M$ 就是最小的正整数解

代码

//中国剩余定理
#include <cstdio>
#include <algorithm>
#define maxn 100
#define ll long long
using namespace std;
ll m[maxn], a[maxn];
void exgcd(ll a, ll b, ll &x, ll &y)
{
    if(!b){x=1,y=0;return;}
    ll xx, yy;
    exgcd(b,a%b,xx,yy);
    x=yy,y=xx-a/b*yy;
}
ll inv(ll a, ll p)
{
    ll x, y;
    exgcd(a,p,x,y);
    return (x+p)%p;
}
ll crt(ll *a, ll *m, ll n)
{
    ll M=1, Mi, x, ans=0, i;
    for(i=1;i<=n;i++)M=M*m[i];
    for(i=1;i<=n;i++)
    {
        Mi=M/m[i];
        x=inv(Mi,m[i]);
        ans=ans+a[i]*Mi*x;
    }
    return ans%M;
}
int main()
{
    ll k, ans;
    scanf("%lld",&k);
    a[1]=2, a[2]=3, a[3]=2;
    m[1]=3, m[2]=5, m[3]=7;
    ans=crt(a,m,3);
    ans+=(k-1)*105;
    printf("%lld\n",ans);
    return 0;
}