Problem
http://poj.org/problem?id=1061
Solution
题外话:右边是西啊。。。
本题是扩欧的必修题,就打打看了看。
我们设跳了k次碰面。
可得
(
x
+
k
∗
m
)
−
(
y
+
k
∗
n
)
=
L
∗
p
(
p
为
整
数
)
(x+k*m)-(y+k*n)=L*p(p为整数)
(x+k∗m)−(y+k∗n)=L∗p(p为整数)
我们将其化简,可得:
(
n
−
m
)
∗
k
+
L
∗
p
=
x
−
y
(n-m)*k+L*p=x-y
(n−m)∗k+L∗p=x−y
其中,
a
=
n
−
m
,
b
=
L
,
c
=
x
−
y
a=n-m,b=L,c=x-y
a=n−m,b=L,c=x−y
哈哈,这不就是扩欧嘛!
a
∗
x
1
+
b
∗
y
1
=
c
a*x1+b*y1=c
a∗x1+b∗y1=c
咚咚打完,输出
x
1
∗
(
c
/
g
c
d
)
x1*(c/gcd)
x1∗(c/gcd)
结果WA了。
后来发现该答案好像不是最优解。
我们可以发现:
如果答案
k
>
L
k>L
k>L的话,我们便可以将
k
−
L
k-L
k−L,然后将
p
+
(
n
−
m
)
p+(n-m)
p+(n−m)。
如果答案
k
<
0
k<0
k<0的话,反之。
这样,求得的答案才会是最优解了。
Code
#include<cstdio>
#define ll long long
using namespace std;
ll x,y,m,n,L,x1,y1;
ll a,b,c,gcd;
void exgcd(ll x,ll y,ll &x1,ll &y1)
{
if (!y) {gcd=x,x1=1,y1=0; return;}
exgcd(y,x%y,y1,x1),y1-=x1*(x/y);
}
int main()
{
freopen("frog.in","r",stdin);
freopen("frog.out","w",stdout);
scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&L);
a=n-m,b=L,c=x-y;
exgcd(a,b,x1,y1);
if (c%gcd) {puts("Impossible"); return 0;}
x1=x1*c/gcd,y1=y1*c/gcd;
printf("%lld\n",(x1%b+b)%b);
return 0;
}