hihocoder 1257 扩展欧几里得

1297 : 数论四·扩展欧几里德

时间限制:10000ms
单点时限:1000ms
内存限制:256MB
描述

小Hi和小Ho周末在公园溜达。公园有一堆围成环形的石板，小Hi和小Ho分别站在不同的石板上。已知石板总共有m块，编号为 0..m-1，小Hi一开始站在s1号石板上，小Ho一开始站在s2号石板上。

小Hi：小Ho，你说我们俩如果从现在开始按照固定的间隔数同时同向移动，我们会不会在某个时间点站在同一块石板上呢？

小Ho：我觉得可能吧，你每次移动v1块，我移动v2块，我们看能不能遇上好了。

小Hi：好啊，那我们试试呗。

一个小时过去了，然而小Hi和小Ho还是没有一次站在同一块石板上。

小Ho：不行了，这样走下去不知道什么时候才汇合。小Hi，你有什么办法算算具体要多久才能汇合么？

小Hi：让我想想啊。。

提示：扩展欧几里德

输入

第1行：每行5个整数s1,s2,v1,v2,m，0≤v1,v2≤m≤1,000,000,000。0≤s1,s2 < m

中间过程可能很大，最好使用64位整型

输出

第1行：每行1个整数，表示解，若该组数据无解则输出-1

样例输入
0 1 1 2 6
样例输出
5

解析跟POJ 1061 一样：http://blog.csdn.net/qq_32866009/article/details/51558587

下面是AC代码：

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;long long exgcd(long long m,long long &x,long long n,long long &y){    long long x1,y1,x0,y0;    x0=1,y0=0;    x1=0,y1=1;    long long r=(m%n+n)%n;    long long q=(m-r)/n;    x=0,y=1;    while(r)    {        x=x0-q*x1;        y-y0-q*y1;        x0=x1;        y0=y1;        x1=x;        y1=y;        m=n;        n=r;        r=m%n;        q=(m-r)/n;    }    return n;}int main(){    long long r,t,yy,x,y,m,n,l,ar,br,cr;    while(~scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l))    {            long long M=exgcd(n-m,ar,l,br);            if((x-y)%M||m==n)            {                printf("-1\n");            }            else            {                long long s=l/M;                ar=ar*((x-y)/M);                ar=(ar%s+s)%s;                printf("%lld\n",ar);            }    }    return 0;}