POJ 2142 The Balance .

最新推荐文章于 2020-07-17 22:27:55 发布

原创最新推荐文章于 2020-07-17 22:27:55 发布 · 274 阅读

0 ·

CC 4.0 BY-SA版权

POJ+百练同时被 2 个专栏收录

120 篇文章

订阅专栏

数论：GCD

5 篇文章

订阅专栏

本文介绍了解决POJ 2142问题的方法，通过分析两种情况得出方程ax + by = d，并运用扩展欧几里得算法求解。详细解释了如何找到x与y的最小非负解，最终确定使|x| + |y|值最小的解。

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

题目地址：http://poj.org/problem?id=2142

有两种可能的情况：

a与d放在一起，即 ax+d=by

b与d放在一起，即 ax=by+d

因为x，y是任取的，所以解下面一个方程

ax+by=d

利用扩展欧几里得定理，可得两组解，

1. x的最小非负解，及其对应的y

2. y的最小非负解，及其对应的x

因为求的是|x|+|y|的最小值，比较一下就好了

#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
int gcdEx(int a,int b,int& x,int& y)
{   //求ax+by=gcd(a,b) 的整数解 返回gcd(a,b); 
	if(b==0) {
		x=1; y=0; return a;
	}
	int x1,y1;
	int gcd=gcdEx(b,a%b,x1,y1);
	x=y1;
	y=x1-a/b*y1;
	return gcd;
}
int main()
{
	int a,b,d; int x,y;
	while(cin>>a>>b>>d)
	{
		if(a==0&&b==0&&d==0) break;
		int gd=gcdEx(a,b,x,y);
		int k=d/gd,s1=b/gd,s2=a/gd;
		int x1=(k*x%s1+s1)%s1,y1=fabs(d-a*x1)/b; 
		int y2=(k*y%s2+s2)%s2,x2=fabs(d-b*y2)/a; 
		if(x1+y1<x2+y2)	cout<<x1<<' '<<y1<<endl;
		else cout<<x2<<' '<<y2<<endl;
	}
	return 0;
}