#include #define int long long #define uint unsigned int #define REP(i, n) for (int i = 0; i < (n); ++i) #define REPR(i, n) for (int i = (n); i >= 0; --i) #define FOR(i, a, b) for (int i = (a); i < (b); ++i) #define FORR(i, a, b) for (int i = (a); i >= (b); --i) #define ALL(a) (a).begin(), (a).end() #define RALL(a) (a).rbegin(), (a).rend() #define SZ(a) ((int)(a).size()) #define PB(a) push_back(a) #define EMP emplace #define EMPB(...) emplace_back(__VA_ARGS__) #define MP(a, b) make_pair(a, b) #define MT(...) make_tuple(__VA_ARGS__) #define Bit(n) (1LL << (n)) using namespace std; using pii = pair; template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } const int MOD = 1000000007; const int INF = 1LL << 30; const double EPS = 1e-10; int gcd(int a,int b) { return b?gcd(b,a%b):a; } int ext_gcd(int a,int b,int& x,int& y) { if(b==0){ x=1;y=0;return a; } int q=a/b; int g=ext_gcd(b,a-q*b,x,y); int z=x-q*y; x=y;y=z; return g; } int inv_mod(int a,int m) { int x,y; ext_gcd(a,m,x,y); x%=m; return (m+x%m)%m; } pair linear_congruence(const vector& A, const vector& B, const vector& M) { int x = 0, m = 1; for (int i = 0; i < A.size(); i++) { int a = A[i] * m, b = B[i] - A[i] * x, d = gcd(M[i], a); if (b % d != 0) return make_pair(0, -1); int t = b / d * inv_mod(a / d, M[i] / d) % (M[i] / d); x = x + m * t; m *= M[i] / d; } return make_pair((m + x % m) % m, m); } signed main() { cin.tie(0); ios_base::sync_with_stdio(false); cout << fixed << setprecision(10); vector A, B, M; REP(i, 3) { int X, Y; cin >> X >> Y; for (int p = 2; p*p <= Y; p++) { int tmp = 1; while (Y % p == 0) { Y /= p; tmp *= p; } if (tmp != 1) { A.PB(1); B.PB(X%tmp); M.PB(tmp); } } if (Y != 1) { A.PB(1); B.PB(X%Y); M.PB(Y); } } auto p = linear_congruence(A, B, M); if (p.second < 0) cout << -1 << endl; else cout << p.first << endl; return 0; }