#include<iostream>
#include<string>
#include<vector>
#include<utility>
#include<algorithm>
#include<map>
#include<cstdlib>
#include<cmath>
#include<numeric>
#include<iomanip>
#include<functional>
#include<cstdlib>
#include<queue>
const double PI = acos(-1);
#define rep(i,n)for(int i=0;i<(n);i++)
using namespace std;
using ll =long long;
const int mod = 1000000007;

//map<string>a;
class mint {
	long long x;
public:
	mint(long long x = 0) : x((x% mod + mod) % mod) {}
	mint operator-() const {
		return mint(-x);
	}
	mint& operator+=(const mint & a) {
		if ((x += a.x) >= mod) x -= mod;
		return *this;
	}
	mint & operator-=(const mint & a) {
		if ((x += mod - a.x) >= mod) x -= mod;
		return *this;
	}
	mint & operator*=(const  mint & a) {
		(x *= a.x) %= mod;
		return *this;
	}
	mint operator+(const mint & a) const {
		mint res(*this);
		return res += a;
	}
	mint operator-(const mint & a) const {
		mint res(*this);
		return res -= a;
	}
	mint operator*(const mint & a) const {
		mint res(*this);
		return res *= a;
	}
	mint pow(ll t) const {
		if (!t) return 1;
		mint a = pow(t >> 1);
		a *= a;
		if (t & 1) a *= *this;
		return a;
	}
	// for prime mod
	mint inv() const {
		return pow(mod - 2);
	}
	mint& operator/=(const mint & a) {
		return (*this) *= a.inv();
	}
	mint operator/(const mint & a) const {
		mint res(*this);
		return res /= a;
	}

	friend ostream& operator<<(ostream & os, const mint & m) {
		os << m.x;
		return os;
	}
};
int ctoi(char c) {
	if (c >= '0' && c <= '9') {
		return c - '0';
	}
	return 0;
}
//組み合わせ
void recursive_comb(vector<int> indexes, int s, int rest, std::function<void(vector<int>)> f) {
	if (rest == 0) {
		f(indexes);
	}
	else {
		if (s < 0) return;
		recursive_comb(indexes, s - 1, rest, f);
		indexes[rest - 1] = s;
		recursive_comb(indexes, s - 1, rest - 1, f);
	}
}

// nCkの組み合わせに対して処理を実行する
void foreach_comb(int n, int k, std::function<void(vector<int>)> f) {
	vector<int>indexes(k);
	recursive_comb(indexes, n - 1, k, f);
}

ll gcd(ll a, ll b) {
	if (a % b == 0) {
		return b;
	}
	else {
		return gcd(b, a % b);
	}
}
ll lcm(ll a, ll b) {
	return a * b / gcd(a, b);
}

int main() {

	/*ll N, M;
	cin >> N >> M;
	ll c, d;
	vector<pair<ll, ll>>A(M);
	rep(i, M) {*/
//		cin >> c >> d;
//		if (c > d) {
//			swap(c, d);
//		}
//		A[i] = make_pair(c, d);
//	}
//	sort(A.begin(), A.end());
//	vector<ll>ans(N,0);
//	ll tmp ;
//	vector<ll>B(N);
//	vector<ll>C(N);
//
//	rep(i, M) {
//		B[i] = A[i].first;
//		C[i] = A[i].second;
//	}
//	ll q,p=1;
//	rep(j, N) {
//		vector<ll>next;
//		if (j == 0) {
//			next.push_back(1);
//		}
//		for (int k = 0; k < next.size(); k++) {
//			tmp = next[k];
//
//			auto iter_lower = lower_bound(B.begin(), B.end(), tmp);
//			auto iter_upper = upper_bound(B.begin(), B.end(), tmp);
//			long idx_lower = distance(B.begin(), iter_lower);  // vec[idx_lower] = vec[7] = 6
//			long idx_upper = distance(B.begin(), iter_upper);  // vec[idx_upper] = vec[8] = 8
//
//			q = idx_lower - 1;
//			for (int i = 0; i < idx_upper - idx_lower; i++) {
//				if (ans[C[i + q]-1] == 0) {
//					ans[C[i + q]-1] = tmp;
//					next.push_back(C[i + q]);
//					p++;
//				}
//			}
//			 iter_lower = lower_bound(C.begin(), C.end(), tmp);
//			 iter_upper = upper_bound(C.begin(), C.end(), tmp);
//			 idx_lower = distance(C.begin(), iter_lower);  // vec[idx_lower] = vec[7] = 6
//			 idx_upper = distance(C.begin(), iter_upper);  // vec[idx_upper] = vec[8] = 8
//			q = idx_lower - 1;
//
//			for (int i = 0; i < idx_upper - idx_lower; i++) {
//				if (ans[B[i + q]-1] == 0) {
//					ans[B[i + q]-1] = tmp;
//					next.push_back(B[i + q]);
//					p++;
//				}
//			}
//		}
//		if (p == N) {
//			break;
//		}
//	}
//	cout << "Yes"<<endl;
//	rep(i, N-1) {
//		cout << ans[i+1]<<endl;
//	}
//}




	ll A,B,C,D;
	cin >> A >> B >> C >> D;
	cout << (B - A + 1) - B / C - B / D + B / lcm(C, D) + (A - 1) / C + (A - 1) / D - (A - 1) / lcm(C, D);

}

Measured memory_bytes: 256