//#define _GLIBCXX_DEBUG #include #define rep(i, n) for(int i=0; i; using vs = vector; using vi = vector; using vvi = vector; template using PQ = priority_queue; template using PQG = priority_queue, greater >; const int INF = 100010001; const ll LINF = (ll)INF*INF*10; template inline bool chmax(T1 &a, T2 b) {return a < b && (a = b, true);} template inline bool chmin(T1 &a, T2 b) {return a > b && (a = b, true);} template istream &operator>>(istream &is, pair &p) { return is >> p.first >> p.second;} template ostream &operator<<(ostream &os, const pair &p) { return os << p.first << ' ' << p.second;} int64_t gcd(int64_t a, int64_t b) {return b?gcd(b, a%b):a;} int64_t lcm(int64_t a, int64_t b) {return a/gcd(a,b)*b;} int64_t extgcd(int64_t a, int64_t b, int64_t &x, int64_t &y) { // ax+by=gcd(a, b)の解を与える int64_t d = a; if(b) { d = extgcd(b, a%b, y, x); y -= (a/b)*x; } else { x = 1; y = 0; } return d; } int main() { ios::sync_with_stdio(false); cin.tie(0); ll n, a[3]; cin >> n; rep(i, 3) cin >> a[i]; int dp[8] = {}; rep(i, 8) { if(!i) continue; ll x = 1; rep(j, 3) if(i>>j&1) { x = lcm(x, a[j]); } dp[i] = n/x; } rep(i, 3) rep(j, 8) { if(j>>i&1) dp[j] -= dp[j-(1<