結果

問題 No.28 末尾最適化
ユーザー furafura
提出日時 2020-01-09 12:57:41
言語 C++14
(gcc 10.2.0 + boost 1.73.0)
結果
AC  
実行時間 210 ms / 5,000 ms
コード長 1,052 Byte
コンパイル時間 1,413 ms
使用メモリ 3,712 KB
最終ジャッジ日時 2021-06-06 18:52:13
合計ジャッジ時間 2,464 ms
ジャッジサーバーID
(参考情報)
judge6 / judge10
このコードへのチャレンジ(β)

テストケース

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

ソースコード

diff #

#include <bits/stdc++.h>

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

using namespace std;
using lint=long long;

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);
	vector<lint> a(n+1);
	a[0]=seed;
	rep(i,n) a[i+1]=1+(a[i]*a[i]+a[i]*12345)%100000009;

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

		int hist[70]={};
		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,70){
			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