题目描述
现有两组数字，每组k个，第一组中的数字分别为：a1，a2，...，ak表示，第二组中的数字分别用b1，b2，...，bk表示。其中第二组中的数字是两两互素的。求最小的非负整数n，满足对于任意的i，n - ai能被bi整除。
输入格式
输入数据的第一行是一个整数k，（1 ≤ k ≤ 10）。接下来有两行，第一行是：a1，a2，...，ak，第二行是b1，b2，...，bk
输出格式
输出所求的整数n。
输入输出样例
输入 #1
3 1 2 3 2 3 5
输出 #1
23
说明/提示
所有数据中，第一组数字的绝对值不超过10^9（可能为负数），第二组数字均为不超过6000的正整数，且第二组里所有数的乘积不超过10^18
这其实只是一道“中国剩余定理”的模板题而已，然鹅出题人真的是丧心病狂 用心良苦，偏要设置几个坑让我们跳，很不幸的，我就中招了。
1 #include<bits/stdc++.h> 2 #define ll long long 3 using namespace std; 4 const int N=10000005; 5 int n; 6 int a[11],b[11]; 7 int read() 8 { 9 int x=0,f=1; 10 char ch=getchar(); 11 while(ch<'0'||ch>'9') 12 { 13 if(ch=='-') f=-1; 14 ch=getchar(); 15 } 16 while(ch>='0'&&ch<='9') 17 { 18 x=x*10+ch-'0'; 19 ch=getchar(); 20 } 21 return x*f; 22 } 23 void write(int x) 24 { 25 if(x<0) 26 { 27 putchar('-'); 28 x=-x; 29 } 30 if(x>9) write(x/10); 31 putchar(x%10+'0'); 32 } 33 void exgcd(ll a,ll b,ll &d,ll &x,ll &y) 34 { 35 if(!b) 36 { 37 d=a; 38 x=1; 39 y=0; 40 } 41 else 42 { 43 exgcd(b,a%b,d,x,y); 44 ll t=x;x=y;y=t-a/b*y; 45 } 46 } 47 int main() 48 { 49 int k;k=read();ll m=1,ans=0; 50 for(int i=1;i<=k;i++)a[i]=read(); 51 for(int i=1;i<=k;i++) 52 { 53 b[i]=read(); 54 m*=b[i]; 55 } 56 for(int i=1;i<=k;i++) 57 { 58 ll mi=m/b[i],d,x,y; 59 exgcd(mi,b[i],d,x,y); 60 ans=(ans+a[i]*mi*x)%m; 61 } 62 printf("%lld",(ans+m)%m); 63 return 0; 64
嗯，
代码敲完后自我感觉良好，
直接Ctrl + c 、 Ctrl + v，
按下提交键。
1 #include<bits/stdc++.h> 2 #define ll long long 3 using namespace std; 4 const int N=10000005; 5 int n; 6 int a[11],b[11]; 7 int read() 8 { 9 int x=0,f=1; 10 char ch=getchar(); 11 while(ch<'0'||ch>'9') 12 { 13 if(ch=='-') f=-1; 14 ch=getchar(); 15 } 16 while(ch>='0'&&ch<='9') 17 { 18 x=x*10+ch-'0'; 19 ch=getchar(); 20 } 21 return x*f; 22 } 23 void write(int x) 24 { 25 if(x<0) 26 { 27 putchar('-'); 28 x=-x; 29 } 30 if(x>9) write(x/10); 31 putchar(x%10+'0'); 32 } 33 void exgcd(ll a,ll b,ll &d,ll &x,ll &y) 34 { 35 if(!b) 36 { 37 d=a; 38 x=1; 39 y=0; 40 } 41 else 42 { 43 exgcd(b,a%b,d,x,y); 44 ll t=x;x=y;y=t-a/b*y; 45 } 46 } 47 int main() 48 { 49 int k;k=read();ll m=1,ans=0; 50 for(int i=1;i<=k;i++)a[i]=read(); 51 for(int i=1;i<=k;i++) 52 { 53 b[i]=read(); 54 m*=b[i]; 55 } 56 for(int i=1;i<=k;i++) 57 { 58 ll mi=m/b[i],d,x,y; 59 exgcd(mi,b[i],d,x,y); 60 ans=(ans+a[i]*mi*x)%m; 61 } 62 printf("%lld",(ans+m)%m); 63 return 0; 64 }
1 #include<bits/stdc++.h> 2 #define ll long long 3 using namespace std; 4 const int N=10000005; 5 int n; 6 int a[11],b[11]; 7 int read() 8 { 9 int x=0,f=1; 10 char ch=getchar(); 11 while(ch<'0'||ch>'9') 12 { 13 if(ch=='-') f=-1; 14 ch=getchar(); 15 } 16 while(ch>='0'&&ch<='9') 17 { 18 x=x*10+ch-'0'; 19 ch=getchar(); 20 } 21 return x*f; 22 } 23 void write(int x) 24 { 25 if(x<0) 26 { 27 putchar('-'); 28 x=-x; 29 } 30 if(x>9) write(x/10); 31 putchar(x%10+'0'); 32 } 33 void exgcd(ll a,ll b,ll &d,ll &x,ll &y) 34 { 35 if(!b) 36 { 37 d=a; 38 x=1; 39 y=0; 40 } 41 else 42 { 43 exgcd(b,a%b,d,x,y); 44 ll t=x;x=y;y=t-a/b*y; 45 } 46 } 47 int main() 48 { 49 int k;k=read();ll m=1,ans=0; 50 for(int i=1;i<=k;i++)a[i]=read(); 51 for(int i=1;i<=k;i++) 52 { 53 b[i]=read(); 54 m*=b[i]; 55 } 56 for(int i=1;i<=k;i++) 57 { 58 ll mi=m/b[i],d,x,y; 59 exgcd(mi,b[i],d,x,y); 60 ans=(ans+a[i]*mi*x)%m; 61 } 62 printf("%lld",(ans+m)%m); 63 return 0; 64
结果居然发现——
然后我赶快回去看了几眼代码，感觉没啥毛病，于是又看了下别人的题解——
喔，我手速加了个快速乘，扫了一遍代码，测了遍样例 样例并没软用，再次Ctrl + c 、 Ctrl + v，按下了提交键，这次肯定没问题的吶~
1 #include<bits/stdc++.h> 2 #define ll long long 3 using namespace std; 4 const int N=10000005; 5 int n; 6 int a[11],b[11]; 7 int read() 8 { 9 int x=0,f=1; 10 char ch=getchar(); 11 while(ch<'0'||ch>'9') 12 { 13 if(ch=='-') f=-1; 14 ch=getchar(); 15 } 16 while(ch>='0'&&ch<='9') 17 { 18 x=x*10+ch-'0'; 19 ch=getchar(); 20 } 21 return x*f; 22 } 23 void write(int x) 24 { 25 if(x<0) 26 { 27 putchar('-'); 28 x=-x; 29 } 30 if(x>9) write(x/10); 31 putchar(x%10+'0'); 32 } 33 void exgcd(ll a,ll b,ll &d,ll &x,ll &y) 34 { 35 if(!b) 36 { 37 d=a; 38 x=1; 39 y=0; 40 } 41 else 42 { 43 exgcd(b,a%b,d,x,y); 44 ll t=x;x=y;y=t-a/b*y; 45 } 46 } 47 ll ff(ll a,ll b,ll m) 48 { 49 ll ans=0; 50 while(b) 51 { 52 if(b&1)ans=(ans+a)%m; 53 a=(a+a)%m; 54 b>>=1; 55 } 56 return ans; 57 } 58 int main() 59 { 60 int k;k=read();ll m=1,ans=0; 61 for(int i=1;i<=k;i++)a[i]=read(); 62 for(int i=1;i<=k;i++) 63 { 64 b[i]=read(); 65 m*=b[i]; 66 } 67 for(int i=1;i<=k;i++) 68 { 69 ll mi=m/b[i],d,x,y; 70 exgcd(mi,b[i],d,x,y); 71 x=(x%b[i]+b[i])%b[i]; 72 ans=(ans+ff(ff(mi,x,m),a[i],m))%m; 73 } 74 printf("%lld",(ans+m)%m); 75 return 0; 76 }
1 #include<bits/stdc++.h> 2 #define ll long long 3 using namespace std; 4 const int N=10000005; 5 int n; 6 int a[11],b[11]; 7 int read() 8 { 9 int x=0,f=1; 10 char ch=getchar(); 11 while(ch<'0'||ch>'9') 12 { 13 if(ch=='-') f=-1; 14 ch=getchar(); 15 } 16 while(ch>='0'&&ch<='9') 17 { 18 x=x*10+ch-'0'; 19 ch=getchar(); 20 } 21 return x*f; 22 } 23 void write(int x) 24 { 25 if(x<0) 26 { 27 putchar('-'); 28 x=-x; 29 } 30 if(x>9) write(x/10); 31 putchar(x%10+'0'); 32 } 33 void exgcd(ll a,ll b,ll &d,ll &x,ll &y) 34 { 35 if(!b) 36 { 37 d=a; 38 x=1; 39 y=0; 40 } 41 else 42 { 43 exgcd(b,a%b,d,x,y); 44 ll t=x;x=y;y=t-a/b*y; 45 } 46 } 47 ll ff(ll a,ll b,ll m) 48 { 49 ll ans=0; 50 while(b) 51 { 52 if(b&1)ans=(ans+a)%m; 53 a=(a+a)%m; 54 b>>=1; 55 } 56 return ans; 57 } 58 int main() 59 { 60 int k;k=read();ll m=1,ans=0; 61 for(int i=1;i<=k;i++)a[i]=read(); 62 for(int i=1;i<=k;i++) 63 { 64 b[i]=read(); 65 m*=b[i]; 66 } 67 for(int i=1;i<=k;i++) 68 { 69 ll mi=m/b[i],d,x,y; 70 exgcd(mi,b[i],d,x,y); 71 x=(x%b[i]+b[i])%b[i]; 72 ans=(ans+ff(ff(mi,x,m),a[i],m))%m; 73 } 74 printf("%lld",(ans+m)%m); 75 return 0;
然鹅打脸就是来的这么突然(￣ε(#￣)☆╰╮(￣▽￣///)
我特么第二个点TLE掉是怎么回事？！！
不是，第一遍提交都莫得问题的鸭！我都忍不住要口吐芬芳(｀へ´)
火速赶到题解区翻到了之前没看完的题解——
╮(╯▽╰)╭这毒瘤题，真拿它没办法┑(￣Д ￣)┍
第三次提交，终于满屏全绿（要想生活过得去，做题就得来点绿~
1 #include<bits/stdc++.h> 2 #define ll long long 3 using namespace std; 4 const int N=10000005; 5 int n; 6 int a[11],b[11]; 7 int read() 8 { 9 int x=0,f=1; 10 char ch=getchar(); 11 while(ch<'0'||ch>'9') 12 { 13 if(ch=='-') f=-1; 14 ch=getchar(); 15 } 16 while(ch>='0'&&ch<='9') 17 { 18 x=x*10+ch-'0'; 19 ch=getchar(); 20 } 21 return x*f; 22 } 23 void write(int x) 24 { 25 if(x<0) 26 { 27 putchar('-'); 28 x=-x; 29 } 30 if(x>9) write(x/10); 31 putchar(x%10+'0'); 32 } 33 void exgcd(ll a,ll b,ll &d,ll &x,ll &y) 34 { 35 if(!b) 36 { 37 d=a; 38 x=1; 39 y=0; 40 } 41 else 42 { 43 exgcd(b,a%b,d,x,y); 44 ll t=x;x=y;y=t-a/b*y; 45 } 46 } 47 ll ff(ll a,ll b,ll m) 48 { 49 ll ans=0; 50 while(b) 51 { 52 if(b&1)ans=(ans+a)%m; 53 a=(a+a)%m; 54 b>>=1; 55 } 56 return ans; 57 } 58 int main() 59 { 60 int k;k=read();ll m=1,ans=0; 61 for(int i=1;i<=k;i++)a[i]=read(); 62 for(int i=1;i<=k;i++) 63 { 64 b[i]=read(); 65 m*=b[i]; 66 } 67 for(int i=1;i<=k;i++) 68 { 69 ll mi=m/b[i],d,x,y; 70 exgcd(mi,b[i],d,x,y); 71 x=(x%b[i]+b[i])%b[i]; 72 ans=(ans+ff(ff(mi,x,m),(a[i]+m)%m,m))%m; 73 } 74 printf("%lld",(ans+m)%m); 75 return 0; 76 }