codevs3990 中国余数定理 2（中国剩余定理）

最新推荐文章于 2021-07-29 23:07:25 发布

Icefox_zhx

最新推荐文章于 2021-07-29 23:07:25 发布

阅读量206

点赞数

分类专栏：其他oj ---------数学--------- CRT 扩展欧几里得算法

本文链接：https://blog.csdn.net/Icefox_zhx/article/details/79815690

版权

其他oj 同时被 3 个专栏收录

336 篇文章 0 订阅

订阅专栏

---------数学---------

209 篇文章 0 订阅

订阅专栏

CRT

7 篇文章 0 订阅

订阅专栏

中国剩余定理（CRT）板子题。

中国剩余定理就是用来求解同余方程组

⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ 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\neq j<=k$
x的最小非负数解。

结论：
令 $M=m_1*m_2*m_3*...*m_k$
则 $x=(\sum\limits_{i=1}^kc_i*\frac{M}{m_i}*inv(\frac{M}{m_i},m_i))\%M$
证明：portal

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define ld long double
#define inf 0x3f3f3f3f
#define N 15
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline ll read(){
    ll x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m[N],c[N];ll L,R,M=1,num=0;
inline void exgcd(ll a,ll b,ll &x,ll &y){
    if(!b){x=1;y=0;return;}
    exgcd(b,a%b,x,y);ll t=x;x=y;y=t-a/b*y;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();L=read();R=read();
    for(int i=1;i<=n;++i){
        M*=(m[i]=read());c[i]=read();
    }ll ans=0,x,y;
    for(int i=1;i<=n;++i){
        ll a=M/m[i],b=m[i];exgcd(a,b,x,y);x=(x%b+b)%b;
        (ans+=c[i]*a%M*x)%=M;
    }if(ans<L) ans+=((L-1-ans)/M+1)*M;
    if(ans>R){printf("0\n0\n");return 0;}
    num=(R-ans)/M+1;
    printf("%lld\n%lld\n",num,ans);
    return 0;
}