博客
关于我
HDU 2669 Romantic(扩展欧几里得算法)
阅读量:719 次
发布时间:2019-03-21

本文共 1394 字,大约阅读时间需要 4 分钟。

为了解决这个问题,我们需要找到满足方程 ax + by = 1 的整数解,其中 x 是非负的,并且尽可能小。只有当 a 和 b 互质时,这个方程才有解。我们可以通过扩展欧几里得算法来找到一组解,并对其进行调整以满足要求。

方法思路

  • 检查互质性:首先检查 a 和 b 是否互质,即它们的最大公约数是否为 1。只有当它们互质时,方程才有解。
  • 扩展欧几里得算法:使用扩展欧几里得算法找到一个初始解 (x0, y0)。这个解可能会有负数的 x。
  • 调整解:将初始解进行调整,使得 x 变为正数。通过调整,我们可以找到所有可能的解中的 x 最小的非负的解。
  • 解决代码

    #include 
    #pragma warning(disable:4996)int a, b, x, y;int gcd(int a, int b) { return b == 0 ? a : gcd(b, a % b);}int extgcd(int a, int b, int x, int y) { int d = a; if (b != 0) { d = extgcd(b, a % b, x, y); y -= (a / b) * x; } else { x = 1; y = 0; } return d;}int main() { while (scanf("%d %d", &a, &b) != EOF) { if (extgcd(a, b, x, y) > 1) { puts("sorry"); continue; } if (y < 0) { y = -y; a = -a; } if (x < 0) { x = -x; b = -b; } if (extgcd(a, b, x, y, 1) != 1) { puts("sorry"); continue; } while (x >= 0) { x -= b; y += a; } while (x < 0) { x += b; y -= a; } printf("%d %d\n", x, y); } return 0;}

    代码解释

  • 输入处理:使用 scanf 读取输入,直到 EOF。
  • 检查互质性:通过扩展欧几里得算法检查 a 和 b 的最大公约数,如果大于 1,则输出 "sorry"。
  • 初始解调整:对于可能的负数解,通过调整 x 和 y 以确保它们是正整数。
  • 调整解:利用扩展欧几里得算法的结果,进一步调整 x 和 y,使得 x 为最小的非负解。
  • 输出结果:打印满足条件的 x 和 y 的值。
  • 通过这种方法,我们能够高效地找到满足条件的解,并确保输出的 x 是最小的非负整数。

    转载地址:http://lotez.baihongyu.com/

    你可能感兴趣的文章
    OSG学习:几何体的操作(二)——交互事件、Delaunay三角网绘制
    查看>>
    OSG学习:几何对象的绘制(一)——四边形
    查看>>
    OSG学习:几何对象的绘制(三)——几何元素的存储和几何体的绘制方法
    查看>>
    OSG学习:几何对象的绘制(二)——简易房屋
    查看>>
    OSG学习:几何对象的绘制(四)——几何体的更新回调:旋转的线
    查看>>
    OSG学习:场景图形管理(一)——视图与相机
    查看>>
    OSG学习:场景图形管理(三)——多视图相机渲染
    查看>>
    OSG学习:场景图形管理(二)——单窗口多相机渲染
    查看>>
    OSG学习:场景图形管理(四)——多视图多窗口渲染
    查看>>
    OSG学习:新建C++/CLI工程并读取模型(C++/CLI)——根据OSG官方示例代码初步理解其方法
    查看>>
    Sql 随机更新一条数据返回更新数据的ID编号
    查看>>
    OSG学习:空间变换节点和开关节点示例
    查看>>
    OSG学习:纹理映射(一)——多重纹理映射
    查看>>
    OSG学习:纹理映射(七)——聚光灯
    查看>>
    OSG学习:纹理映射(三)——立方图纹理映射
    查看>>
    OSG学习:纹理映射(二)——一维/二维/简单立方图纹理映射
    查看>>
    OSG学习:纹理映射(五)——计算纹理坐标
    查看>>
    OSG学习:纹理映射(六)——灯光
    查看>>
    OSG学习:纹理映射(四)——三维纹理映射
    查看>>
    OSG:从源码看Viewer::run() 一
    查看>>