One Person Game ZOJ (pintia.cn)
分析:
-
求解 a x + b y = ∣ B − A ∣ ax+by=|B-A| ax+by=∣B−A∣ ,
-
设 x 0 , y 0 x_0 ,y_0 x0,y0 是解,则通解为 x 0 + k ∗ b / g c d 和 y 0 − k ∗ a / g c d x_0+k*b/gcd\ 和\ y_0-k*a/gcd x0+k∗b/gcd 和 y0−k∗a/gcd , k ∈ N k\in N k∈N
-
根据题意,要求的是 ∣ x ∣ + ∣ y ∣ |x|+|y| ∣x∣+∣y∣ , 且 x x x 和 y y y 同号时,可以共有一部分(可以 a + b = c a+b=c a+b=c ),此时 a n s = m a x ( ∣ x ∣ , ∣ y ∣ ) ans=max(|x|,|y|) ans=max(∣x∣,∣y∣)
要使 a n s ans ans 最小,即公共部分最大,最优解是 x = y x=y x=y 的情况
-
若 x = y x=y x=y ,即 x 0 + k ∗ b / g c d = y 0 − k ∗ a / g c d x_0+k*b/gcd=y_0-k*a/gcd x0+k∗b/gcd=y0−k∗a/gcd
故, k = ( y 0 − x 0 ) / ( a / g c d + b / g c d ) k=(y0-x0)/(a/gcd+b/gcd) k=(y0−x0)/(a/gcd+b/gcd) , 再偏移一下,遍历求出答案
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+5;
void exgcd(ll a,ll b,ll &x,ll &y,ll &d)
{
if(!b) { x=1; y=0; d=a; return; }
exgcd(b,a%b,x,y,d);
int t=y;
y=x-a/b*y; x=t;
}
void solve(ll a,ll b,ll c)
{
ll x,y,d;
exgcd(a,b,x,y,d);
if(c%d!=0)
{
cout<<"-1"<<endl;
return ;
}
x*=c/d; y*=c/d;
a/=d; b/=d;
ll mid=(y-x)/(a+b), ans=1e18;
for(ll i=mid-1;i<=mid+1;i++)
{
ll xn=x+i*b, yn=y-i*a;
if(xn*yn>0) // 同号
{
ans=min(ans,max(abs(xn),abs(yn)));
}
else ans=min(ans,abs(xn)+abs(yn)); // 异号
}
cout<<ans<<endl;
}
signed main()
{
ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
int T;
cin>>T;
while(T--)
{
ll a,b,an,bn;
cin>>a>>b>>an>>bn;
ll c=abs(a-b);
solve(an,bn,c);
}
return 0;
}