Luogu P3868 猜数字 题解报告

题目传送门

【题目大意】

【思路分析】

是中国剩余定理板子题嗷QwQ

中间会爆$long\ long$，简单一点就用$\_int128$。我为了练手打了个龟速乘，结果$T$了QAQ，改了好久才过。

【代码实现】

 1 #include
 2 #include
 3 #include
 4 #include
 5 #include
 6 #include
 7 #define g() getchar()
 8 #define rg register
 9 #define go(i,a,b) for(rg int i=a;i<=b;i++)
10 #define back(i,a,b) for(rg int i=a;i>=b;i--)
11 #define db double
12 #define ll long long
13 #define il inline
14 #define pf printf
15 #define mem(a,b) memset(a,b,sizeof(a))
16 using namespace std;
17 int fr(){
18     int w=0,q=1;
19     char ch=g();
20     while(ch<'0'||ch>'9'){
21         if(ch=='-') q=-1;
22         ch=g();
23     }
24     while(ch>='0'&&ch<='9') w=(w<<1)+(w<<3)+ch-'0',ch=g();
25     return w*q;
26 }
27 int n,a[12],b[12];
28 ll B=1,ans,x,y;
29 il ll mul(ll x,ll y){
30     ll as=0;
31     while(y){
32         if(y&1) as=(as+x)%B;
33         x=(x+x)%B;y>>=1;
34     }
35     return (as+B)%B;
36 }
37 il void exgcd(ll aa,ll bb,ll &x,ll &y){
38     if(!bb) {x=1;y=0;return;}
39     exgcd(bb,aa%bb,x,y);
40     ll z=x;x=y;y=z-y*(aa/bb);
41     return;
42 }
43 int main(){
44     //freopen("","r",stdin);
45     //freopen("","w",stdout);
46     n=fr();
47     go(i,1,n) a[i]=fr();
48     go(i,1,n) b[i]=fr(),B*=b[i];
49     go(i,1,n){
50         a[i]=(a[i]%b[i]+b[i])%b[i];
51         ll bb=B/b[i];
52         exgcd(bb,(ll)b[i],x,y);
53         x=(x%b[i]+b[i])%b[i];
54         ans=(ans+mul(bb,mul(x,a[i]))+B)%B;
55     }
56     pf("%lld\n",ans);
57     return 0;
58 }

代码戳这里