結果

問題 No.2365 Present of good number
ユーザー a01sa01toa01sa01to
提出日時 2023-06-30 22:10:31
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,337 bytes
コンパイル時間 2,239 ms
コンパイル使用メモリ 214,784 KB
実行使用メモリ 4,500 KB
最終ジャッジ日時 2023-09-21 16:13:01
合計ジャッジ時間 3,786 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 WA -
testcase_03 WA -
testcase_04 AC 1 ms
4,376 KB
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 WA -
testcase_33 WA -
testcase_34 WA -
testcase_35 WA -
testcase_36 WA -
testcase_37 WA -
testcase_38 WA -
testcase_39 WA -
testcase_40 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#ifdef LOCAL
  #include "settings/debug.cpp"
  #define _GLIBCXX_DEBUG
#else
  #define Debug(...) void(0)
#endif
using ll = long long;
#define rep(i, n) for (int i = 0; i < (n); ++i)

#include <atcoder/modint>
using mint = atcoder::modint1000000007;

int main() {
  ll n, k;
  cin >> n >> k;
  vector<pair<ll, ll>> pf(0);
  for (int i = 2; i * i <= n; ++i) {
    if (n % i == 0) {
      pf.push_back({ i, 0 });
      while (n % i == 0) {
        n /= i;
        pf.back().second++;
      }
    }
  }
  if (n > 1) pf.push_back({ n, 1 });
  while (pf.size() > 2 && k > 0) {
    --k;
    map<ll, ll> new_pf;
    new_pf[2] = new_pf[3] = 0;
    for (auto [p, e] : pf) {
      if (p == 2)
        new_pf[2] += e;
      else if (p == 3)
        new_pf[3] += e;
      else {
        assert(p % 2 == 1);
        ll new_p = p + 1;
        while (new_p % 2 == 0) {
          new_p /= 2;
          new_pf[2] += e;
        }
        assert(new_p % 2 == 1);
        if (new_p != 1) {
          new_pf[new_p] += e;
        }
      }
    }
    pf.clear();
    for (auto [p, e] : new_pf) {
      pf.push_back({ p, e });
    }
    Debug(k, pf);
  }
  mint ans = 1;
  for (auto [p, e] : pf) {
    ans *= mint(p + (k % p)).pow((mint(e) * mint(2).pow(k / p)).val());
  }
  cout << ans.val() << endl;
  return 0;
}
0