设x∈Z,存在
{ L = r 1 ( m o d p 1 ) L = r 2 ( m o d p 2 ) L = r 3 ( m o d p 3 ) \left\{ \begin{aligned} L & = r_1 \ (mod \ p1) \\ L & = r_2 \ (mod \ p2) \\ L & = r_3 \ (mod \ p3) \end{aligned} \right. ⎩⎪⎨⎪⎧LLL=r1 (mod p1)=r2 (mod p2)=r3 (mod p3)
令
{ P = p 1 p 2 p 3 . . . p n P i = P p i t i 是 P i 关 于 p i 的 逆 元 , 即 P i t i ≡ 1 ( m o d p i ) \left\{ \begin{aligned} P & = p_1p_2p_3...p_n \\ P_i & = \frac {P}{p_i} \\ t_i & 是P_i关于p_i的逆元,即P_it_i ≡ 1 \ (mod \ p_i) \end{aligned} \right. ⎩⎪⎪⎪⎨⎪⎪⎪⎧PPiti=p1p2p3...pn=piP是Pi关于pi的逆元,即Piti≡1 (mod pi)
则
L = ∑ i = 1 n a i P i t i L = \sum_{i=1}^{n}a_iP_it_i L=∑i=1naiPiti
扩展欧几里得求L即可
例题:https://vjudge.net/problem/LightOJ-1319
板子:
const int N=12+7;
ll p[N],r[N];
int n;
ll P=1;
ll exgcd(ll a,ll b,ll &x,ll &y)//gcd(a,b)=ax+by;返回值gcd
{
if(b==0)
{
x=1;y=0;
return a;
}
ll g=exgcd(b,a%b,x,y);
int z=x;
x=y;
y=z-y*(a/b);
return g;
}
void solve()
{
ll ans=0;
for(int i=1;i<=n;i++)
{
ll pi=P/p[i];
ll ti,y;
exgcd(pi,p[i],ti,y);
ans+=r[i]*pi*ti;
}
cout<<(ans%P+P)%P<<endl;
return ;
}