結果

問題 No.2744 Power! or +1
ユーザー cho435
提出日時 2024-04-20 04:05:57
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,257 bytes
コンパイル時間 2,492 ms
コンパイル使用メモリ 202,872 KB
最終ジャッジ日時 2025-02-21 06:01:52
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 8 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define rep(i,n) for(int i=0;i<(int)(n);i++)

bool ck(ll x,ll k,ll m){
	double dx=x;
	rep(i,k) dx*=x;
	return dx>m;
}

int main(){
	ll n,a,b,c;
	cin>>n>>a>>b>>c;
	vector<ll> fac(n,1);
	ll ud=-1;
	for(ll i=2;i<n;i++){
		if(ud==-1){
			fac.at(i)=fac.at(i-1)*i;
			if(fac.at(i)>=n){
				ud=i;
				fac.at(i)%=n;
			}
		}else{
			fac.at(i)=fac.at(i-1)*i%n;
		}		
	}
	ll mxc=n*a;
	vector<ll> cs(n+1,1e18);
	priority_queue<pair<ll,ll>,vector<pair<ll,ll>>,greater<pair<ll,ll>>> pq;
	pq.push({0,1});
	while(!pq.empty()){
		auto[sc,v]=pq.top();
		pq.pop();
		if(cs.at(v)!=1e18) continue;
		cs.at(v)=sc;
		if(v==n) break;
		if(v==0){
			pq.push({sc+c,n});
			continue;
		}
		if(cs.at(v+1)==1e18){
			pq.push({sc+a,v+1});
		}
		{
			ll nx=1;
			ll pw=1;
			int flg=0;
			for(int k=1;1;k++){
				nx=(nx*v);
				if(nx>n){
					nx%=n;
					flg=1;
				}
				if(ck(b,k,mxc)) break;
				pw*=b;
				if(cs.at(nx)==1e18){
					pq.push({sc+pw,nx});
				}
				if(flg&&cs.at(0)==1e18){
					pq.push({sc+pw,0});
				}
			}
		}
		{
			int nx=fac.at(v);
			if(nx==0) nx=n;
			if(cs.at(nx)==1e18){
				pq.push({sc+c,nx});
			}
			if(ud<=v&&cs.at(0)==1e18){
				pq.push({sc+c,0});
			}
		}
	}
	cout<<cs.at(n)<<endl;
}
0