結果

問題 No.28 末尾最適化
ユーザー assy1028assy1028
提出日時 2015-03-22 19:11:04
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 374 ms / 5,000 ms
コード長 1,854 bytes
コンパイル時間 853 ms
コンパイル使用メモリ 98,008 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-11 09:25:37
合計ジャッジ時間 1,813 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ（β）

テストケース

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

ソースコード

diff # 

#include <algorithm>
#include <iostream>
#include <cstdio>
#include <map>
#include <numeric>
#include <cmath>
#include <set>
#include <sstream>
#include <string>
#include <vector>
#include <queue>
#include <complex>
#include <string.h>
#include <unordered_set>
#include <unordered_map>
using namespace std;

#define endl '\n'
#define all(v) (v).begin(), (v).end()
#define rall(v) (v).rbegin(), (v).rend()
#define uniq(v) (v).erase(unique((v).begin(), (v).end()), (v).end())

typedef long long ll;
typedef pair<int, int> P;
typedef unsigned int uint;
typedef unsigned long long ull;
struct pairhash {
public:
    template<typename T, typename U>
    size_t operator()(const pair<T, U> &x) const {
	size_t seed = hash<T>()(x.first);
	return hash<U>()(x.second) + 0x9e3779b9 + (seed<<6) + (seed>>2);
    }
};

const int inf = 1000000009;
const ll mod = (ll)1e8 + 9;

map<int, int> prime_factor(int n) {
    map<int, int> res;
    for (int i = 2; i * i <= n; i++) {
	while (n % i == 0) {
	    res[i]++;
	    n /= i;
	}
    }
    if (n != 1) res[n] = 1;
    return res;
}

int calc(ll num, int k) {
    int res = 0;
    while (num % k == 0) {
	res++;
	num /= k;
    }
    return res;
}

int solve(ll seed, int n, int k, int b) {
    ll x[n+1];
    x[0] = seed;
    for (int i = 1; i <= n; i++) {
	x[i] = 1 + (x[i-1] * x[i-1] + x[i-1] * 12345) % mod;
    }
    map<int, int> fact = prime_factor(b);

    int res = inf;
    for (auto& p : fact) {
	vector<int> t(n+1);
	for (int i = 0; i <= n; i++) {
	    t[i] = calc(x[i], p.first);
	}
	sort(all(t));
	int sum = 0;
	for (int i = 0; i < k; i++)
	    sum += t[i];
	res = min(res, sum / p.second);
    }
    return res;
}

int main() {
    int q;
    scanf("%d", &q);
    for (int i = 0; i < q; i++) {
	ll seed;
	int n, k, b;
	scanf("%lld%d%d%d", &seed, &n, &k, &b);
	printf("%d\n", solve(seed, n, k, b));
    }
}
0