/**
[exgcd]hdu 4180 给出分数(a/b)求另一个分数(c/d){d < b),且满足fabs(a/b - c/d)最小。
ab互质时
最小时满足 bc+1 = ad || bc = ad+1,正符合扩展欧几里得。
求出两个d后比较。
因为下载了数据,充分见证了精度损失的严重性,比较两个实数的距离时一直WA(用计算器发现小数点后20位才不同)
比较两个d的大小即可
*/
#include <stdio.h>
#define EPS 1e-20
__int64 exgcd(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
if(b == 0LL)
{
x = 1LL;
y = 0LL;
return a;
}
__int64 d = exgcd(b,a%b,x,y);
__int64 t = x;
x = y;
y = t - a / b * y;
return d;
}
int main()
{
__int64 t,a,b,x,y,c1,d1,c2,d2,d;
scanf("%I64d",&t);
while(t--)
{
scanf("%I64d/%I64d",&a,&b);
x =
[exgcd] hdu 4180
最新推荐文章于 2025-04-17 23:24:36 发布