HDU 2669 Romantic(扩展欧几里德裸题)

裸的扩展欧几里德解不定方程,话不多说,看代码中间
#include<iostream>
#include<cstdio>

using namespace std;
int x, y, c, r;
int exgcd(int a, int b, int &x, int &y)//知道扩展欧几里德此处不难写,裸模版;
{
    if(b == 0)
    {
        x = 1;
        y = 0;
        return a;
    }
    int r = exgcd(b, a%b, x, y);
    int t = x;
    x = y;
    y = t - a/b*y;
    return r;
}
bool judge_exgcd(int a, int b, int c, int &x, int &y)
{
    int r = exgcd(a, b, x, y);
    if(c % r != 0)//即 ax + by = c; C mod GCD(a ,b)方程才有整数解,否则没有整数解,返回false表示无解,输出sorry;
    {
        return false;
    }
    int k = c / r; //用来求ax + by = c; 的特解,只不过c = 1; 进行到此处a,b一定互质,也就是r = 1; 即此及一下两行均可省略。 
    x *= k;
    y *= k;
    return true;
}
int main()
{
    int a, b;
    while(scanf("%d%d",&a, &b) != EOF)
    {
        c = 1;
        bool flag = judge_exgcd(a, b, c, x, y);
        if(!flag)
        {
            printf("sorry\n");
            continue;
        }
        while(x < 0)//x有小于零的情况, 别问我为什么,就像妈妈一定有孩子一样,一定有,非要找原因的话,到16, 17行找原因去,把过程中的数字给打印出来看看就知道了
        {
            x += b*10;
            y -= a*10;
        }
        int k = x/b;
        x = x%b;
        y += a*k;
        printf("%d %d\n",x, y);
    }
    return 0;
}

的注释即可:
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值