結果

問題 No.2744 Power! or +1
ユーザー hiikunZhiikunZ
提出日時 2024-04-20 22:12:41
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 1,958 bytes
コンパイル時間 2,241 ms
コンパイル使用メモリ 214,112 KB
実行使用メモリ 283,796 KB
最終ジャッジ日時 2024-10-12 20:30:17
合計ジャッジ時間 10,657 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 TLE -
testcase_01 -- -
testcase_02 -- -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
権限があれば一括ダウンロードができます

ソースコード

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 < 60;i++){
                x *= now;
                x %= N;
                if(MAX / 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(MAX / c < B) break;
                c *= B;
                que.emplace(cost + c,min(x,N));
            }
            que.emplace(cost + C,fac_b[now]);
        }
    }
    cout << min(memo[0],memo2[N] + C) << endl;
}
0