#include using namespace std; using ll = long long; #define rep(i, s, t) for (ll i = s; i < (ll)(t); i++) #define all(x) begin(x), end(x) template bool chmin(T& x, T y) { return x > y ? (x = y, true) : false; } template bool chmax(T& x, T y) { return x < y ? (x = y, true) : false; } void solve() { ll n, a, b, c; cin >> n >> a >> b >> c; vector f(n), r(n); f[0] = 1; rep(i, 1, n) { r[i] = r[i - 1] | (f[i - 1] * i >= n); f[i] = f[i - 1] * i % n; } vector cost(n + 1, 1e18); cost[1] = 0; priority_queue, vector>, greater>> pq; pq.push({0, 1}); while(!pq.empty()) { auto [d, x] = pq.top(); pq.pop(); if(x == 0) break; if(x == n) { if(chmin(cost[0], d + c)) pq.push({d + c, 0}); continue; } if(cost[x] > d) continue; if(chmin(cost[(x + 1) % n], d + a)) pq.push({d + a, (x + 1) % n}); { ll cx = x; for(ll cb = b; cb < n * a; cb *= b) { if(cx >= n) if(chmin(cost[n], d + cb)) pq.push({d + cb, n}); cx %= n; if(chmin(cost[cx], d + cb)) pq.push({d + cb, cx}); cx = cx * x; } } if(chmin(cost[f[x]], d + c)) pq.push({d + c, f[x]}); if(r[x]) if(chmin(cost[n], d + c)) pq.push({d + c, n}); } cout << cost[0] << '\n'; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(15); int t = 1; // cin >> t; while (t--) solve(); }