結果

問題 No.2744 Power! or +1
ユーザー hiikunZ
提出日時 2024-04-20 22:17:16
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,996 bytes
コンパイル時間 8,317 ms
コンパイル使用メモリ 268,004 KB
最終ジャッジ日時 2025-02-21 07:27:58
ジャッジサーバーID
（参考情報） 		judge3 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 2
other TLE * 1 -- * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
//#include<atcoder/all>
//using namespace atcoder;
using ll = long long int;
using ull = unsigned long long int;
using ld = long double;
constexpr ll MAX = 2000000000000000000;
constexpr ld PI = 3.14159265358979;
constexpr ll MOD = 0;//2024948111;
ld dotorad(ld K){ return PI * K / 180.0; }
ld radtodo(ld K){ return K * 180.0 / PI; }
mt19937 mt;
void randinit(){ srand((unsigned)time(NULL));mt = mt19937(rand()); }


int main(){
    ll N,A,B,C;
    cin >> N >> A >> B >> C;
    vector<ll> fac(N + 2,1),fac_b(N + 2,1);
    for(ll i = 1;i <= N;i++) fac[i] = (fac[i - 1] * i) % N;
    for(ll i = 1;i <= N;i++){
        fac_b[i] = (fac_b[i - 1] * i);
        if(fac_b[i] > N) fac_b[i] = N;
    }

    priority_queue<tuple<ll,ll>,vector<tuple<ll,ll>>,greater<tuple<ll,ll>>> que;
    que.emplace(0,1);
    vector<ll> memo(N + 2,MAX);
    while(!que.empty()){
        auto [cost,now] = que.top();
        que.pop();
        if(memo[now] == MAX){
            memo[now] = cost;
            que.emplace(cost + A,(now + 1) % N);
            ll x = 1,c = 1;
            for(ll i = 1;i < 35;i++){
                x *= now;
                x %= N;
                if(A * N / c < B) break;
                c *= B;
                que.emplace(cost + c,x);
            }
            que.emplace(cost + C,fac[now]);
        }
    }
    vector<ll> memo2(N + 2,MAX);
    que.emplace(0,1);
    while(!que.empty()){
        auto [cost,now] = que.top();
        que.pop();
        if(memo2[now] == MAX){
            memo2[now] = cost;
            que.emplace(cost + A,(now + 1));
            ll x = 1,c = 1;
            for(ll i = 1;i < 60;i++){
                x *= now;
                x = min(x,N);
                if(A * N / c < B) break;
                c *= B;
                que.emplace(cost + c,min(x,N));
                if(x >= N) break;
            }
            que.emplace(cost + C,fac_b[now]);
        }
    }
    cout << min(memo[0],memo2[N] + C) << endl;
}
0