結果

問題 No.2744 Power! or +1
ユーザー Kude
提出日時 2024-04-21 12:38:53
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 277 ms / 3,000 ms
コード長 1,975 bytes
コンパイル時間 3,355 ms
コンパイル使用メモリ 277,740 KB
実行使用メモリ 6,908 KB
最終ジャッジ日時 2024-10-13 11:09:27
合計ジャッジ時間 4,733 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 9
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include<bits/stdc++.h>
namespace {
#pragma GCC diagnostic ignored "-Wunused-function"
#include<atcoder/all>
#pragma GCC diagnostic warning "-Wunused-function"
using namespace std;
using namespace atcoder;
#define rep(i,n) for(int i = 0; i < (int)(n); i++)
#define rrep(i,n) for(int i = (int)(n) - 1; i >= 0; i--)
#define all(x) begin(x), end(x)
#define rall(x) rbegin(x), rend(x)
template<class T> bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } else return false; }
template<class T> bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } else return false; }
using ll = long long;
using P = pair<int,int>;
using VI = vector<int>;
using VVI = vector<VI>;
using VL = vector<ll>;
using VVL = vector<VL>;
using mint = modint;
} int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  int n, a, b, c;
  cin >> n >> a >> b >> c;
  int fact_int[10];
  fact_int[0] = 1;
  for (int i = 1; i < 10; i++) fact_int[i] = fact_int[i-1] * i;
  mint::set_mod(n);
  vector<mint> fact(n);
  fact[0] = 1;
  for (int i = 1; i < n; i++) fact[i] = fact[i-1] * i;
  int ub = 2 * a + 4 * c;
  VI dp(n, ub);
  priority_queue<P, vector<P>, greater<P>> q;
  auto add = [&](int u, int cost) {
    if (chmin(dp[u], cost)) q.emplace(cost, u);
  };
  add(1, 0);
  while (q.size()) {
    auto [cost, u] = q.top(); q.pop();
    add((u + 1) % n, cost + a);
    {
      int bk = b;
      mint uk = u;
      while (!__builtin_mul_overflow(bk, b, &bk) && bk < ub) {
        uk *= u;
        add(uk.val(), cost + bk);
      }
    }
    add(fact[u].val(), cost + c);
  }
  int ans = dp[0];
  rep(i, n) if (i >= 10 || fact_int[i] >= n) chmin(ans, dp[i] + 2 * c);
  for (int x = 2; x < n; x++) {
    int bk = b;
    int xk = x;
    while (true) {
      if (__builtin_mul_overflow(xk, x, &xk)) xk = n;
      if (__builtin_mul_overflow(bk, b, &bk) || bk >= ub) break;
      if (xk >= n) {
        chmin(ans, dp[x] + bk + c);
        break;
      }
    }
  }
  cout << ans << '\n';
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0