[TJOI2009]猜数字（洛谷 3868）

题目描述

现有两组数字，每组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 }

 //参考：lahlah 的博客