这是一道标准的孙子定理的题,题意浅显,思路明确
然后我就交了整整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;
}