poj 1061カエルのデート(拡張ユークリッドアルゴリズム)

1319 ワード

【テーマ】
poj 1061カエルのデート
【題意】
2匹のカエルは所定の長さの数軸上で運動し、所定の初期位置とジャンプの毎回の長さを与え、いつ2匹のカエルが出会うことができるかを尋ねた.
【考え方】
題意によりx+mt=y+nt+kl,tはジャンプ回数を表す.式x−y=(n−m)*t+k*lを変化させてax+by=cのような形にしたので,拡張ユークリッドアルゴリズムを用いて最小ジャンプ回数を求めることを考慮した.
【コード】
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#define mem(a) memset(a,0,sizeof(a))
#define go(i,a,b) for(int i=a;i<=b;i++)
#define og(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e6;
const int inf=0x3f3f3f3f;
typedef long long ll;
typedef unsigned long long ull;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
    if(!b)
    {
        x=1;y=0;
        return a;
    }
    int r=exgcd(b,a%b,x,y);
    int t=x; x=y; y=t-a/b*y;
    return r;
}
int main()
{
    ll x,y,m,n,l,a,b,c,d;
    scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l);
    b=n-m; a=x-y; c=exgcd(b,l,x,y);
    if(a%c) printf("Impossible
"); else { x=x*a/c; x=(x%l+l)%l; printf("%d
",x); } return 0; }