#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; typedef pair p_ii; const int INF = 1e9; const double PI = acos(-1.0); const ll MOD = 1e9 + 7; ll gcd(ll a, ll b){ if(b==0)return a; else return gcd(b, a%b); } ll lcm(ll a, ll b){ ll d = gcd(a, b); return a*b/d; } int main() { ll N; cin>>N; vector a(3); cin>>a[0]>>a[1]>>a[2]; ll res=0; res+=N/a[0]+N/a[1]+N/a[2]; ll d1 = lcm(a[0], a[1]), d2=lcm(a[1], a[2]), d3=lcm(a[2], a[0]), d=lcm(lcm(a[0], a[1]), a[2]); res-=N/d1+N/d2+N/d3; res+=N/d; cout<