扩展欧几里得算法——例题4: 最大公约数问题2

本文介绍了解最大公约数方程AX+BY=C的一种方法,给出了一组整数解X和Y的求解思路及代码实现。适用于1<=A,B<=1000且1<=C<=10^9的情况。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

最大公约数问题2

题目描述

输入正整数A,B,C,求一组X, Y,使得方程:AX+BY=C,保证有解。

输出任何一组解即可。

输入

第1行:1个整数T,表示测试数据的组数 (1 <= T <= 100)

接下来T行,每行2个整数,表示A和B和C (1 <= A, B <= 1000, 1 <= C <= 10^9)

输出

输出T行,每行2个整数,表示X和Y

样例输入

1
14 8 4

样例输出

-2 4

解析

这和最大公约数问题1有相同之处。我们先求出 AX + BY = GCD(A,B)的整数解X,Y。

∵ AX + BY = GCD(A,B)

∴ A(CX) + B(CY) = GCD(A,B) * C

∴ A(XC / GCD(A,B)) + B(YC / GCD(A,B)) = C

若有整数解的话,就必须满足 C / GCD(A,B) 也是整数,即 C mod GCD(A,B) = 0,题目保证有解,就不管了,直接除。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
using namespace std;
int read() {
    int f = 1,x = 0;char s = getchar();
    while(s < '0' || s > '9') {if(s == '-')f = -1;s = getchar();}
    while(s >= '0' && s <= '9') {x = x * 10 + s - '0';s = getchar();}
    return x * f;
}
bool f[2005];
int p[2005];
int n,m,i,j,k,s,o,num,num1;
int exgcd(int a,int b,int &x,int &y) {
    if(!b) {
        x = 1,y = 0;
        return a;
    }
    else {
        int r = exgcd(b,a % b,y,x);
        y -= x * (a / b);
        return r;
    }
}
void answer(int a,int b,int c,int &x,int &y) {
    int cu = c / exgcd(a,b,x,y);
    x *= cu;y *= cu;
}
int main() {
    n = read();
    while(n --) {
        s = read();o = read();k = read();
        answer(s,o,k,i,j);
        printf("%d %d\n",i,j);
    }
    return 0;
}
 

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值