码迷,mamicode.com
首页 > 其他好文 > 详细

【扩展欧几里得】BZOJ1477-青蛙的约会

时间:2016-11-13 22:05:59      阅读:152      评论:0      收藏:0      [点我收藏+]

标签:using   long   最大   namespace   gcd   can   坐标   坐标轴   include   

一直在WA,后来我发现我把东西看反了……

【题目大意】

给出一个长度为L的环状坐标轴,两个点开始时位于(X,0)、(Y,0)。每次两点分别往右边移动m和n,问能否相遇?

【思路】

由题意,可得:

X+mt=Y+nt(mod L)

(X+mt)-(Y+nt)=L*k

(n-m)t+L*k=X-Y。

可以用扩展欧几里得来做。具体来说,显然要满足n-m和L的最大公约数(记为d)要整除X-Y,否则无解。这个可以在扩展欧几里得中求出。

式子可以化简为:[(n-m)/d]*t+(L/d)*k=(X-Y)/d。因此通解为 t‘=t0+(L/d)*n(n∈Z)

由于扩展欧几里得杰出的结果其实是(n-m)*t+L*k=d的解,因此要乘以(X-Y)/d才是答案。然后往上加得到第一个大于0的解即可。

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 using namespace std;
 6 typedef long long ll;
 7 
 8 ll extgcd(ll a,ll b,ll &x,ll &y)
 9 {
10     ll ret;
11     if (b==0)
12     {
13         x=1,y=0;
14         return a;
15     }
16     ret=extgcd(b,a%b,x,y);
17     ll tmp=x;
18     x=y;
19     y=tmp-(a/b)*y;
20     return ret;
21 }
22 
23 int main()
24 {
25     ll X,Y,m,n,L,x,y;
26     scanf("%lld%lld%lld%lld%lld",&X,&Y,&m,&n,&L);
27     ll d=extgcd(n-m,L,x,y);
28     if ((X-Y)%d!=0) puts("Impossible");else
29     { 
30         ll delta=abs(L/d);
31         x=(x*(X-Y)/d+delta)%delta;
32         while (x<0) x=(x+delta)%delta;
33         printf("%lld",x);
34 
35     }
36     return 0;
37 } 

 

【扩展欧几里得】BZOJ1477-青蛙的约会

标签:using   long   最大   namespace   gcd   can   坐标   坐标轴   include   

原文地址:http://www.cnblogs.com/iiyiyi/p/6059950.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!