結果

問題 No.28 末尾最適化
コンテスト
ユーザー chocorusk
提出日時 2019-01-26 23:08:59
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 307 ms / 5,000 ms
コード長 1,138 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,307 ms
コンパイル使用メモリ 103,340 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-19 10:38:57
合計ジャッジ時間 1,650 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
const ll MOD=100000009;
vector<P> f[37];
int main()
{
	for(int i=2; i<=36; i++){
		int i1=i;
		for(int j=2; j<=31; j++){
			if(i1%j==0){
				int e=0;
				while(i1%j==0){
					i1/=j;
					e++;
				}
				f[i].push_back(P(j, e));
			}
		}
	}
	int q; cin>>q;
	for(int t=0; t<q; t++){
		ll s; cin>>s;
		int n, k, b; cin>>n>>k>>b;
		int ans=1e9;
		for(int j=0; j<f[b].size(); j++){
			ll p=f[b][j].first; int e=f[b][j].second;
			vector<int> v;
			ll x=s;
			for(int i=0; i<=n; i++){
				int e1=0;
				ll x1=x;
				while(x1%p==0){
					x1/=p;
					e1++;
				}
				v.push_back(e1);
				x=1+(x*x+x*12345)%MOD;
			}
			sort(v.begin(), v.end());
			int ct=0;
			for(int i=0; i<k; i++) ct+=v[i];
			ans=min(ans, ct/e);
		}
		cout<<ans<<endl;
	}
	return 0;
}
0