码迷,mamicode.com
首页 > 编程语言 > 详细

POJ 1061 青蛙绕地球约会-数论-(解一元一次同余方程+扩展欧几里得算法)

时间:2015-08-14 06:30:49      阅读:180      评论:0      收藏:0      [点我收藏+]

标签:poj

题意:两只青蛙同向跳,起点是x,y,每次分别跳m,n米,地球周长是L,求最少跳几次相遇。

分析:

把式子写好就发现是一个一元一次同余方程。用扩展欧几里得算法来求。这题很基本得会。

代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<string>
#include<queue>
#define INF 1000000007
using namespace std;
long long x,y,m,n,l;
long long a,b,c,d;
long long p,q,ans;
long long gcd(long long a,long long b)
{
	if(b==0) return a;
	else return gcd(b,a%b);
}
void exgcd(long long a,long long b,long long &p,long long &q)
{
	if(b==0){
		p=1;q=0;return;
	}
	else{
//		exgcd(b,a%b,p,q);  //注释掉的部分也是可行的
//		long long tmp=p;
//		p=q;
//		q=tmp-a/b*q;
		exgcd(b,a%b,q,p);
		q-=a/b*p;
	}
}
int main()
{
	while(cin>>x>>y>>m>>n>>l){
		a=m-n;b=l;
		c=y-x;d=gcd(a,b);
		if(c%d){
			cout<<"Impossible"<<endl;
		}
		else{
			a/=d;b/=d;c/=d;
			exgcd(a,b,p,q);//cout<<p<<" "<<q<<endl;
			p*=c;
			ans=p%l;
			if(ans<0) ans+=l;
			cout<<ans<<endl;
		}
	}
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

POJ 1061 青蛙绕地球约会-数论-(解一元一次同余方程+扩展欧几里得算法)

标签:poj

原文地址:http://blog.csdn.net/ac_0_summer/article/details/47638239

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