实数即使用double定义也会损失一定精度,在精度要求很高时,就需要用到分数
0.最大公约数&最小公倍数
#define LL long long
LL gcd(LL a,LL b) //最大公约数
{
if(b==0)
return a;
return gcd(b,a%b);
}
LL lcm(LL a,LL b) //最小公倍数
{
return a*b/gcd(a,b);
}
1.分数存储
struct Fraction //分数
{
LL a,b;
void update() //约分
{
LL k=gcd(a,b);
a/=k;
b/=k;
return;
}
};
2.“+”
Fraction operator+(Fraction x,Fraction y)
{
Fraction z;
z.a=x.a*y.b+x.b*y.a;
z.b=x.b*y.b;
z.update();
return z;
}
3.“ - ”
Fraction operator-(Fraction x,Fraction y)
{
Fraction z;
z.a=x.a*y.b-x.b*y.a;
z.b=x.b*y.b;
z.update();
return z;
}
4.“ * ”
Fraction operator*(Fraction x,Fraction y)
{
Fraction z;
z.a=x.a*y.a;
z.b=x.b*y.b;
z.update();
return z;
}
5." / "
Fraction operator/(Fraction x,Fraction y)
{
Fraction z;
z.a=x.a*y.b;
z.b=x.b*y.a;
z.update();
return z;
}
6.大小判断
bool operator<(Fraction x,Fraction y)
{
return x.a*y.b<y.a*x.b;
}
bool operator<=(Fraction x,Fraction y)
{
return x.a*y.b<=y.a*x.b;
}
bool operator>(Fraction x,Fraction y)
{
return x.a*y.b>y.a*x.b;
}
bool operator>=(Fraction x,Fraction y)
{
return x.a*y.b>=y.a*x.b;
}
bool operator==(Fraction x,Fraction y)
{
return x.a*y.b==y.a*x.b;
}
Code
//比例简化
#include<iostream>
#include<cstdio>
#define LL long long
#define INF 0x3f3f3f3f3f
using namespace std;
const int L=100;
LL gcd(LL a,LL b) //最大公约数
{
if(b==0)
return a;
return gcd(b,a%b);
}
LL lcm(LL a,LL b) //最小公倍数
{
return a*b/gcd(a,b);
}
struct Fraction //分数
{
LL a,b;
void update()
{
LL k=gcd(a,b);
a/=k;
b/=k;
return;
}
};
Fraction operator+(Fraction x,Fraction y)
{
Fraction z;
z.a=x.a*y.b+x.b*y.a;
z.b=x.b*y.b;
z.update();
return z;
}
Fraction operator-(Fraction x,Fraction y)
{
Fraction z;
z.a=x.a*y.b-x.b*y.a;
z.b=x.b*y.b;
z.update();
return z;
}
Fraction operator*(Fraction x,Fraction y)
{
Fraction z;
z.a=x.a*y.a;
z.b=x.b*y.b;
z.update();
return z;
}
Fraction operator/(Fraction x,Fraction y)
{
Fraction z;
z.a=x.a*y.b;
z.b=x.b*y.a;
z.update();
return z;
}
bool operator<(Fraction x,Fraction y)
{
return x.a*y.b<y.a*x.b;
}
bool operator<=(Fraction x,Fraction y)
{
return x.a*y.b<=y.a*x.b;
}
bool operator>(Fraction x,Fraction y)
{
return x.a*y.b>y.a*x.b;
}
bool operator>=(Fraction x,Fraction y)
{
return x.a*y.b>=y.a*x.b;
}
bool operator==(Fraction x,Fraction y)
{
return x.a*y.b==y.a*x.b;
}
Fraction t;
Fraction x,y;
Fraction ans;
Fraction Min,tmp;
int l;
int main()
{
// freopen("1.in","r",stdin);
scanf("%lld%lld",&t.a,&t.b);
t.update();
scanf("%d",&l);
Min.a=INF;
Min.b=1;
for(int i=1;i<=l;i++) {
x.a=i;
for(int j=1;j<=l;j++) {
x.b=j;
if(x<t) continue;
tmp=x-t;
if(tmp<Min)
Min=tmp,ans=x;
}
}
ans.update();
printf("%lld %lld\n",ans.a,ans.b);
return 0;
}
来源:CSDN
作者:GBLoi
链接:https://blog.csdn.net/qq_45646484/article/details/104195893