結果

問題 No.28 末尾最適化
ユーザー msm1993msm1993
提出日時 2020-03-22 12:45:30
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 115 ms / 5,000 ms
コード長 1,030 bytes
コンパイル時間 1,526 ms
コンパイル使用メモリ 172,496 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-26 03:36:01
合計ジャッジ時間 2,123 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 115 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#define rep(i,n) for(int i=0;i<(n);i++)

using namespace std;

const int INF=1<<29;

map<long long,int> prime_factorize(long long a){
	map<long long,int> res;
	if(a%2==0){
		int cnt=0;
		do{ cnt++; a/=2; }while(a%2==0);
		res.emplace(2,cnt);
	}
	for(long long p=3;p*p<=a;p+=2) if(a%p==0) {
		int cnt=0;
		do{ cnt++; a/=p; }while(a%p==0);
		res.emplace(p,cnt);
	}
	if(a>1) res.emplace(a,1);
	return res;
}

void solve(){
	int seed,n,k,base; scanf("%d%d%d%d",&seed,&n,&k,&base);

	int a[10001];
	a[0]=seed;
	rep(i,n) a[i+1]=1+(1LL*a[i]*a[i]+a[i]*12345LL)%100000009;

	int ans=INF;
	for(auto q:prime_factorize(base)){
		int p=q.first;

		int hist[30]={};
		rep(i,n+1){
			int cnt=0;
			while(a[i]%p==0) a[i]/=p, cnt++;
			hist[cnt]++;
		}

		int sum=0,rest=k;
		rep(x,30){
			if(rest>hist[x]){
				sum+=hist[x]*x;
				rest-=hist[x];
			}
			else{
				sum+=rest*x;
				break;
			}
		}
		ans=min(ans,sum/q.second);
	}
	printf("%d\n",ans);
}

int main(){
	int q; scanf("%d",&q); rep(_,q) solve();
	return 0;
}
0