//#include #include using namespace std; using lint = long long; constexpr lint mod = 1e9 + 7; #define all(x) begin(x), end(x) #define bitcount(n) __builtin_popcountll((lint)(n)) #define fcout cout << fixed << setprecision(15) #define highest(x) (63 - __builtin_clzll(x)) #define rep(i, n) for(int i = 0; i < int(n); i++) #define rep2(i, l, r) for(int i = int(l); i < int(r); i++) #define repr(i, n) for(int i = int(n) - 1; i >= 0; i--) #define repr2(i, l, r) for(int i = int(r) - 1; i >= int(l); i--) #define accumulate you cannot use this word! constexpr int inf9 = 1e9; constexpr lint inf18 = 1e18; inline void Yes(bool condition){ if(condition) cout << "Yes" << endl; else cout << "No" << endl; } template void array_output(itr start, itr goal){ for(auto i = start; i != goal; i++) cout << (i == start ? "" : " ") << (*i); cout << endl; } template void cins(itr first, itr last){ for(auto i = first; i != last; i++){ cin >> (*i); } } template T gcd(T a, T b){ if(b) return gcd(b, a % b); else return a; } template T lcm(T a, T b){ return a / gcd(a, b) * b; } template bool chmax(T &a, const T &b){ if(a < b){ a = b; return 1; } return 0; } template bool chmin(T &a, const T &b){ if(b < a){ a = b; return 1; } return 0; } inline int has(lint i, int j){ return (i >> j) & 1; } int dy[4] = {1, 0, -1, 0}; int dx[4] = {0, 1, 0, -1}; bool is_inside(lint y, lint x, lint H, lint W){ return (0 <= y && y < H && 0 <= x && x < W); } struct io_init { io_init() { cin.tie(nullptr); cout.tie(nullptr); std::ios::sync_with_stdio(false); } } io_init; tuple extgcd(lint a, lint b){ // Need: a, b >= 0 if(b == 0){ return {1, 0, a}; } auto [x, y, g] = extgcd(b, a % b); return {y, x - (a / b) * y, g}; } lint signed_mod(lint x, lint modulo) { return (x % modulo + modulo) % modulo; } lint modinv(lint a, lint modulo){ // Need: gcd(x, modulo) = 1 auto [x, y, g] = extgcd(a, modulo); return signed_mod(x, modulo); } pair CRT(lint a1, lint m1, lint a2, lint m2){ // return (a, m) a1 = signed_mod(a1, m1); a2 = signed_mod(a2, m2); auto [b, c, g] = extgcd(m1, m2); if(a1 % g != a2 % g) return {-1, -1}; lint m = m1 / g * m2; lint tmp = ((a2 - a1) / g) * b % (m2 / g); lint a = tmp * m1 + a1; a = signed_mod(a, m); return {a, m}; } pair Garner_simple(vector a, vector m, lint MOD = LONG_LONG_MAX){ // Need: m はそれぞれが互いに素 / return (a, m) int n = int(a.size()); rep(i, n) a[i] = signed_mod(a[i], m[i]); m.push_back(MOD); lint coefs[n + 1], val_sums[n + 1]; rep(i, n + 1) coefs[i] = 1; rep(i, n + 1) val_sums[i] = 0; rep(i, n){ lint t = a[i] - val_sums[i]; t = signed_mod(t, m[i]); t = t * modinv(coefs[i], m[i]) % m[i]; // cout << t << endl; rep2(j, i + 1, n + 1){ (val_sums[j] += t * coefs[j]) %= m[j]; (coefs[j] *= m[i]) %= m[j]; } } // rep(i, n + 1) cout << val_sums[i] << " " << coefs[i] << endl; return {val_sums[n], coefs[n]}; } pair Garner(vector a, vector m, lint MOD = LONG_LONG_MAX){ // return (a, m) int n = int(a.size()); rep(i, n) rep(j, i){ lint g = gcd(m[i], m[j]); if(a[i] % g != a[j] % g) return {-1, -1}; m[i] /= g; m[j] /= g; lint gi = gcd(m[i], g), gj = g / gi; while(true){ lint move = gcd(gi, gj); if(move == 1) break; gj /= move; gi *= move; } m[i] *= gi; m[j] *= gj; a[i] %= m[i], a[j] %= m[j]; } // rep(i, n) cout << a[i] << " " << m[i] << endl; // cout << endl; return Garner_simple(a, m, MOD); } int main(){ vector a(3), m(3); rep(i, 3){ lint x, y; cin >> x >> y; a[i] = x; m[i] = y; } auto [x, l] = Garner(a, m); if(x == -1){ cout << -1 << endl; }else{ cout << (x > 0 ? x : l) << endl; } }