poj1061 青蛙的约会扩展欧几里得

时间：2015-05-09 17:33:38 阅读：122 评论：0 收藏：0 [点我收藏+]

标签：

//(x + km)%l == (y + kn)%l
//s*l + k(m - n) = y - x
//用扩展欧几里得求出k
//注意一下k需要为最小正数
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std ;
typedef __int64 ll ;
ll exgcd(ll a , ll b , ll &x , ll &y)
{
if(b == 0)
{
x = 1;
y = 0;
return a ;
}
ll r = exgcd(b , a%b , x , y) ;
ll temp = x ;
x = y ;
y = temp - a/b*y ;
return r;
}
int main()
{
ll x_1 , y_1 , m , l , n ;
while(~scanf("%I64d%I64d%I64d%I64d%I64d" , &x_1 , &y_1 , &m ,&n ,&l))
{
ll x , y ;
ll r = exgcd(l , m-n , x , y) ;
if((x_1 - y_1) % r != 0)
puts("Impossible") ;
else
{
y *= ((y_1 - x_1)/r);
while(y > 0){y -= l;}
while(y < 0){y += l;}
printf("%I64d\n" ,y) ;
}
}
}

poj1061 青蛙的约会扩展欧几里得

标签：

原文地址：http://blog.csdn.net/cq_pf/article/details/45602553

踩

(0)

评论一句话评论（0）

分享档案

更多>

2021年07月29日 (22)
2021年07月28日 (40)
2021年07月27日 (32)
2021年07月26日 (79)
2021年07月23日 (29)
2021年07月22日 (30)
2021年07月21日 (42)
2021年07月20日 (16)
2021年07月19日 (90)
2021年07月16日 (35)

周排行

poj1061 青蛙的约会 扩展欧几里得

poj1061 青蛙的约会扩展欧几里得