#include using namespace std; struct Initializer { Initializer() { cin.tie(0); ios::sync_with_stdio(0); cout << fixed << setprecision(15); } } initializer; template inline istream& operator>>(istream &s, vector &v) { for (T &t : v) s >> t; return s; } template inline ostream& operator<<(ostream &s, const vector &v) { for (const T &t : v) s << t << endl; return s; } template inline T min(vector& v) {return *min_element(v.begin(), v.end());} template inline T max(vector& v) {return *max_element(v.begin(), v.end());} template inline int min_element(vector& v) {return min_element(v.begin(), v.end()) - v.begin();} template inline int max_element(vector& v) {return max_element(v.begin(), v.end()) - v.begin();} template inline void sort(vector& v) {sort(v.begin(), v.end());} template inline void sort(vector& v, Function func) {sort(v.begin(), v.end(), func);} template inline void rsort(vector& v) {sort(v.rbegin(), v.rend());} template inline void reverse(vector& v) {reverse(v.begin(), v.end());} template inline void unique(vector& v) {v.erase(unique(v.begin(), v.end()), v.end());} template inline void nth_element(vector& v, int n) {nth_element(v.begin(), v.begin() + n, v.end());} template inline bool next_permutation(vector& v) {return next_permutation(v.begin(), v.end());} template inline int find(vector& v, T t) {return find(v.begin(), v.end(), t) - v.begin();} template inline int in(vector v, T t) {return find(v, t) != (int)v.size();} template inline int lower_bound(vector& v, T t) {return lower_bound(v.begin(), v.end(), t) - v.begin();} template inline int upper_bound(vector& v, T t) {return upper_bound(v.begin(), v.end(), t) - v.begin();} template inline T accumulate(const vector& v, function func = plus()) {return accumulate(v.begin(), v.end(), T(), func);} template inline void adjacent_difference(vector& v) {adjacent_difference(v.begin(), v.end(), v.begin());} template inline void adjacent_difference(vector& v, vector& u) {adjacent_difference(v.begin(), v.end(), u.begin());} template inline void partial_sum(vector& v, vector& u) {partial_sum(v.begin(), v.end(), u.begin());} template inline T inner_product(vector& v, vector& u) {return inner_product(v.begin(), v.end(), u.begin(), T(0));} template inline int count(const vector& v, T t) {return count(v.begin(), v.end(), t);} template inline int count_if(const vector& v, Function func) {return count_if(v.begin(), v.end(), func);} template inline void remove_if(vector& v, Function func) {v.erase(remove_if(v.begin(), v.end(), func), v.end());} template inline bool any_of(vector v, Function func) {return any_of(v.begin(), v.end(), func);} template inline vector subvector(vector& v, int a, int b) {return vector(v.begin() + a, v.begin() + b);} template inline int kinds(const vector& v) {return set(v.begin(), v.end()).size();} template inline T gcd(T t) { return t; } template inline T gcd(T t, S... s) { return __gcd(t, gcd(s...)); } template inline T lcm(T t) { return t; } template inline T lcm(T t, S... s) { T l = lcm(s...); return t / gcd(t, l) * l; } template inline T floor(T a, T b) { return a / b * b <= a ? a / b : a / b - 1; } template<> inline float floor(float a, float b) { return floor(a / b) * b <= a ? floor(a / b) : floor(a / b) - 1; } template<> inline double floor(double a, double b) { return floor(a / b) * b <= a ? floor(a / b) : floor(a / b) - 1; } template<> inline long double floor(long double a, long double b) { return floor(a / b) * b <= a ? floor(a / b) : floor(a / b) - 1; } template inline T ceil(T a, T b) { return floor(a + b - 1, b); } template inline T round(T a, T b) { return floor(a + b / 2); } template inline T mod(T a, T b) { return a - floor(a, b) * b; } template inline T factorial(T n) { return n <= 1 ? 1 : factorial(n - 1) * n; } template inline T square(T n) { return n * n; } template inline T cube(T n) { return n * n * n; } template inline T norm(T x1, T y1, T x2, T y2) { return square(x1 - x2) + square(y1 - y2); } inline long long sqrt(long long n) { return sqrt((long double)n); } int main() { int n, res = 0; vector a(3); cin >> n >> a; for (int i : a) res += n / i; for (int i = 0; i < 3; ++i) res -= n / lcm(a[i], a[(i + 1) % 3]); cout << res + n / lcm(a[0], a[1], a[2]) << endl; }