拓展欧几里得解方程ax+by=gcd(a,b)

最新推荐文章于 2024-08-05 21:08:45 发布

Triose

最新推荐文章于 2024-08-05 21:08:45 发布

阅读量2.1k

点赞数

分类专栏：解题报告文章标签：解题报告

本文链接：https://blog.csdn.net/Triose/article/details/43953507

版权

本文介绍了如何求解方程ax+by=gcd(a,b)的初始解，提供了具体的输入输出示例，并分享了一段高效的代码实现。文章还提及了后续看到的更优解法，适合对数论和算法感兴趣的读者。" 91879510,8493834,配置STP协议实验与拓扑详解,"['网络协议', '交换机配置', 'STP', '网络实验']

摘要由CSDN通过智能技术生成

好吧，过年在家还是看了点点书，还自己给自己出了两道坑爹题。

这道题目一开始想了半天敲出来了，但是只对一组数据得出的答案是对的。。。。后来又写了一遍就对了。。。。

描述：

求出方程：ax+by=gcd(a,b)的初始解,输入a,b（int），输出x,y（空格隔开）;

输入：

60 22

22 60

12453 2467

输出：

-4 11

11 -4

304 -1613

代码：

#include<stdio.h>
void swap(__int64 &a,__int64 &b)
{
     a^=b;
     b^=a;
     a^=b;
}
__int64 gcd(__int64 a,__int64 b)
{
     return b?gcd(b,a%b):a;
}
void func(__int64 a,__int64 b)
{
     int f;
     f=0;
     if(a<b)
    {
        swap(a,b);
        f=1;
     }
     __int64 x1=1,x2,x,y1=1,y2,y,p,q,r1,r2,r;
     r1=a%b;r2=b%r1;
    q=a/b;p=b/r1;r=r1%r2;
     x1=1;y1=-q;
     x2=-p;y2=1+p*q;
 //    printf("1\n x1:%d y1:%d\n x2:%d y2:%d\nq:%d\n",x1,