这是一道标准的孙子定理的题,题意浅显,思路明确

然后我就交了整整16遍啊,欺负人啊,题解暴力就能过,我就TLE

。。悲惨的提交记录







下面是题面

题目描述

自从曹冲搞定了大象以后,曹操就开始捉摸让儿子干些事业,于是派他到中原养猪场养猪,可是曹冲满不高兴,于是在工作中马马虎虎,有一次曹操想知道母猪的数量,于是曹冲想狠狠耍曹操一把。举个例子,假如有16头母猪,如果建了3个猪圈,剩下1头猪就没有地方安家了。如果建造了5个猪圈,但是仍然有1头猪没有地方去,然后如果建造了7个猪圈,还有2头没有地方去。你作为曹总的私人秘书理所当然要将准确的猪数报给曹总,你该怎么办?

输入输出格式

输入格式:

第一行包含一个整数n (n <= 10) – 建立猪圈的次数,解下来n行,每行两个整数ai, bi( bi <= ai <= 1000), 表示建立了ai个猪圈,有bi头猪没有去处。你可以假定ai,aj互质.

输出格式:

输出包含一个正整数,即为曹冲至少养母猪的数目。

输入输出样例

输入样例

3

3 1

5 1

7 2

输出样例

16

题目已经描述的非常清晰了

先上90分代码,根据代码进行讲解

#include<iostream>
#include<cstdio>
#define ll __int128 //实际上开long long也能过,我后来试过了,当时因为怀疑long long不够,多此一举
#define re register //可耻的register也并没有使我A掉最后一个点
using namespace std;
const int maxn=15;
int n,a[maxn],b[maxn];
ll tot=1,ans; //tot即所有数乘积,使用时除以当前数即可,避免多次计算浪费时间
void read(int& a){ //毫无意义的快读
a=0;re char p;
p=getchar();
while(p>'9'||p<'0')
p=getchar();
while(p>='0'&&p<='9'){
a=(a<<3)+(a<<1)+(p^48);
p=getchar();
}
return;
}
void print(ll a){ //为int_128写的输出,一般用不到
if(a==0){
printf("0");
return;
}
string pp="";
while(a){
pp+=a%10+'0';
a/=10;
}
for(int i=pp.size()-1;i>=0;i--)printf("%c",pp[i]);
}
int main(){
read(n);
for(re int i=1;i<=n;++i){
read(a[i]),read(b[i]);
tot*=a[i];
}
for(re int i=1;i<=n;++i){
re ll p=tot/a[i],sum=p;
while(sum%a[i]!=1) //求其余数的公倍数,使其取模当前数余1
sum+=p;
ans=(ans+b[i]*sum)%tot;
}
print(ans);
return 0;
}

因为有暴力那一块,无论怎么优化都A不了(如果有哪位大佬能优化A掉,请私信我,情愿被打脸)

于是乎,将暴力推改为直接扩欧解同余方程,一次AC

下放满分代码

#include<iostream>
#include<cstdio>
#define ll __int128
#define re register
using namespace std;
const int maxn=15;
int n,a[maxn],b[maxn];
ll tot=1,ans;
void read(int& a){
a=0;re char p;
p=getchar();
while(p>'9'||p<'0')
p=getchar();
while(p>='0'&&p<='9'){
a=(a<<3)+(a<<1)+(p^48);
p=getchar();
}
return;
}
void print(ll a){
if(a==0){
printf("0");
return;
}
string pp="";
while(a){
pp+=a%10+'0';
a/=10;
}
for(int i=pp.size()-1;i>=0;i--)printf("%c",pp[i]);
}
ll e_gcd(ll a,ll b,ll &x,ll &y){ //就是这里,解同余方程
if(b==0){
x=1;y=0;
return a;
}
ll ans=e_gcd(b,a%b,x,y);
ll temp=x;
x=y;
y=temp-a/b*y;
return ans;
}
ll cal(ll a,ll b){
ll x,y;
ll gcd=e_gcd(a,b,x,y);
if(1%gcd!=0)return -1;
x*=1/gcd;
return (x%b+b)%b;
}
int main(){
read(n);
for(re int i=1;i<=n;++i){
read(a[i]),read(b[i]);
tot*=a[i];
}
for(re int i=1;i<=n;++i){
re ll p=tot/a[i];
ans=(ans+b[i]*p*cal(p,a[i]))%tot;
}
print(ans);
return 0;
}

最新文章

  1. tensorflow学习笔记四:mnist实例--用简单的神经网络来训练和测试
  2. json格式的优点
  3. linux 目录
  4. cf486A Calculating Function
  5. java中字符串的操作
  6. jquery +ajax 上传加预览
  7. POST调用WCF方法-项目实践
  8. leetcode-83.删除排序链表中的重复元素
  9. ASP.NET MVC5 + EF6 + LayUI实战教程,通用后台管理系统框架(3)
  10. javascript 面向对象之路.1 - 小蜜蜂
  11. SD-WAN供应商列表
  12. 在ubuntu下运行python脚本
  13. Mac安装WineHQ
  14. PgSQL &middot; 源码分析&middot; pg_dump分析
  15. brew - 安装gradle
  16. 1.node.js下载
  17. BMP格式详解
  18. 4-fiddler抓包中文乱码:
  19. 前端页面 script 事件总结
  20. 【互动问答分享】第11期决胜云计算大数据时代Spark亚太研究院公益大讲堂

热门文章

  1. CF 1064 D. Labyrinth
  2. L013-linux基础正则表达式手把手实战讲解小节
  3. 【JUC源码解析】Semaphore
  4. ubuntu 16.04 安装php 5 6等版本
  5. AtCoder Grand Contest 026 D - Histogram Coloring
  6. EXE模块说明
  7. 从零开始的Python学习Episode 10——函数
  8. tendermint 跟tikv结合
  9. How to pass an Amazon account review
  10. sqlserver-查阻塞