結果

問題 No.28 末尾最適化
ユーザー kyuridenamidakyuridenamida
提出日時 2016-02-11 08:41:46
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 239 ms / 5,000 ms
コード長 1,005 bytes
コンパイル時間 2,474 ms
コンパイル使用メモリ 169,288 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-21 22:52:21
合計ジャッジ時間 2,336 ms
ジャッジサーバーID
(参考情報)
judge14 / judge13
このコードへのチャレンジ(β)

テストケース

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

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for(int (i) = 0 ; (i) < (int)(n) ; (i)++)
#define REP(i,a,b) for(int (i) = a ; (int)(i) <= (int)(b) ; (i)++)
#define all(n) (n).begin(),(n).end()
typedef vector<int> Vi;
typedef vector<Vi> VVi;
typedef pair<int,int> Pii;
typedef vector<Pii> VPii;


void solve(vector<int> A,int K,int B){
	vector< pair<int,int> > p;
	REP(i,2,B){
		if( B % i == 0 ){
			int c = 0;
			while( B % i == 0 ) B /= i, c++;
			p.push_back({i,c});
		}
	}
	int r = 1e9;
	for( auto target : p ){
		vector<int> C;
		rep(i,A.size()){
			int t = A[i];
			int c = 0;
			while( t % target.first == 0 ) t /= target.first, c++;
			C.push_back(c);
		}
		sort(all(C));
		int tot = 0;
		rep(i,K) tot += C[i];
		r = min(r,tot / target.second );
	}
	cout << r << endl;
	
}
int main(){
	int Q;
	cin >> Q;
	rep(i,Q){
		int s,N,K,B;
		cin >> s >> N >> K >> B;
		Vi A(N+1);
		A[0] = s;
		REP(j,1,N) A[j] = 1 + (1ll * A[j-1] * A[j-1] + A[j-1] * 12345ll) % 100000009;
		solve(A,K,B);
	}
}
0