题目描述

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

这其实只是一道“中国剩余定理”的模板题而已，然鹅出题人真的是丧心病狂 用心良苦，偏要设置几个坑让我们跳，很不幸的，我就中招了。

 #include<bits/stdc++.h>

 #define ll long long

 using namespace std;

 const int N=;

 int n;

 int a[],b[];

 int read()

 {

     int x=,f=;

     char ch=getchar();

     while(ch<''||ch>'')

     {

         if(ch=='-') f=-;

         ch=getchar();

     }

     while(ch>=''&&ch<='')

     {

         x=x*+ch-'';

         ch=getchar();

     }

     return x*f;

 }

 void write(int x)

 {

     if(x<)

     {

         putchar('-');

         x=-x;

     }

     if(x>) write(x/);

     putchar(x%+'');

 }

 void exgcd(ll a,ll b,ll &d,ll &x,ll &y)

 {

     if(!b)

     {

         d=a;

         x=;

         y=;

     }

     else

     {

         exgcd(b,a%b,d,x,y);

         ll t=x;x=y;y=t-a/b*y;

     }

 }

 int main()

 {

     int k;k=read();ll m=,ans=;

     for(int i=;i<=k;i++)a[i]=read();

     for(int i=;i<=k;i++)

     {

         b[i]=read();

         m*=b[i];

     }

     for(int i=;i<=k;i++)

     {

         ll mi=m/b[i],d,x,y;

         exgcd(mi,b[i],d,x,y);

         ans=(ans+a[i]*mi*x)%m;

     }

     printf("%lld",(ans+m)%m);

     return ;

嗯，

代码敲完后自我感觉良好，

直接Ctrl + c 、 Ctrl + v，

按下提交键。

 #include<bits/stdc++.h>

 #define ll long long

 using namespace std;

 const int N=;

 int n;

 int a[],b[];

 int read()

 {

     int x=,f=;

     char ch=getchar();

     while(ch<''||ch>'')

     {

         if(ch=='-') f=-;

         ch=getchar();

     }

     while(ch>=''&&ch<='')

     {

         x=x*+ch-'';

         ch=getchar();

     }

     return x*f;

 }

 void write(int x)

 {

     if(x<)

     {

         putchar('-');

         x=-x;

     }

     if(x>) write(x/);

     putchar(x%+'');

 }

 void exgcd(ll a,ll b,ll &d,ll &x,ll &y)

 {

     if(!b)

     {

         d=a;

         x=;

         y=;

     }

     else

     {

         exgcd(b,a%b,d,x,y);

         ll t=x;x=y;y=t-a/b*y;

     }

 }

 int main()

 {

     int k;k=read();ll m=,ans=;

     for(int i=;i<=k;i++)a[i]=read();

     for(int i=;i<=k;i++)

     {

         b[i]=read();

         m*=b[i];

     }

     for(int i=;i<=k;i++)

     {

         ll mi=m/b[i],d,x,y;

         exgcd(mi,b[i],d,x,y);

         ans=(ans+a[i]*mi*x)%m;

     }

     printf("%lld",(ans+m)%m);

     return ;

 }
 #include<bits/stdc++.h>

 #define ll long long

 using namespace std;

 const int N=;

 int n;

 int a[],b[];

 int read()

 {

     int x=,f=;

     char ch=getchar();

     while(ch<''||ch>'')

     {

         if(ch=='-') f=-;

         ch=getchar();

     }

     while(ch>=''&&ch<='')

     {

         x=x*+ch-'';

         ch=getchar();

     }

     return x*f;

 }

 void write(int x)

 {

     if(x<)

     {

         putchar('-');

         x=-x;

     }

     if(x>) write(x/);

     putchar(x%+'');

 }

 void exgcd(ll a,ll b,ll &d,ll &x,ll &y)

 {

     if(!b)

     {

         d=a;

         x=;

         y=;

     }

     else

     {

         exgcd(b,a%b,d,x,y);

         ll t=x;x=y;y=t-a/b*y;

     }

 }

 int main()

 {

     int k;k=read();ll m=,ans=;

     for(int i=;i<=k;i++)a[i]=read();

     for(int i=;i<=k;i++)

     {

         b[i]=read();

         m*=b[i];

     }

     for(int i=;i<=k;i++)

     {

         ll mi=m/b[i],d,x,y;

         exgcd(mi,b[i],d,x,y);

         ans=(ans+a[i]*mi*x)%m;

     }

     printf("%lld",(ans+m)%m);

     return ;

结果居然发现——

然后我赶快回去看了几眼代码，感觉没啥毛病，于是又看了下别人的题解——

喔，我手速加了个快速乘，扫了一遍代码，测了遍样例 样例并没软用，再次Ctrl + c 、 Ctrl + v，按下了提交键，这次肯定没问题的吶~

 #include<bits/stdc++.h>

 #define ll long long

 using namespace std;

 const int N=;

 int n;

 int a[],b[];

 int read()

 {

     int x=,f=;

     char ch=getchar();

     while(ch<''||ch>'')

     {

         if(ch=='-') f=-;

         ch=getchar();

     }

     while(ch>=''&&ch<='')

     {

         x=x*+ch-'';

         ch=getchar();

     }

     return x*f;

 }

 void write(int x)

 {

     if(x<)

     {

         putchar('-');

         x=-x;

     }

     if(x>) write(x/);

     putchar(x%+'');

 }

 void exgcd(ll a,ll b,ll &d,ll &x,ll &y)

 {

     if(!b)

     {

         d=a;

         x=;

         y=;

     }

     else

     {

         exgcd(b,a%b,d,x,y);

         ll t=x;x=y;y=t-a/b*y;

     }

 }

 ll ff(ll a,ll b,ll m)

 {

     ll ans=;

     while(b)

     {

         if(b&)ans=(ans+a)%m;

         a=(a+a)%m;

         b>>=;

     }

     return ans;

 }

 int main()

 {

     int k;k=read();ll m=,ans=;

     for(int i=;i<=k;i++)a[i]=read();

     for(int i=;i<=k;i++)

     {

         b[i]=read();

         m*=b[i];

     }

     for(int i=;i<=k;i++)

     {

         ll mi=m/b[i],d,x,y;

         exgcd(mi,b[i],d,x,y);

         x=(x%b[i]+b[i])%b[i];

         ans=(ans+ff(ff(mi,x,m),a[i],m))%m;

     }

     printf("%lld",(ans+m)%m);

     return ;

 }
 #include<bits/stdc++.h>

 #define ll long long

 using namespace std;

 const int N=;

 int n;

 int a[],b[];

 int read()

 {

     int x=,f=;

     char ch=getchar();

     while(ch<''||ch>'')

     {

         if(ch=='-') f=-;

         ch=getchar();

     }

     while(ch>=''&&ch<='')

     {

         x=x*+ch-'';

         ch=getchar();

     }

     return x*f;

 }

 void write(int x)

 {

     if(x<)

     {

         putchar('-');

         x=-x;

     }

     if(x>) write(x/);

     putchar(x%+'');

 }

 void exgcd(ll a,ll b,ll &d,ll &x,ll &y)

 {

     if(!b)

     {

         d=a;

         x=;

         y=;

     }

     else

     {

         exgcd(b,a%b,d,x,y);

         ll t=x;x=y;y=t-a/b*y;

     }

 }

 ll ff(ll a,ll b,ll m)

 {

     ll ans=;

     while(b)

     {

         if(b&)ans=(ans+a)%m;

         a=(a+a)%m;

         b>>=;

     }

     return ans;

 }

 int main()

 {

     int k;k=read();ll m=,ans=;

     for(int i=;i<=k;i++)a[i]=read();

     for(int i=;i<=k;i++)

     {

         b[i]=read();

         m*=b[i];

     }

     for(int i=;i<=k;i++)

     {

         ll mi=m/b[i],d,x,y;

         exgcd(mi,b[i],d,x,y);

         x=(x%b[i]+b[i])%b[i];

         ans=(ans+ff(ff(mi,x,m),a[i],m))%m;

     }

     printf("%lld",(ans+m)%m);

     return ;

然鹅打脸就是来的这么突然(￣ε(#￣)☆╰╮(￣▽￣///)

我特么第二个点TLE掉是怎么回事？！！

不是，第一遍提交都莫得问题的鸭！我都忍不住要口吐芬芳(｀へ´)

火速赶到题解区翻到了之前没看完的题解——

╮(╯▽╰)╭这毒瘤题，真拿它没办法┑(￣Д ￣)┍

第三次提交，终于满屏全绿（要想生活过得去，做题就得来点绿~

 #include<bits/stdc++.h>

 #define ll long long

 using namespace std;

 const int N=;

 int n;

 int a[],b[];

 int read()

 {

     int x=,f=;

     char ch=getchar();

     while(ch<''||ch>'')

     {

         if(ch=='-') f=-;

         ch=getchar();

     }

     while(ch>=''&&ch<='')

     {

         x=x*+ch-'';

         ch=getchar();

     }

     return x*f;

 }

 void write(int x)

 {

     if(x<)

     {

         putchar('-');

         x=-x;

     }

     if(x>) write(x/);

     putchar(x%+'');

 }

 void exgcd(ll a,ll b,ll &d,ll &x,ll &y)

 {

     if(!b)

     {

         d=a;

         x=;

         y=;

     }

     else

     {

         exgcd(b,a%b,d,x,y);

         ll t=x;x=y;y=t-a/b*y;

     }

 }

 ll ff(ll a,ll b,ll m)

 {

     ll ans=;

     while(b)

     {

         if(b&)ans=(ans+a)%m;

         a=(a+a)%m;

         b>>=;

     }

     return ans;

 }

 int main()

 {

     int k;k=read();ll m=,ans=;

     for(int i=;i<=k;i++)a[i]=read();

     for(int i=;i<=k;i++)

     {

         b[i]=read();

         m*=b[i];

     }

     for(int i=;i<=k;i++)

     {

         ll mi=m/b[i],d,x,y;

         exgcd(mi,b[i],d,x,y);

         x=(x%b[i]+b[i])%b[i];

         ans=(ans+ff(ff(mi,x,m),(a[i]+m)%m,m))%m;

     }

     printf("%lld",(ans+m)%m);

     return ;

 }

//参考：lahlah 的博客

[TJOI2009]猜数字（洛谷 3868）的更多相关文章

  1. P3868 [TJOI2009]猜数字

    [TJOI2009]猜数字 中国剩余定理 求解i=1 to n : x≡a[i] (mod b[i])的同余方程组 设 t= ∏i=1 to n b[i] 我们先求出 i=1 to n : x≡1 ( ...

  2. 洛谷P3868 [TJOI2009]猜数字（中国剩余定理，扩展欧几里德）

    洛谷题目传送门 90分WA第二个点的看过来! 简要介绍一下中国剩余定理 中国剩余定理,就是用来求解这样的问题: 假定以下出现数都是自然数,对于一个线性同余方程组(其中\(\forall i,j\in[ ...

  3. 洛谷 P3868 [TJOI2009]猜数字

    题意简述 给定\(a[1],a[2],\cdots,a[n]\) 和 \(b[1],b[2],\cdots,b[n]\),其中\(b\)中元素两两互素. 求最小的非负整数\(n\),满足对于任意的\( ...

  4. USACO1.4 1.5 搜索剪枝与数字 洛谷OJ P1214 P1215 P1217 P1218

    USACO1.4 题解 Arithmetic Progressions 题意 让你求长为n的由小于2*m*m的双平方数组成的等差数列有几个 双平方数:形如 B=P*P+Q*Q,p,q>0的数 题 ...

  5. [TJOI2009]猜数字

    题目描述 现有两组数字,每组k个,第一组中的数字分别为:a1,a2,...,ak表示,第二组中的数字分别用b1,b2,...,bk表示.其中第二组中的数字是两两互素的.求最小的非负整数n,满足对于任意 ...

  6. CRT【p3868】[TJOI2009]猜数字

    Description 现有两组数字,每组k个,第一组中的数字分别为:a1,a2,...,ak表示,第二组中的数字分别用b1,b2,...,bk表示.其中第二组中的数字是两两互素的.求最小的非负整数n ...

  7. [Luogu3868] [TJOI2009]猜数字

    题目描述 现有两组数字,每组k个,第一组中的数字分别为:a1,a2,...,ak表示,第二组中的数字分别用b1,b2,...,bk表示.其中第二组中的数字是两两互素的.求最小的非负整数n,满足对于任意 ...

  8. [TJOI2009] 猜数字 - 中国剩余定理

    现有两组数字,每组k个,第一组中的数字分别为:a1,a2,...,ak表示,第二组中的数字分别用b1,b2,...,bk表示.其中第二组中的数字是两两互素的.求最小的非负整数n,满足对于任意的i,n ...

  9. Luogu P3868 [TJOI2009]猜数字

    题目链接 \(Click\) \(Here\) 中国剩余定理的板子.小心取模. #include <bits/stdc++.h> using namespace std; const in ...

随机推荐

  1. [转帖]linux基础知识大纲

    linux基础知识大纲 https://blog.csdn.net/CSDN___LYY/article/details/80810403 1.Linux操作系统概述Linux操作系统的发展过程.创始 ...

  2. 轻量级搜索工具【Everything】的设置

    下面是我对 轻量级搜索工具[Everything]的设置 Everything 下载地址:http://www.voidtools.com/downloads/ 1. 排除 Windows => ...

  3. GitBook基本使用

    1.安装 Node.js:https://nodejs.org/en/ 2.安装 GitBook: npm install gitbook-cli -g gitbook -V  #查看gitbook是 ...

  4. java虚拟机JVM基础

    (1)内存模型以及分区,需要详细到每个区放什么. 堆:Java虚拟机管理内存中最大的一块,线程共享区域.所有对象实例和数组都在堆上分配内存空间. 栈:线程私有,每个线程都会创建一个虚拟机栈,生命周期与 ...

  5. 在Nginx服务器上安装SSL证书

    配置nginx 1.下载证书文件 2.在nginx的conf目录中创建目录cert目录,并将证书文件拷贝进去. 3.配置nginx.conf,完整的nginx.conf如下: #user nobody ...

  6. Mybatis一对一，一对多，多对多代码

    一对一 <!-- 关系映射 --> <!-- 1-1:自动映射 --> <select id="oneToOne" resultType=" ...

  7. logger.info占位符的使用

    {}表示占位符,使用方法如下: package org.pine.controller; import javax.annotation.Resource; import org.pine.servi ...

  8. weblogic漏洞

    1 weblogic控制台地址暴露 ²  整改方法: 禁用weblogic控制台.在weblogic域(sguap-domain和base-domain)的config下的config.xml中 &l ...

  9. APS系统帮助寻找企业最优库存

    零库存模式的实施要有深厚的民族文化和企业文化为支点.随着对零库存管理研究的深入,就会发现它不仅仅是一种运营管理技术,更是一种文化.一种哲学. 当这种认同文化体现在企业与企业之间时,就会表现出彼此的认同 ...

  10. MYSQL的修改表结构SQL语句

    更多java学习资料>>> 1.背景 使用sql语句对表结构进行修改 2.案例演示 案例:表结构 CREATE TABLE `login_user` ( `id` ) NOT NUL ...