ZOJ - 3593 One Person Game （扩展欧几里得技巧）

最新推荐文章于 2019-10-13 00:24:09 发布

_wjunjie

最新推荐文章于 2019-10-13 00:24:09 发布

阅读量302

点赞数

CC 4.0 BY-SA版权

分类专栏：数学

本文链接：https://blog.youkuaiyun.com/qq_41661919/article/details/81456627

数学专栏收录该内容

29 篇文章

订阅专栏

本文通过一个具体的编程实例，深入介绍了扩展欧几里德算法的应用，特别是如何解决类似于poj2142-TheBalance这样的问题。文章详细展示了算法的实现过程，包括求最大公约数及解线性方程组等内容。

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

类似题目：poj2142-The Balance(扩展欧几里德算法)

#include<iostream>
#include<cmath>
#define ll long long
#define INF 0x7fffffff
using namespace std;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
    if(!b)
    {
        x=1;
        y=0;
        return a;
    }
    ll gcd=exgcd(b,a%b,x,y),t;
    t=x;
    x=y;
    y=t-a/b*y;
    return gcd;
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        ll A,B,a,b;
        ll x,y;
        cin>>A>>B>>a>>b;
        ll C=B-A;
        ll gcd=exgcd(a,b,x,y);
        if(C%gcd==0)
        {
            x*=C/gcd;
            y*=C/gcd;//通解形式
            ll mid=(y-x)*gcd/(a+b);
            ll minn=-1,ans=(ll)(INF)*(ll)(INF);
            for(ll i=mid-1;i<=mid+1;++i)
            {
                if((x+b/gcd*i)*(y-a/gcd*i)>=0)
                {
                    minn=max(fabs(x+b/gcd*i),fabs(y-a/gcd*i));
                }
                else
                {
                    minn=fabs(x+b/gcd*i)+fabs(y-a/gcd*i);
                }
                ans=min(ans,minn);

            }
            cout<<ans<<endl;

        }
        else cout<<"-1"<<endl;

    }
    return 0;
}