首页 > 代码库 > POJ 1601 青蛙的约会
POJ 1601 青蛙的约会
题意是两只青蛙,一开始分别在x,y点,他们分别每单位时间可以移动m,n个点位时间,路径是一个单位长度为l的圆环
通过题意很容易得到一个方程:
(n-m)*k+t*l=x-y;
其中k和t是未知数,(n-m)以及l我们可以通过输入确定
看到这个方程我们联想到ax+by=d;这个二元一次方程。d=gcd(a,b);
我们不妨另a=n-m b=l
则原式就为a*k+b*t=(x-y);。。。。1
已知1有整数解的条件为(d|(x-y)))
所以如果(x-y)%d!=0 那么必然无整数解,则输出impossible
另c= (x-y)
我们不妨给a/=d b/=d c/=d
利用欧几里得扩展 excgcd(a,b,x,y)
计算出x; 有欧几里得扩展的意义可以知道 x *=c;
然后计算最小值即可
#include <iostream> #include <cstring> #include <cstdio> using namespace std; long long d;//代表gcd(a,l) long long gcd(long long a,long long b) { if(b==0) { return a; } else{ gcd(b,a%b); } } void extgcd(long long a,long long b,long long &x,long long &y) { if(b==0) { x=1; y=0; } else{ extgcd(b,a%b,y,x); y -= (a/b)*x; } } long long x,y,m,n,l; int main() { while(cin>>x>>y>>m>>n>>l) { long long a = n-m; long long b = x-y; long long k,y; d = gcd(a,l); if(b%d!=0) { printf("Impossible\n"); continue; } else{ a/=d; l/=d; b/=d; extgcd(a,l,k,y); k *= b; long long t =k%l;//????? while(t<0) { t+=l; } cout<<t<<endl; } } return 0; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。