#include using namespace std; using ll = long long; int main() { ll N, A, B, C; cin >> N >> A >> B >> C; vector fac(N, array()); fac[0] = {1, 0}; for (int i = 1; i < N; i++) { fac[i][1] = fac[i - 1][1]; fac[i][0] = fac[i - 1][0] * i % N; if (fac[i - 1][0] * i >= N) fac[i][1] = 1; } vector dist(N, array({(ll)1e18, (ll)1e18})); dist[1][0] = 0; priority_queue>, vector>>, greater>>> pq; pq.push({0, {1, 0}}); ll mCOST = (N - 1) * A; while (pq.size()) { auto [d, v] = pq.top(); pq.pop(); // 操作1 if (v[0] + 1 == N) { if (dist[(v[0] + 1) % N][1] > d + A) { dist[(v[0] + 1) % N][1] = d + A; pq.push({d + A, {(v[0] + 1) % N, 1}}); } } else { if (dist[(v[0] + 1) % N][v[1]] > d + A) { dist[(v[0] + 1) % N][v[1]] = d + A; pq.push({d + A, {(v[0] + 1) % N, v[1]}}); } } // 操作2 { ll b = B, v2 = v[0]; int ov_N = v[1]; while (b < mCOST) { if (dist[v2][ov_N] > d + b) { dist[v2][ov_N] = d + b; pq.push({d + b, {v2, ov_N}}); } b *= B; v2 *= v[0]; if (v2 >= N) { v2 %= N; ov_N = 1; } } } // 操作3 if (v[1]) { if (dist[0][1] > d + C) { dist[0][1] = d + C; pq.push({d + C, {0, 1}}); } } else { if (dist[fac[v[0]][0]][fac[v[0]][1]] > d + C) { dist[fac[v[0]][0]][fac[v[0]][1]] = d + C; pq.push({d + C, {fac[v[0]][0], fac[v[0]][1]}}); } } } cout << min(dist[0][0], dist[0][1]) << '\n'; }