結果

問題 No.28 末尾最適化
コンテスト
ユーザー ytft
提出日時 2022-11-03 21:55:17
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 350 ms / 5,000 ms
コード長 827 bytes
コンパイル時間 2,909 ms
コンパイル使用メモリ 249,192 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-18 04:32:40
合計ジャッジ時間 3,248 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
long long mod=100000009;
long long primes[]={2,3,5,7,11,13,17,19,23,29,31};
void solve(){
	long long seed,N,B,K;
	scanf("%lld%lld%lld%lld",&seed,&N,&K,&B);
	long long X[N+1],p[N+1],ans=INT_MAX,count,elem[11];
	X[0]=seed;
	for(int i=0;i<11;++i){
		elem[i]=0;
		while(B%primes[i]==0){
			elem[i]+=1;
			B/=primes[i];
		}
	}
	for(int i=0;i<N;++i){
		X[i+1]=1+(X[i]*(X[i]+12345ll))%mod;
	}
	for(int i=0;i<11;++i){
		if(elem[i]==0){
			continue;
		}
		for(int j=0;j<N+1;++j){
			p[j]=0;
			long long temp=X[j];
			while(temp%primes[i]==0){
				temp/=primes[i];
				p[j]+=1;
			}
		}
		sort(p,p+N+1);
		count=0;
		for(int j=0;j<K;++j){
			count+=p[j];
		}
		ans=min(ans,count/elem[i]);
	}
	printf("%d\n",ans);
}
int main(){
	int T;
	scanf("%d",&T);
	for(int i=0;i<T;++i){
		solve();
	}
}
0