CodeForces 7C Line Exgcd(数论)

C. Line
time limit per test
1 second
memory limit per test
256 megabytes
input
standard input
output
standard output

A line on the plane is described by an equation Ax + By + C = 0. You are to find any point on this line, whose coordinates are integer numbers from  - 5·1018 to 5·1018 inclusive, or to find out that such points do not exist.

Input

The first line contains three integers AB and C ( - 2·109 ≤ A, B, C ≤ 2·109) — corresponding coefficients of the line equation. It is guaranteed that A2 + B2 > 0.

Output

If the required point exists, output its coordinates, otherwise output -1.

Examples
input
2 5 3
output
6 -3

传送门:http://codeforces.com/problemset/problem/7/C


题意&&思路:找出方程一组通解,扩展欧几里得算法就是解方程(欧几里得算法详解


先附上扩展欧几里得算法模板:

ll gcd(ll a,ll b,ll& x,ll& y)
{
	if(!b)
	{
		x=1;
		y=0;
		return a;
	}	
		ll ans=gcd(b,a%b,x,y);
		int temp=x;
		x=y;
		y=temp-a/b*y;
		return ans;
}

既可以求出最大公约数,还可以顺带求解出使得:a*x + b*y = gcd (a,b) 的通解 x 和 y

算法的返回值为 gcd(a,b) 求出的 x,y 为负

代入此题时因为解为gcd(a,b),因此最后要乘上c/gcd(a,b)

同时拓展一下:求乘法逆元时,等价于求 ax + by = 1 的解,即条件为 gcd(a,b)==1 ;


AC代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
#define INF 5*1e18
using namespace std;
ll gcd(ll a,ll b,ll& x,ll& y)
{
	if(!b)
	{
		x=1;
		y=0;
		return a;
	}	
		ll ans=gcd(b,a%b,x,y);
		int temp=x;
		x=y;
		y=temp-a/b*y;
		return ans;
}
int main()
{
	ll a,b,c,d,x,y;
	while(scanf("%lld%lld%lld",&a,&b,&c)!=EOF)
	{
		d=gcd(a,b,x,y);
		if(c%d!=0)
			printf("-1\n");
		else
			printf("%lld %lld\n",-x*(c/d),-y*(c/d));
	}
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值