#include using namespace std; __int128_t safemod(__int128_t a,__int128_t m){ a %= m; a += m; return a>=m?a-m:a; } __int128_t floor(__int128_t a,__int128_t b){ if(b < 0) a *= -1,b *= -1; return a<0?(a+1)/b-1:a/b; } __int128_t ceil(__int128_t a,__int128_t b){ if(b < 0) a *= -1,b *= -1; return a>0?(a-1)/b+1:a/b; } __int128_t floorsum(__int128_t N,__int128_t M,__int128_t a,__int128_t b){ //return Σ(i=0~N-1)floor((a*i+b)/M). __int128_t ret = 0; if(a >= M) ret += a/M*(N*(N-1)/2),a %= M; if(b >= M) ret += b/M*N,b %= M; __int128_t maxy = a*N+b; if(maxy >= M) ret += floorsum(maxy/M,a,M,maxy%M); return ret; } __int128_t floorsum2(__int128_t N,__int128_t M,__int128_t a,__int128_t b){ if(N <= 0) return 0; //a < 0 or b < 0 or N <= 0; __int128_t ret = 0; if(a < 0 || a >= M) ret += floor(a,M)*(N*(N-1)/2),a = safemod(a,M); if(b < 0 || b >= M) ret += floor(b,M)*N,b = safemod(b,M); __int128_t maxy = a*N+b; if(maxy >= M) ret += floorsum(maxy/M,a,M,maxy%M); return ret; } int main(){ ios_base::sync_with_stdio(false); cin.tie(nullptr); long long A,B,C,D; cin >> A >> B >> C >> D; __int128_t g = gcd(A,B); A /= g,B /= g; g = gcd(C,D); C /= g,D /= g; if(A == C && B == D){cout << "-1\n"; return 0;} if(A*D > C*B) swap(A,C),swap(B,D); __int128_t lcm = B*D,a = A*D,c = C*B; __int128_t n = (lcm+c-a-1)/(c-a); __int128_t answer = n; answer += floorsum(n,2*B,2*A,2*A+B); answer -= floorsum(n,2*D,2*C,2*C+D); cout << (long long)answer << endl; }