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

hihocoder 1297

描述

小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

Solution

通过扩展欧几里得求得方程(v1-v2)*t+k*m=(s2-s1)的解t，其中k为超过的圈数

Code

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <stack>
#include <map>
#include <vector>
#include <queue>
#define LL long long
#define L 1000010
using namespace std;

LL s1, s2, v1, v2, m, x, y;

inline int gcd(LL a, LL b) {
  return (a % b == 0) ? b : gcd(b, a % b);
}

inline void exgcd(LL a, LL b, LL& x, LL& y) {
  if (!b) {x = 1, y = 0;}
  else {
    exgcd(b, a % b, y, x);
    y -= (a / b) * x;
  }
}

int main() {
  scanf("%lld %lld %lld %lld %lld", &s1, &s2, &v1, &v2, &m);
  LL a = v1 - v2, b = m, c = s2 - s1;
  if (a < 0) a += m;
  LL d = gcd(a, b);
  if (c % d) printf("-1\n");
  else {
    a /= d, b /= d, c /= d;
    exgcd(a, b, x, y);
    x = x * c % b;
    if (x < 0) x += b;
    printf("%lld\n", x);
  }
  return 0;
}

Summary

通过扩展欧几里得求得的是经过变形后的方程的一组解，需要先扩展解系后在回归到原先的方程，求符合要求的解