結果

問題 No.28 末尾最適化
ユーザー 0w10w1
提出日時 2020-09-16 21:53:12
言語 C++17(1z)
(gcc 10.2.0 + boost 1.73.0)
結果
AC  
実行時間 1,568 ms / 5,000 ms
コード長 1,132 Byte
コンパイル時間 2,161 ms
使用メモリ 7,024 KB
最終ジャッジ日時 2021-04-07 09:33:50
合計ジャッジ時間 4,562 ms
ジャッジサーバーID
(参考情報)
judge8 / judge10
このコードへのチャレンジ(β)

テストケース

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

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;

int main() {
  ios::sync_with_stdio(false);

  int Q;
  cin >> Q;

  while (Q--) {
    int Seed, N, K, B;
    cin >> Seed >> N >> K >> B;
    vector<int> X(N + 1); { 
      X[0] = Seed;
      for (int i = 1; i <= N; ++i) {
        X[i] = (1 + ((int64_t) X[i - 1] * X[i - 1] + (int64_t) X[i - 1] * 12345)) % 100000009;
      }
    }
    vector<pair<int, int>> pe; {
      int b = B;
      for (int i = 2; i * i <= b; ++i) {
        if (b % i == 0) {
          pe.emplace_back(0, i);
          while (b % i == 0) pe.back().first++, b /= i;
        }
      }
      if (b > 1) pe.emplace_back(1, b);
    }
    int res = 1 << 30;
    for (auto b : pe) {
      int q = b.second;
      function<int(int, int)> dcnt = [&](int v, int p) { return v % p ? 0 : 1 + dcnt(v / p, p); };
      sort(X.begin(), X.end(), [&](int a, int b) { return dcnt(a, q) < dcnt(b, q); });
      for (auto c : pe) {
        int r = 0;
        for (int i = 0; i < K; ++i) {
          r += dcnt(X[i], c.second);
        }
        res = min(res, r / c.first);
      }
    }
    cout << res << endl;
  }

  return 0;
}

0