P3868 [TJOI2009]猜数字

最新推荐文章于 2022-03-19 08:24:15 发布

dingqiongliang0363

最新推荐文章于 2022-03-19 08:24:15 发布

阅读量97

点赞数

原文链接：http://www.cnblogs.com/kafuuchino/p/9772530.html

版权

[TJOI2009]猜数字

中国剩余定理

求解i=1 to n : x≡a[i] (mod b[i])的同余方程组

设 t= ∏i=1 to n b[i]

我们先求出 i=1 to n : x≡1 (mod b[i]) ; j=1 to n,j≠i : x≡0 (mod b[j])的解

我们可以把 t/b[i] y 代入 x 求解

ans= ∑i=1 to n : x[i]a[i]

中间过程可能爆long long所以用快速乘

（当然模数不互质的话就只能excrt了）

#include<iostream>
#include<cctype>
using namespace std;
typedef long long ll;
ll n,A[12],B[12],lcm,t=1,X,Y,ans,m;
inline ll kmul(ll x,ll y){ //快速乘
    ll tmp=x*y-(ll)((long double)x/t*y+1.0e-8)*t;
    return tmp>0 ? tmp:tmp+t;
}
void exgcd(ll a,ll b,ll &x,ll &y){
    if(!b) x=1,y=0;
    else exgcd(b,a%b,y,x),y-=x*(a/b);
}
int main(){
    cin>>n;
    for(int i=1;i<=n;++i) cin>>A[i];
    for(int i=1;i<=n;++i) cin>>B[i],t*=B[i];
    for(int i=1;i<=n;++i){
        exgcd(B[i],m=t/B[i],X,Y);
        ans=(ans+kmul(kmul(m,Y),A[i])%t+t)%t;
    }cout<<ans;
}