//#define _GLIBCXX_DEBUG #include #define rep(i, n) for(int i=0; i; using vs = vector; using vi = vector; using vvi = vector; template using PQ = priority_queue; template using PQG = priority_queue, greater >; const int INF = 0xccccccc; const ll LINF = 0xcccccccccccccccLL; template inline bool chmax(T1 &a, T2 b) {return a < b && (a = b, true);} template inline bool chmin(T1 &a, T2 b) {return a > b && (a = b, true);} template istream &operator>>(istream &is, pair &p) { return is >> p.first >> p.second;} template ostream &operator<<(ostream &os, const pair &p) { return os << p.first << ' ' << p.second;} struct barrett { unsigned int _m; uint64_t im; barrett(unsigned int m):_m(m), im(uint64_t(-1)/m + 1) {} unsigned int umod() const {return _m;} unsigned int mul(unsigned int a, unsigned int b) const { uint64_t z = a; z *= b; uint64_t x = uint64_t(((unsigned __int128)(z)*im) >> 64); unsigned int v = (unsigned int)(z-x*_m); if(_m <= v) v += _m; return v; } }; constexpr int64_t safe_mod(int64_t x, int64_t m) { x %= m; if(x < 0) x += m; return x; } constexpr pair inv_gcd(int64_t a, int64_t b) { a = safe_mod(a, b); if(a == 0) return {b, 0}; int64_t s = b, t = a; int64_t m0 = 0, m1 = 1; while(t) { int64_t u = s/t; s -= t*u; m0 -= m1*u; int64_t tmp = s; s = t; t = tmp; tmp = m0; m0 = m1; m1 = tmp; } if(m0 < 0) m0 += b/s; return {s, m0}; } int64_t pow_mod(int64_t x, int64_t n, int m) { if(m == 1) return 0; barrett bt((unsigned int)(m)); unsigned int r = 1, y = (unsigned int)(safe_mod(x, m)); while(n) { if(n&1) r = bt.mul(r, y); y = bt.mul(y, y); n >>= 1; } return r; } int64_t inv_mod(int64_t x, int64_t m) { pair z = inv_gcd(x, m); assert(z.first == 1); return z.second; } //(y, z) answer mod z = y, 無理な時は(-1, -1) pair crt(const vector &r, const vector &m) { int n = int(r.size()); int64_t r0 = 0, m0 = 1; for(int i = 0; i < n; i++) { assert(1 <= m[i]); int64_t r1 = safe_mod(r[i], m[i]), m1 = m[i]; if(m0 < m1) { swap(r0, r1); swap(m0, m1); } if(m0 % m1 == 0) { if(r0 % m1 != r1) return {-1, -1}; continue; } int64_t g, im; tie(g, im) = inv_gcd(m0, m1); int64_t u1 = (m1/g); if((r1-r0)%g) return {-1, -1}; int64_t x = (r1-r0)/g%u1 * im % u1; r0 += x*m0; m0 *= u1; if(r0 < 0) r0 += m0; } return {r0, m0}; } //[0, n) (a*i+b)/m int64_t floor_sum(int64_t n, int64_t m, int64_t a, int64_t b) { int64_t ans = 0; if(a >= m) { ans += ((n-1)*n>>1)*(a/m); a %= m; } if(b >= m) { ans += n*(b/m); b %= m; } int64_t y_max = (a*n+b)/m, x_max = y_max*m-b; if(y_max == 0) return ans; ans += (n - (x_max+a-1)/a) * y_max; ans += floor_sum(y_max, a, m, (a - x_max%a)%a); return ans; } //head ll ans(pair x) { return x.first?x.first:x.second; } int main() { ios::sync_with_stdio(false); cin.tie(0); vector r(3), m(3); rep(i, 3) cin >> r[i] >> m[i]; cout << ans(crt(r, m)) << endl; }